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We study the charge transport of the noninteracting electron gas in a two-dimensional 
quantum Hall system with Anderson-type impurities at zero temperature. We prove that 
there exist localized states of the bulk order in the disordered-broadened Landau bands 
whose energies are smaller than a certain value determined by the strength of the uniform 
magnetic field. We also prove that, when the Fermi level lies in the localization regime, 
the Hall conductance is quantized to the desired integer and shows the plateau of the bulk 
order for varying the filling factor of the electrons rather than the Fermi level. 
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1 Introduction 

The two most remarkable facts of the integral quantum Hall effect [1] are the integrality 
of the Hall conductance and its robustness for varying the parameters such as the filling 
factor of the electrons and the strength of the disorder. The integrality is explained by the 
topological nature [2, 3] of the Hall conductance. The constancy of the Hall conductance 
is due to the Anderson localization of the wavefunctions of the electrons [4]. 

First of all we shall survey recent mathematical analysis of the quantum Hall effect. 
As for justification of the conductance formula leading to the topological invariant, sat- 
isfactory results have been obtained in the recent papers within the linear response ap- 
proximation or an adiabatic limit of slowly applying an electric field [5, 6, 7, 8, 9]. Avron, 
Seiler and Yaffe [5] proved that a flux averaged charge transport 1 is quantized to an in- 
teger in the adiabatic limit under the assumption of a nonvanishing spectral gap above a 
non-degenerate ground state for a finite- volume interacting electron gas. In [6], a static 
electric field with a regularized boundary condition was used as an external force to derive 
an electric current for a finite-volume interacting electron gas under the assumption of 
a nonvanishing spectral gap above the sector of the ground state(s). The resulting Hall 
conductance is equal to the universal conductance multiplied by the filling factor of the 
electrons in the infinite- volume limit. When the Fermi level lies in a spectral gap for a 
noninteracting electron gas on the whole plane R 2 , Elgart and Schlein [7] justified the 
Hall conductance formula which is written in terms of switch functions in the adiabatic 
limit. This formula was first introduced by Avron, Seiler and Simon [10]. Without relying 
on the gap assumption, a general conductance formula was obtained for finite-volume in- 
teracting electron gases [8]. For the whole plane, Bouclet, Germinet, Klein and Schenker 
[9] obtained a Hall conductance formula for a random noninteracting electron gas with 
translation ergodicity under the assumption that the Fermi level falls into a localization 
regime. 

As to the localization and the related conductance plateaus, we refer only to a class of 
noninteracting electron gases because the localization of interacting electrons is still an un- 
solved problem. The existence of the localization at the edges of the disordered-broadened 
Landau bands was proved within a single-band approximation [11], for a sufficiently strong 
magnetic field [12, 13, 14], or for a low density of the electrons at the band edges [15]. 
The existence of the quantized Hall conductance plateaus was first proved by Kunz [16] 
under assumptions on a linear response formula of the conductance and on the band edge 
localization. The latter assumption on the localization can be removed for a tight-binding 
model. Namely, the constancy of the quantized Hall conductance was proved within the 
tight-binding approximation for varying the Fermi level [17, 18], or the strength of the 
potential [19]. We should remark that, without relying on the translation ergodicity of 
the Hamiltonian, Elgart, Graf and Schenker [20] proved the constancy of the quantized 
Hall conductance for a tight-binding case. For continuous models, Nakamura and Bellis- 
sard [21] proved that the states at the bottom of the spectrum do not contribute to the 
Hall conductance. Quite recently, Germinet, Klein and Schenker [22] proved that the Hall 
conductance formula [10] which is written in terms of switch functions shows a plateau 

1 This is a non-trivial charge transport which is intrinsically different from the response to a static 
external field. 
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for a random Landau Hamiltonian with translation ergodicity. In order to determine the 
integer of the quantized value of the Hall conductance, they further required the condition 
that the disordered-broadened Landau bands are disjoint, i.e., there exists a nonvanish- 
ing spectral gap between two neighboring Landau bands. However, the existence of the 
localized states at the band edges does not necessarily implies the appearance of the Hall 
conductance plateaus for varying the filling factor of the electrons because the density of 
the localized states may be vanishing in the infinite volume. In order to show the existence 
of such plateaus, we need to prove the existence of localized states of the bulk order. In 
passing, we remark that Wang [23] obtained the asymptotic expansion for the density of 
states in the large magnetic field limit. 2 

In this paper, we focus on the issue of proving the bulk order plateaus, and consider 
a noninteracting electron gas with Anderson-type impurities in a magnetic field in two 
dimensions at zero temperature. The centers of the bumps of the impurities form the 
triangular lattice. First we prove that there exist localized states of the bulk order in 
the disordered-broadened Landau bands whose energies are smaller than a certain value 
determined by the strength of the magnetic field. In order to obtain the Hall conductance 
as a linear response coefficient to an external electric field, we apply a time-dependent 
vector potential A ex (t) = (0, cx(t)), where the function a(t) of time t is given by (2.18) in 
the next section. For t G [— T, 0] with a large positive T, the corresponding electric field is 
adiabatically switched on, and for t > 0, the electric field becomes (0, F) with the constant 
strength F. First we consider the finite, isolated system of an L x x L y rectangular box, and 
impose periodic boundary conditions for the wavefunctions with the help of the magnetic 
translation (2.17), and then we take the infinite- volume limit. The explicit expression of 
the conductance formula which we will use is given in ref. [8]. We should remark that, 
when the system is translationally invariant, the constant Hall current flows on the torus 
without dissipation of energy as in ref. [8]. We prove that, when the Fermi level lies in the 
localization regime, the Hall conductance is quantized to the desired integer and shows the 
plateau of the bulk order for varying the filling factor of the electrons. In our approach, we 
require neither the disjoint condition for the Landau bands nor the translation ergodicity 3 
of the Hamiltonian which were assumed in [22] as mentioned above. Instead of these 
condition, we need the "covering condition" that the whole plane R 2 is covered by the 
supports of the bumps of the impurity potentials so that the sum of the bumps is strictly 
positive on the whole plane R 2 . This "covering condition" is not required in [22]. 

The present paper is organized as follows. In Section 2, we describe the model, and 
state our main theorems. As preliminaries, we study the spectrum of the Hamiltonian 
without the random potential in Section 3 and the site percolation on the triangular 
lattice of the impurities in Section 4. In Section 5, we obtain a decay bound for the 
resolvent (Green function) of a finite volume. This bound becomes the initial data for the 

2 In general, an asymptotic series does not give us any information for a fixed finite value of the 
parameter because the asymptotic series is not necessarily convergent. See, for example, Section XII. 3 of 
the book [24] . Thus the result of [23] dose not imply the existence of localized states of the bulk order 
for a fixed finite value of the magnetic field. See also the recent paper [25] for the difficulty of obtaining 
a lower bound for the density of states. 

3 In a generic, realistic situation that there exist one- or two-dimensional objects such as dislocations 
in crystals and interfaces in semiconductors, we cannot expect that the system has translation ergodicity. 
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multi-scale analysis [26, 27, 28] which is given in Section 6. In order to prove constancy of 
the Hall conductance, we further need the fractional moment bound [29] for the resolvent. 
The bound is given in Section 7. As preliminaries for proving the integrality and constancy 
of the Hall conductance, we study the finite volume Hall conductance in Section 8. The 
integrality of the Hall conductance is proved within the framework of "noncommutative 
geometry" [10, 17, 18] in Section 9, and the constancy is proved by using the homotopy 
argument [17, 19] in Section 10. The widths of the Hall conductance plateaus and the 
corrections to the linear response formula are estimated in Sections 11 and 12, respectively. 
Appendices A-H are devoted to technical estimates. The standard Hall conductance 
formula which is given in Section 9 is written in terms of the position operator of the 
electron. In Appendix I, we give a proof that this Hall conductance is equal to another 
Hall conductance [10] which is written in terms of switch functions for a class of continuous 
models. 



2 Model and the main results 

Consider a two-dimensional electron system with Anderson-type impurities in a uniform 
magnetic field (0, 0, B) perpendicular to the x-y plane in which the electron is confined. 
For simplicity we assume that the electron does not have the spin degrees of freedom. The 
Hamiltonian is given by 

H U = H + V U (2.1) 

with the unperturbed Hamiltonian, 

H = -^-(p + eA) 2 + V , (2.2) 
2m e 

and with a random potential V u , where p := —ihV with the Planck constant H, and — e 
and m e are, respectively, the charge of electron and the mass of electron; A and Vo are, 
respectively, a vector potential and an electrostatic potential. The system is defined on a 
rectangular box 

A sys := [-LJ2, L x /2] x [-L y /2, L y /2] C R 2 (2.3) 

with the periodic boundary conditions. The vector potential A = (A x , A y ) consists of two 
parts as A = Ap + A , where A (r) = {—By, 0) which gives the uniform magnetic field 
and the vector potential A P satisfies the periodic boundary condition, 

A P (x + L x , y) = A P (x, y + L y ) = A P (x, y). (2.4) 

This condition for A P implies that the corresponding magnetic flux piercing the rectan- 
gular box A sys is vanishing. Therefore the total magnetic flux is given by BL x L y from the 
vector potential A only. We assume that the components of the vector potential Ap are 
continuously differentiable on R 2 . Further we assume that the electrostatic potential Vo 
satisfies the periodic boundary condition, 



V (x + L x , y) = V (x, y + L y ) = V (x, y), 



(2.5) 
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Figure 1: The parallelogram in the triangular lattice with its dual, hexagonal lattice. 



and 1 1 1 1 oo + 1 1 IV | |oo < Vq < oc with some positive constant v which is independent of 
the system sizes L x , L y . Here Vq 1 = max{±Vo, 0}. As a random potential V^, we consider 
an Anderson-type impurity potential, 

V u (r) = £ A,H«(r-z), (2.6) 

zGL 2 

for r := (x,y) G R 2 . The constants {A z (cj)| z G L 2 } form a family of independent, 
identically distributed random variables on the two-dimensional triangular lattice L 2 C R 2 
with the lattice constant a > 0. The common distribution of the random variables has a 
density g > which has compact support, i.e., supp g C [A min , A max ] with A min < < A max . 
Further the density g satisfies the following conditions: 

g G L°°(R) n C(R) and f + g(X)d\ > 1/2 (2.7) 

with two positive numbers A + and A_. We consider two cases: (i) a small A_ and (ii) 
a small A+. We assume that the condition (2.7) holds for both of the two cases. If the 
density g is an even function of A and is concentrated near A = 0, this requirement holds. 
We take 

L 2 ={z=ma 1 +na 2 |(m,n)GZ 2 } (2.8) 

with the two primitive translation vectors, a x = (a, 0) and a 2 = (a/2, y/3a/2 s j. The trian- 
gular lattice is imbedded in R 2 such that each face is an equilateral triangle as described 
in Fig. 1. We also consider its dual, hexagonal lattice which is defined as follows. Choose 
a vertex of the dual lattice at the center of gravity of each triangle, i.e., the intersection 
of the bisectors of the sides of the triangle. For the edges of the dual lattice, take the 
line segments along these same bisectors, and connecting the centers of gravity of adjacent 
triangles. 
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We assume that the bump u of the single-site potential in (2.6) satisfies the following 
conditions: < u G L°°(R 2 ), 

u(r) = for |r| > r u with a constant r u G (VSa/3,VSa/2j , (2.9) 

and 

u(r) > u > for r in the face of the hexagon with the center r = of gravity. (2.10) 

Here uq is a positive constant. The first condition (2.9) implies that the single-site po- 
tentials u has compact support, and overlap with only nearest neighbor u. The next 
condition (2.10) implies that the whole space R 2 is covered by the supports of the bumps 
{u{- — z)} zeL 2 of the impurities so that 

^2 u ( r - z ) > u o for any r G R 2 . (2-11) 
zeL 2 

We should remark that this "covering condition" is needed for estimating the number of 
the localized states and for applying the fractional moment method [29], which yields a 
decaying bound for the resolvent. 

Clearly the random potential V u of (2.6) does not necessarily satisfy the periodic bound- 
ary condition, 

V u (x + L x , y) = V u (x, y + L y ) = V u (x, y), (2.12) 

without a special relation between the lattice constant a and the system sizes L x ,L y . 
Therefore we will replace the random potential V u with which is slightly different from 
V w in a neighborhood of the boundaries so that V u satisfies the periodic boundary condition 
(2.12). Before proceeding further, we check that the boundary effect due to this procedure 
is almost negligible and does not affect the following argument. Write 

Ll/2 = N x a and Ly/2 = N y ■ V3a/2 with positive integers N x , N y . (2.13) 

When we take the sizes to be L x = I? x and L y = L y , the periodic boundary condition 
is automatically satisfied without replacing the random potential. However, for a given 
lattice constant a, the sizes do not necessarily satisfy the flux quantization condition, 
L^Ly = 2irM£ 2 B , which we need in the following argument. Here M is a positive integer, 

and £ B is the so-called magnetic length defined as £b '■= \Jh/(eB). In a generic situation, 
we have 

2nM£% < L^Ly < 2vr(M + l)f B with some positive integer M. (2.14) 

In order to recover the flux quantization condition, we change the size a little bit in the y 
direction. Namely we choose the sizes as L x = I? x and L y = L y — SL y with a small 5L y . 
Substituting this into L x L y = 2irMl 2 B , we have 

< SL y < 2n£ 2 B /L x . (2.15) 

Notice L y < Ly, and consider the triangles which overlap with the upper boundary of A sys . 
We replace these equilateral triangles with the isosceles triangles of the height (\/3a/2 — 
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SL y ). From the above bound (2.15) for SL y , the height of the isosceles triangles is slightly 
shorter than the height y/3a/2 of the equilateral triangles for a sufficiently large L x . In 
the same way as in (2.6), we put the impurity potentials on the center of the gravity of 
the isosceles triangles. Because of the bound (2.15) for 5L y , the effect of this procedure 
is negligibly small for a sufficiently large size L x . When the boundary effect plays an 
essential role as in Section 9 below, we denote by V Ui \ this resulting random potential for 
a region A. Otherwise, we will often use the same notation V u for short. 

When Ap = 0, we require the differentiability, V — Vq + V w e C 2 , in addition to 
the above conditions, in order to obtain the exponential decay bound for the resolvent 
(H w — z)^ 1 in Appendix D.2. 

As mentioned above, we require the flux quantization condition, L x L y = 2nMl 2 Bl with 
a sufficiently large positive integer M. The number M is exactly equal to the number 
of the states in a single Landau level of the single-electron Hamiltonian in the simple 
uniform magnetic field with no electrostatic potential. This condition L x L y = 2irM£ 2 B for 
the sizes L x , L y is convenient for imposing the following periodic boundary conditions: For 
an electron wavefunction <p, we impose periodic boundary conditions, 

&\L x )<p(r) = <p(r) and t^(L y )ip(r) = ip(r), (2.16) 

where t^ x \- ■ •) and t^(- ■ •) are magnetic translation operators [30] defined as 

&\x')f(x,y) = f(x-x',y), t^(y')f(x,y) = ex P [ty'x/£ 2 B }f(x,y-y') (2.17) 

for a function / on R 2 . 

In order to measure the conductance, we introduce the time-dependent vector field 
A cx (t) = (0,a(t)) with 

«W=-«x{^ (2.18) 

Here F is the strength of the electric field, and rj > is a small adiabatic parameter. The 
y-component of the corresponding external electric field is given by 

The time-dependent Hamiltonian is given by 

H u (t) = -^[p + eA + eA cx (t)] 2 + V + V u . (2.20) 
The velocity operator is 

v(f) = (v x (t),v y (t)) = — [p + eA + eA ex (0]. (2.21) 

Tfl e 

Let U(t,to) be the time evolution operator. We choose the initial time t = to — —T with 
a large T > 0. Then the total current density is given by 

jtot(t) = (jtot,o ; (t),Jtot, s/ (t)) = - T ^Tr[/ t (t,to)v(t)[/(t,to)PF for t>0, (2.22) 
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where Pp is the projection on energies smaller than the Fermi energy E F . This total current 
density is decomposed into [8] the initial current density jo and the induced current density 
jind(0 as j t ot(0 = jo +jind(*)- Here the initial current density j is given by 

jo = -T4- Tr — [P + eA ] P F- ( 2 - 23 ) 
±j x L y m e 

Further the induced current density jind(^) is decomposed into the linear part and the 
nonlinear part in the strength F as 

jind(t) = (^tot^(t), a tot , yy (t))F + j^(t) with j/ d (0 = o(F), (2.24) 

where the coefficients, a to t,sy{t) for s — x,y, of the linear term are the total conductance 
which are written 

(t) = o sy + 6a sy (t) (2.25) 

with the small corrections, Sa sy (t), due to the initial adiabatic process, and o(F) denotes 
a quantity q satisfying q/F — > as F — + 0. Since the present system has no electron- 
electron interaction, the order estimate for the nonlinear part y^it) in (2.24) holds also 
in the infinite volume limit [8, 9, 31] with the same form of the linear part of the induced 
current. But the nonlinear part jind(0 depends on the adiabatic parameters 77 and T. 
Therefore we cannot take the adiabatic limit T | 00 and 77 J, for the nonlinear part 
JindW °f the induced current. 

Now we describe our main theorems. Let v = N/ M be the filling factor of the electrons 
for a finite volume, where N is the number of the electrons, and write u c = eB/m e for the 
cyclotron frequency. Consider first the case with Ap = in the infinite volume limit. 

Theorem 2.1 Assume that the filling factor v satisfies n — 1 < v < n with a positive 
integer n. Then there exist positive constants, B (n) andv , such that there appear localized 
states of the bulk order around the energy £ n -i — (n — l/2)hu c , i.e., the n-th Landau 
band center, for any magnetic field B > B (n) and for any potential Vo E C 2 satisfying 
||Vo + ||oo + || II 00 < vq. Further, when the Fermi level lies in the localization regime, the 
conductances are quantized as 

e 2 f n for the upper localization regime 



0~ X y — ^ x 



r n for the upper localization regime _ ,^ 2Q) 

\ (n — 1) for the lower localization regime yy ' 



and exhibit the plateaus for varying the filling factor. With probability one, there exist 
positive constants, Cj{uS) < 00, j = 1,2,3, such that the corrections 5a sy (t) due to the 
initial adiabatic process satisfy 

\Sa sy (t)\ < [CiH + C 2 {uj)T]e^ T + C 3 (^ 1/13 , (2.27) 

and that the expectation E[Cj] of the positive constants Cj(uj) is finite for j = 1,2,3, i.e., 
E[Cj] < 00. Here the constant Cj(cu) itself without the expectation may depend on the 
random event uj of the random potential V u . 
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In the case with A P 7^ in the infinite volume limit, we require a strong disorder 
because of a technical reason. (See Appendix D for details.) We take u = Tiuj c u with 
a fixed, dimensionless function u for the random potential V u of (2.6). This potential 
behaves as ||m||oo ~ Const, x B for a large B. For this random potential, we have: 

Theorem 2.2 Assume that the filling factor v satisfies n — 1 < v < n. Then there 
exist positive constants, B Q (n),ao(n) and wq, such that there appear localized states of 
the bulk order around the energy £ n -i — (n — l/2)hu c for any magnetic field B > B (n) 
and for any vector potential A P satisfying HIAplH^ < a^njB 1 / 2 and for the function u 
satisfying ||w||oo < w o- Further, when the Fermi level lies in the localization regime, the 
conductances, a xy and a yy , are quantized as in (2.26) and exhibit the plateaus for varying 
the filling factor. With probability one, the corrections Sa sy (t) due to the initial adiabatic 
process satisfy a bound having the same form as that of the bound (2.27). 

Remark: 1. In the conditions in the first theorem, we can also take u = hu c u for the 
random potential with a small norm ||w||oo- 

2. In the conditions in the second theorem, we can take V which behaves as ||V ||oo ~ 
Const, x B for a large B with a small positive constant, instead of a fixed potential V . 

3. These two theorems do not necessarily state that both of the upper and lower band 
edges exhibit the localized states of the bulk order. Namely a localization region of the 
bulk order may appears only at one side for a single Landau band. 

4. We do not require the assumption that the n-th Landau band is separated from the 
rest of the spectrum by two spectral gaps. 

5. In the previous analyses [12, 13, 14, 22], they considered the Hamiltonian having 
the form with Ap = and Vq = 0. The analyses rely on the special properties of the 
unperturbed Hamiltonian. For example, they use the explicit forms of the integral kernel 
of the projections onto the Landau levels. The extension to the case with Ap 7^ and 
Vb 7^ needs additional, non-trivial analyses for localization. In addition, we do not require 
a periodicity of the potentials A P , V with a finite period. Therefore the Hamiltonian 
does not need to be translation ergodic. 

6. The widths of the plateaus can be estimated as we will show in Section 11 below. In 
particular, when Ap = and Vq = 0, the ratio of the localized states to the total number 
M of the states in the single Landau level tends to one as the strength B of the magnetic 
field goes to infinity. This implies that our estimate for the widths of the plateaus shows 
the optimal, expected value in this limit. 



In order to show the localization for the disordered-broadened Landau bands, we first 
need to check the condition for the appearance of the spectral gaps in the spectrum of the 
Hamiltonian H of (2.2) with a generic, bounded potential V$. 

First we recall the simplest Landau Hamiltonian for a single electron only in the uniform 
magnetic field. The Hamiltonian is given by 



3 Spectral gaps of the Hamiltonian Hq 




(3.1) 
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We assume that the electron is confined to the same finite rectangular box A sys of (2.3) 
as the box for the Hamiltonian H w of (2.1), and impose the periodic boundary conditions 
(2.16) for the wavefunctions with the flux quantization condition L x L y = 2ttM£ 2 b . Then 
the energy eigenvalues £ n of are given by 4 



S n := (^n + %uj c for n — 0, 1, 2, 

The Hamiltonian H of (2.2) on the finite box A sys is written 
1 



(3.2) 



Hn 



2m f 



(p + eA + eA P ) 2 + V 



= ^ + 2^ Ap ' (P + eAo) 
Using the Schwarz inequality one has 



2m f 



(p + eA ) • A P 



2m P 



\A P \ 2 + V . (3.3) 



| (^, A P • (p + eA )^) | < || | A P | IL Jty, (p + eA ) 2 ^) 



(3.4) 



for the normalized vector ip in the domain of the Hamiltonian. From this inequality the 
energy expectation can be evaluated as 



V2e 



2m P 



+ 



Vn 



and 



V2e 



'm P 



Vn 



(3.5) 



(3.6) 



where = max{±Vo,0}. Let us denote by £^ d + e and S^ d - e , respectively, the upper and 
lower edges of the n+ 1-th Landau band which is broadened by the potentials Vq and Ap. 
From the standard argument about the min-max principle, 5 one has 



ocdgc < c , 



\[2t 



■"n,+ 



I At 



- v Sn+ 2m P 



IA T 



+ 



(3.7) 



for n — 0, 1, 2, For the lower edge, we assume 



y / 2m^ 



< 



-hu r . 



(3.8) 



Then the right-hand side of the bound (3.6) is a strictly monotone increasing function of 
the expectation (ip,H L ip). Therefore, the same argument yields 



,7 edge ^ p 



V2e 



(3.9) 



4 See, for example, Refs. [6, 32]. 

5 See, for example, Section XIII. 1 of the book [24] by M. Reed and B. Simon. 
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for n = 0,1,2, If this right-hand side with the index n + 1 is strictly larger than 

the right-hand side of (3.7) with the index n, then there exists a spectral gap above the 
Landau band with the index n, i.e., 

ClT,- > £n d f- This gap condition can be written as 



flUJ c > 



V2e 



(v^+V^) + ^-|||A r 



+ 



v: 



+ 



(3.10) 



Clearly this is stronger than the condition (3.8) for the vector potential Ap. Therefore we 
have no need to take into account the condition (3.8). 



4 Site percolation on the triangular lattice 

The classical motion of the electron is forbidden in the regions that the strength of the 
potential is smaller than the deviation of the energy of the electron from the Landau ener- 
gies S n of (3.2). In those regions, the Green function of the electron decays exponentially. 
In order to get the decay bound for the Green function, we study the distribution of those 
classically forbidden regions. We reformulate this problem as a site percolation problem 
on the triangular lattice. The idea of using percolation is due to Combes and Hislop [12] 
or Wang [13]. But both of their random potentials are different from the present potential 
which we require for estimating the number of the localized states. 

We begin with setting up site percolation on the triangular lattice L 2 for the present 
random potential. We say that the site z 6 L 2 is occupied if \ z (lu) G (— A_,A+). The 
probability p that a site z is occupied is given by 

P 



I + g(X)d\. (4.1) 

J — A 



The assumption (2.7) implies p > p c = 1/2. Here p c is the critical probability which 
equals 1/2 for the present site percolation on the triangular lattice [33, 34]. A path of 
L 2 is a sequence z , z 1; . . . , z n of sites Zj such that all of the adjacent two site Zj, Zj +1 are 
corresponding to a side of a unit triangle. If z = z„, then we say that the path is closed, 
and we call a closed path a circuit. If all of the site Zj of the path are occupied, then 
we say that the path is occupied. Similarly we define an unoccupied path, an occupied 
circuit, etc. We denote by P p (A) the probability that an event A occurs. 

Let lie/' be a parallelogram with the lengths £a, £'a of the sides in the triangular lattice. 
See Fig. I. More precisely, it is given by 



n^/ := I mai + hsl 2 



t2 



\m\<-,\n\<-,(m,n)eZ 2 . (4.2) 



Here we take £, £' even integers for simplicity. Consider the event A^p that there exists a 
unoccupied path in the parallelogram H^v joining a site on the lower side with the length 
t to a site on the upper side. Since the connectivity between two sites with an unoccupied 
path decays exponentially for p > p c , the probability P P {A^') that the event A^p occurs 
is bounded as 

P P (A~e,e) < Const, x £exp[-m p £'} for p > p c = 1/2, (4.3) 
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where m p is a positive function of p. Let Apj be the event that there exists an occupied 
path in the parallelogram IL^/ joining a site on the left side with the length £' to a site 
on the right side. Then this event Ay^ is the complementary event of A^> because of the 
structure of the triangular lattice. Immediately, 

P p (A e/ )+P p (A e/ ,) = l. (4.4) 

Combining this with the above inequality, one has 

P p (A e >,e) > 1 - Const, x £exp[-m p £'} for p > p c = 1/2. (4.5) 

Consider a parallelogram-shaped region consisting of hexagons such that 

A^, +1 (z ) := U h z+zo , (4.6) 

where h z is the region of the face of the hexagon with the center z, including the six sides 
of the hexagon, and z is the center of the region A^ r ] a ^, +1 (z ). The region has the jagged 
boundary as in Fig. 1. Moreover we define an annular region as 

A^W := A§S,(zo)\A£>o), (4.7) 

where both £ and £' take an odd integer. 

Let us consider the event in A^3",(z ) flL 2 that there exists an occupied circuit C 
encircling the inside region A^ a (z ). Then, from the inequality (4.5) and FKG inequality, 6 
the probability P P (D^>) for this event satisfies 

P P (D et e) > [P P {A VM )f [Pp(Ae,3t>)] 2 

> 1 - Const, x {£ exp[-m p £'} + £' exp[-m p £]} 

for p>p c = 1/2. (4.8) 

Therefore the event of an occupied circuit occurs with the probability nearly equal to one 
for large £, £'. 

We denote by bjj+i the side zjzj^j" of a unit triangle, i.e., 

bjj+i ■= {r = Az, + (1 - A)z i+1 | A e [0, 1]}. (4.9) 
We define the region TZjj+i including the side bjj+i as 

Kj, j+ i := {r|dist(r, b jJ+1 ) < rj with r l := —a - r u , (4.10) 

where r u is given in the condition (2.9) for the bump u. Further we define the ribbon 
region IZc associated a circuit C by 

^c:= (J *W 

Clearly r\ is strictly positive from the condition r u E (y/3a/3, \/3a/2) of (2.9), and the 
ribbon region IZc has a nonzero width 2r x . 

6 See, for example, the book [34]. 
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Proposition 4.1 There appears an occupied circuit C in the annular region A|^,(z ) 
with a probability larger than 

ppcrc ;= x _ ^pcrc {£ exp [_ m ^'] + g e x P [-m p £]} , (4.12) 

where C pcrc is the positive constant in the right-hand side of the above bound (4-8), and 
p > p c = 1/2 is given by (4-1)- Further the following hold for the ribbon region IZc 
associated with the occupied circuit C: 

dist (72c, 0A£$,(z o )) =r u -^j-a=:r 2 >0 (4.13) 

and 

-A_«i < K,(r) < A+mi for r e TZ C . (4.14) 
Here dA^^/^o) is the boundary of the annular region A 3 " 3 £/(z ), and u\ : = 211-ujloo. 

Proof: The lower bound (4.12) of the probability is nothing but the right-hand side 
of (4.8). The positivity (4.13) of r 2 follows from the condition (2.9) of the bump u of 
the random potential V^, and the bound (4.14) follows from the condition (2.9) and the 
definition (2.6) of the random potential V u . ■ 



5 Initial decay estimate for the resolvent 

Now let us estimate the decay of the resolvent (Green function) for a finite parallelogram- 
shaped region. The resulting decay bound in Proposition 5.2 below will become the initial 
data for the multi-scale analysis to obtain the decay bounds for the resolvent in larger scales 
in the next section. However, by using the multi-scale analysis, we cannot get a similar 
decay bound for the resolvent for two arbitrary points in the infinite-volume limit. On 
the other hand, the fractional moment method leads us to a decay bound for a fractional 
moment of the resolvent in the infinite- volume limit. Actually, as we will see in Section 7, 
the initial decay estimate of this section yields such a decay bound. But the resolvent 
itself without taking a fractional moment cannot be evaluated by the method [29]. Due 
to technical reason related to these observations, we need both multi-scale analysis and 
fractional moment analysis, in order to prove the existence of the conductance plateaus 
with a bulk order width. 

Although the method in this section is basically the same as in the previous papers 
[12, 13] as mentioned in the preceding section, we need more detailed analysis about the 
magnetic field dependence of the decay bounds, in order to estimate the number of the 
localized states which yield the Hall conductance plateau with a bulk order width. 

Fix the random variables A z+Z0 with z e L 2 \n 3 ^_i )3 ^_ 1 . Here £ and £' are odd integers 
larger than 1. Consider the parallelogram-shaped region A^ r 3 ^,(z ) centered at z , and 
assume A^^^zq) C A sys with a sufficiently large box A sys of (2.3). We write A 3 ^/ = 
A^3^,(z ) for short. Further we consider the Hamiltonian H w restricting to the region 
A 3 £ t3 £i with the Dirichlet boundary conditions. The Hamiltonian is written as 

H **u* = 2^"(P + eA ) 2 + V&Ia*,*, + VW* + SV UJ3e ,ae, (5-1) 
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where we have decomposed the random potential V w into two parts, 



V u ,3e,ae(r)= A z+Z0 (cj)w(r - z - z). 

zen 



(5.2) 



and 5V Ut3e *e = Vj 



3l-l,Zl'-l 

Vuj,3£,3£'- Clearly, the first part V w $t$t> of the random potential 



is determined by only the random variables A z+Zo (o;) at the sites z + zo lying in the 
parallelogram-shaped region A|^,(z ), and so it is independent of the outside random 
variables. Further we have 



Y u ( r - z - z) > m for any r G A 3 £,3£' 

zen 3£-l,3£'-l 



(5.3) 



from the assumption (2.10). This condition will be useful to obtain the Wegner estimate 
[35] for the density of the states. (See Appendix A for the Wegner estimate.) 

On the other hand, the random potential SV^^^e which is supported by only the region 
near the boundary of A 3 ^ i3 ^, depends on the the outside random variables. Following 
[12], we absorb this term into the operator W(x) of (5.7) below which will appear in 
the geometric resolvent equation. (See Appendix F for details.) Thus we consider the 
Hamiltonian, 

1 



jj^ 



(5.4) 



without the potential SV^^e^e, instead of the Hamiltonian H\ , of (5.1). 
Assume that the energy E G R satisfies the condition, 

£n d + + < E < S^i _ — \-Ui with u\ = 2 lliiHoQ. 



H> 



31,31' 



E-ie 



(5.5) 
with e G R. 



We write the resolvent as R^w = Rse,3e'(E + iz) '- 
For 5 G (0,r 2 ), consider the region, A^^, := {r G A 3 ^/| dist(r, dAse,3e) > where r 2 
is given by (4.13), and dh.u,u' is the boundary of the region A 3 ^ 3 £/. Let xle,w De a C 2 , 
positive cut-off function which satisfies 



= 1 and supp Vxit,3t> 

31,31' 



C A 3 ^ 3 ^\A^ i3 ^. 



(5.6) 



We denote by xe,e> the characteristic function of the region Aj^ a (z ). 

The purpose of this section is to estimate the decay of W(x3£ 3£')R3£,3e'Xi,t'i where 



W( X ) 

for a C 2 function x- Note that 



(p + eA) 2 /(2m e ),x 



(5.7) 



,31'XLt' 



< 



h 2 



2m P 



\{^xitM')Ru,u'Xt,t' 



+ — \\( d iXli,3£>)(Pi + eA i) R 3e,3e>X£,e> 



m 



e i=x,y 



(5.8) 
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where we have written V = (d x , d y ), and used 



ih 



ti 2 



ih 



W( X ) = (p + eA) • Vx + 7T- A X = Vx • (p + eA) 



m P 



2m P 



m P 



2m, 



(5.9) 



We write 1Z for the ribbon region TZc in Proposition 4.1 for short. Let e be a small 
positive number, and let C e := {r G 7?.|dist(r, C) < e/2}, so that the region C e has the 
width e, where C is the occupied, closed path in Proposition 4.1. Let Xi be a C 2 , positive 
cut-off function satisfying Xi A ee , =1 an d supp|Vxi| C C t , where we have written A^/ 
for A^ a (z ) for short. Since we can take Xi to satisfy (dixle,3i')Xi = from the definition 
of xleje', one has 

( d iXit,3t')(Pi + eA i) R u,wXi,i' = D i R u,M'XiXw = - D i R u,u'W{xi)Ru,wXt,i'^ ( 5 - 10 ) 

where we have written Di = (c^XiL 3^')(Pi + eAi). In order to estimate this right-hand 
side, we define the e border of the ribbon region 1Z by 7Z e := {r e 72-|dist(r, cflZ) < e}, 
and define r 3 := dist(7?. e , C e ) > 0. We choose a small parameter e so that the distance 
r 3 becomes strictly positive. Further we introduce two C 2 , positive cut-off functions, Xn 
and Xn-i which satisfy the following conditions: 



~e/2 

Xn 



Tz\n e 



/2 



1, supp Vxn 



C 11, 



e/2, 



and 

Xn\n\n t = 
Consider the Hamiltonian, 



Xk\n e/a = °> and su PP|Vx^l C n t \n t/2 . 



H 



K 



2m P 



eA) 2 + ^, 



(5.11) 
(5.12) 

(5.13) 



on the finite rectangular box A sys of (2.3), where we impose the periodic boundary con- 
ditions (2.16) with the flux quantization condition L x L y = 2nM£ 2 B , and the poten- 
tial is given by Vr = X-r^Y® + K>)- Then one has the geometric resolvent equation, 
R3e,3£'Xn = Xn R Ti ~ R3i,3i'W{xii) R Ti with R n := (H n - E - ie)' 1 , where we have used 
XnXn = Xn which is easily obtained from the definitions. Using this equation and 
DiXii = 0; the right-hand side of (5.10) is written as 

—D i R 3 i i se , W(xi)R3e,3e'Xt,e' — ~ DiRu,3i'XvW{Xi)R3i,3i'Xt,t' 

= DiR u ^W(xn)Rn w i.Xi)R3i,3i'Xt,i'- (5.14) 

Consequently, one obtains 

DiRu&rXw = D i R 3e,3e'W(xn) R Ti w (Xi)R3e,3e'Xe,e'- (5.15) 

In the same way, 

(Axle,3e') R 3e,3e'Xe,e' = { A XuM') R ^M' w {Xn) R nW{xi) R 3t,3t'Xi,e- (5-16) 



Widths of the Hall Conductance Plateaus 



17 



If the random potential V u satisfies the condition (4.14), then the energy E satisfying 
the condition (5.5) is in the spectral gap of the Hamiltonian H-n- Therefore we can apply 
the Combes-Thomas method [36] to evaluate decay of the resolvent R-jz in the right-hand 
sides of (5.15) and (5.16). We write A = supp|Vx^| and B = C e , and denote by \A an d 
Xb the characteristic function of A and B, respectively. The resolvent R-ji decays as 



IXaR-kXbW < C{ 'e 



(") c -/3r 3 



(5.17) 



with a probability larger than pp crc of (4.12) in Proposition 4.1, where C[ n ^ and (3 are 
positive constants. The derivation of this decay bound and the explicit parameter depen- 
dence of C[ n ^ and of (3 for the present model are given in Appendix D. In the application 
of the decay bound, we choose the cut-off function Xn so that the region supp|Vx^.| has 
a smooth boundary. 

Lemma 5.1 Let R = (H UJ — E—ie)~ 1 with generic, bounded potentials Vo,V u andE,eEH 
satisfying E ai^Hj) or e ^ ; and let a = (a x , a y ) be a vector-valued C l function. Then 



\a ■ (p + eA)i*|| < 2V2hT e \\R\\ 1 / 2 {l + / E ,*) 1/2 max^HU}, 



i=x,y 



\\(pi + eAi)R(p + eA) • a\\ < 2m e \\ |a| ||oo(l + Se,r) 



and 



\\R(p + eA) • a|| < V^h\\ \a\ || oo || 1/2 ( 1 + f E , R ) 1/2 , 
where we have written 



fE*=[\E\ + \\(V - + V ll 



w 7l|ooJ 



\R\ 



(5.18) 
(5.19) 

(5.20) 

(5.21) 



The proof is given in Appendix E. From (5.9) and these bounds of Lemma 5.1, one has 

\\W{ X i)RuM\ < h(\E\,\\RuM\), ( 5 - 22 ) 

\\Ru,wW{xn)\\ < f2(\E\, WRseMl) (5-23) 

and 

\\(Pi + eA)Rse,se'W(x £ n )\\ < f 3 (\E\,\\R 3e , 3e ,\\), (5.24) 

for the operators in the right-hand sides of (5.15) and (5.16), where the functions, /i,/2 
and f 3 , are given by 

-2 



h(\E\,\\R\ 



2m, 



l|Ax 



1 oo 



\R\ 




+ 2hJ—{\\R\\ + [\E\ + \\(V - + KT)lloo] ||i?|| 2 } 1/2 max{||^ X i||oo}, 



(5.25) 



MINIMI) = 



2m, 



■IIA^IIooPII 



+ hJ^- {\\R\\ + [\e\ + \\(v - + KT)lloo] l|i?ll 2 } V2 |||Vx 



VR\ Woo 



(5.26) 
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and 



h(\E\A\R\\) = ^IIA^iufii^ + ^i + iK^ + K;)^]^!! 2 } 
+ 2ft|||v x y|oc {1 + [\e\ + \\(Yq + KDIloo] l|i?ll} . 



(5.27) 



The norm ||-R3^'|| of the resolvent in these upper bounds can be evaluated by using 
the Wegner estimate. See Appendix A for details. From the resulting Theorem A. 2, we 
have that, for any 5E > 0, 

\\Rat,3t'\\ < W 1 ( 5 - 28 ) 

with a probability larger than 



1 - CyjKzWgW^SE 



(5.29) 



where C\y is a positive constant, and the positive constant K 3 is given by (A. 39) in the 
theorem. 

Proposition 5.2 For any E satisfying the gap condition (5.5), and for any SE > 0, the 
following bound is valid: 



sup 

£^0 



W(xteMeME + ie) W < Cf : ] e~^ h(\E\, (SB)' 1 



x 



1,2 -\\A X 5 M , 3e \U2(\E\, (SE)' 1 ) + — max \\d iX ie, 3e \Us(\E\, (SE)' 1 ) 



(5.30) 



_2m e 

with probability at least 
P hli : = 1 - {C pcrc {£exp[-m/] + £'exp[-m/]} + C^K^gW^E |Aj££,(z )|} . (5.31) 

Proof: Combining (5.15), (5.16), (5.17), (5.22), (5.23), (5.24) and (5.28), the right-hand 
side of (5.8) is estimated. Using (4.12), (5.29) and the inequality Prob(An5) > Prob(A) + 
Prob(-B) — 1, the probability is estimated. ■ 

Since we can take the ribbon region 1Z satisfying supp SV^^t^e H 1Z C dlZ, we can 
obtain a similar bound for \\8V^^£,wRu,w{E + ie)xe,e'\\ to (5.30). Fix the ratio £'/£. For 
simplicity we take £' = £. Fix £ > 4. We choose £ = £ to satisfy 



2C perc £ exp[-m p £ ] < £ ^/2, 
and choose 5E in (5.28) so that 

C w K 3 \\g\UE A^ () (z ) =tf/2. 



(5.32) 



(5.33) 



Clearly these implies P im > 1 — £ ^. Therefore, if we can find a large (3 in (5.30) so that 
the right-hand side of (5.30) with £ = £' = £ becomes small, then we have 



Prob 



sup 



(E + ie) X e , 



< e -70^0 



>i-4 ? 



(5.34) 
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with some 7o > 0, where := W{xlt,u) - 5v ^,u,^Xu,u- 



For the norm 



XwRm,w{E + iz)W \xu,w) i one can obtain a similar bound to (5.30) 



in the same way. As a result, we have 



Prob 



sup 



< e 



-70 4) 



> 1 



l 



(5.35) 



under the same assumption. 

Let us study the condition which realizes such a large (3. Consider first the case with 
A P = 0, and fix the size £ = £ of the box to satisfy the condition (5.32). Further, fix 
AS > and 5± > 0. Assume that the energy E G R satisfies 



£ n -i + IIK) + lloo + K u i + 5-h^ c < E <£„ 



Wr 



lloo 



A_ Ml - AS 



(5.36) 



or 

S n + 1 1 Vq + 1 1 oo + A + mi + AS < E < S n+ i - llV^Hoo - X-Ux - 25 + hoj c , (5.37) 

where S-\ = — oo. We call the energy interval satisfying the condition (5.36) the lower 
localization regime around the band center S n , and call the interval for the condition (5.37) 
the upper localization regime. These conditions imply that we cannot treat the energy E 
near the Landau level S n . From the assumption A P = 0, we have 

£^ i+ <£ n _i + ||y + || (X> , and -Halloo- (5.38) 

From these bounds, one has 

£n- g i,+ + A+mi + Lhu c <E< S c n dgc - A_ Ml - AS (5.39) 

or 

£n d f + Kui + AS < E < S^l_ - \-U\ - 25+huj c , (5.40) 

and so the energy E satisfies the condition (5.5). Since the random potential V w satis- 
fies the condition (4.14) on the ribbon region 1Z, the energy E satisfies the condition in 
Theorem D.7 as 

Sn-i + \\(V + K;) + ||oc + Lhu c < E < S n - \\ (V + V^'W^ - AS (5.41) 

or 

S n + \\(V + \4) + || 00 + AS < E < S n+1 - \\(V + V u )-\\ 0O - 25 + hw c (5.42) 

As a result, we can take (3 = k n t^ oc \f~B for B > Bq^, and the constant C[ n ^ in (5.30) 

is independent of B. Here k n and Bq^ are positive constants which depend only on the 
index n of the Landau level. 

On the other hand, one can choose SE to satisfy (SE)^ 1 ~ Const, x B for a large B 
from the condition (5.33) and K 3 ~ Const, x B for a large B. The asymptotic behavior 
of K 3 is easily derived from the expression (A. 39). See the remark below Theorem A. 2. 
Combining these observations with the bound (5.30), we reach the result that there exists 
a large, positive B^ such that the statement (5.34) holds for all B > B^\ The positive 
constant B^ depends only on the index n of the Landau level. 
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We have fixed the size £ of the box to t = £ at the first, and chosen a large strength 
B of the magnetic field. However, fixing the initial size of the box is not convenient 
for applying the multi-scale analysis to the present system. In fact, we must choose a 
sufficiently large £q for the initial size to satisfy a certain condition which depends on the 
strength B of the magnetic field in the analysis. Clearly, changing the initial size £q to a 
larger value is not allowed for a fixed B because the right-hand side of the bound (5.30) 
depends on the size £ of the box through oc £q +2 . In order to avoid this technical 

difficulty, we take £ to be a function of B as 



i = £ (B) = to 



B/Bf ] 

odd 



< 



(5.43) 



where [x]f dd denotes the largest odd integer which is smaller than or equal to x, and B^ 
is the lower bound for the strength B of the magnetic field which is determined in the 
above; the odd integer £ is chosen so that £ (B^) satisfies the condition (5.32). 
From the definition (5.43) and (3 = k n £~^ ', one has 



with % = k n r 3 /(£ R (n)£ ). Using this, (5.43) and (SE)- 1 ~ Const, x B£l +2 for a large B, 
the right-hand side of the bound (5.30) is bounded from above by 

C^ +11 e-^° = exp [- - (logC)/4 - (2f + II )(\og£ )/£ } 4>] (5.45) 

for a large B. Here C is the positive constant. Thus there exists a large, positive B^ such 
that the statement (5.34) holds for any B > The positive constant B^ depends 

only on the index n of the Landau level. In addition to this, we can choose the constant 
7o — 7cj n ^ so fh & f lo^ is independent of B and of the initial size £q = £o{B). Actually 7^ 
depends only on the index n of the Landau level. 

Next consider the case with A P 7^ 0. In this case, we also take £ = £q(B) of (5.43). The 
decay bound (D.85) for the resolvent in Theorem D.7 was the key to the above argument. 
However, for Ap 7^ 0, we must rely on the different, weaker bound (D.25) in Theorem D.2. 
In fact, we cannot obtain a similar bound to (D.85) because of a technical reason. To 
begin with, let us see the difference between the two bounds. Let E be the energy in 
the spectrum a(H U) ) of the Hamiltonian of the whole system, and let a(H n ) be the 
spectrum of the Hamiltonian H n of (5.13) having the local potential Vr, supported by the 
ribbon region 71. Then the distance between E and a(H-ji) is at most of order of HV^Hoo. 
Namely, 

dist(a(H n ),E) =min{\E+-E\,\E-E-\} < H^IU, (5.46) 

where (E_,E + ) is the spectral gap of H-R. Substituting this into the expression (D.24) of 
(5 in the bound (D.25) for the resolvent, we have that the parameter (5 is at most of 0(1) 
for a large strength B of the magnetic field. Thus we cannot realize a large (5 by taking 
only a large B. 

In order to realize a large f3, we require a strong disorder, together with the strong 
magnetic field. To this end, we take u = Tluj c u with a fixed, dimensionless function u 
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for the random potential V w of (2.6), and choose C = hu c C with a fixed, dimensionless 
constant C as the constant in (D.6). Fix 5± > 0. Assume that the energy E satisfies 



£5 e + hw c L <E< S^t- ~ W+, 



edge 



where 



and 



ccdgc _ c I 
n,+ °n T 



V2e 



L P| I |oo V £-n 



2m P 



|2 



+ ||V r + ||oc + A+ui 



pcdgc 



s, 



V2e 



n+l 



L P|||ooV^-n+l 



lloo 



A_Mi. 



(5.47) 
(5.48) 

(5.49) 



The condition (5.47) implies that we cannot treat the energy E in the interval 

near the Landau level £ n . From (3.7), (3.9) and (4.14), we have 



e. < + < s c n d r and e + > e£»_ - \_ Ul > £:X- 



(5.50) 



Therefore we obtain E + — E > huj c 5 + and E — E_ > hu c 5^. Substituting these into the 
expression (D.24) of (3, we have 



y/2m e 



(hu c ) 3 C 5 + 8- 



h \C (E-E_) + 16(E + + C )(E„ + C ) 



0(B 1/2 ) for a large B, (5.51) 



where we have used Co = Tiu c Cq with a fixed constant Co- In this case, one can easily 
have K 3 = 0(1) for a large B, from (A.3), (A.4), (A.7), (A.25) and (A.39). Therefore we 
can choose {SE)' 1 ~ Const, x £^ +2 to satisfy the condition (5.33). Consequently the same 
statement (5.34) holds for a strong magnetic field and for a strong potential u = 0(B). 



6 Multi-scale analysis 

Starting from the initial decay estimates (5.34) and (5.35) for the resolvent, we derive 
similar estimates for larger scales without losing too much. The main results of this 
section are given in Lemmas 6.2 and 6.3 below. The proofs are given in Appendix F. 
We stress again that these results for large but finite volumes never yield a similar decay 
bound for two arbitrary points in the whole plane R 2 . As to the resolvent in the infinite 
volume, we will rely on the fractional moment method in the next section. The multi-scale 
analysis given here is a simplified version [28, 37, 38] of [26]. Although the method itself is 
well known, we must carefully handle the magnetic field dependence of the decay bound 
for the resolvent again. 

Let i be an odd integer larger than 1, and denote by A^(z) = A^ ra (z) the parallelogram 
box with sidelength I and with center z G T( := £L 2 = {misLi + Tn,n G Z}. The 

distance between two lattice sites in Te is defined by |z| = max{|m|£, \n\£}. Fix a small 
5 > 0. Let Xe( z ) be the characteristic function of the region A^(z), and let xie( z ) be a C 3 , 
positive cut-off function satisfying 



xUz) A4 =1, and supp V^(z) C A 3 ,(z)\A^(z), (6.1) 
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where A^(z) := {r G A^(z)| dist(r, dA 3e (z) > 5}. We write 

R 3e ,.(E+ie) = (H A3l M-E-ie)~\ and W$e(z) = W()&(z))-6V u>3t , 3t )&(z). (6.2) 
Here the Hamiltonian #a m (z) is given by (5.4). 

Definition 6.1 A parallelogram box A 3 ^(z) is called ^-good for some 7 > if the following 
two bounds hold: 



sup 

6^0 



W 5 u {z)R u , x {E + is) Xl {z) 



< e 



and 



sup 



Xi(z)R 3e , z (E + te) (W$(z))*|| <e-*. 



(6.3) 
(6.4) 



Remark: The probability Prob[A 3 £(z) is 7-good] is independent of the center z. 



Lemma 6.2 Let £, be odd integers larger than 1 such that £' is a multiple of £ and 
satisfies £' > 4£. Assume Prob[A 3 ^(- • •) is j-good] > 1 — 77 with a small rj > 0. T/ien 



with 

and with 



Prob[A 3 £/(- • •) is j'-good] > 1 — rf 
r)' = (5e'/£) 4 r) 2 + (£')-* /2 



y = 7 (1 - 4£/f ) - r 1 log^Hf?!) - (f ) _1 (2s + 7) bgf . 

i/ere c a positive constant. 

We define a sequence of monotone increasing length scales £k as 



(6.5) 

(6.6) 
(6.7) 



71/2 



odd 



for A; = 0,1, 2,.... 



(6.* 



where [x]^ dd denotes the smallest odd integer which is larger than or equal to x. Clearly 
we have £k+i > £ 3 J 2 and £ k+ i > 4£ k for all k if the initial scale £ is large enough. 

Lemma 6.3 Take £q = £o(B) which is given by (5.43), i.e., the function of the strength 
B of the magnetic field. Then there exists a minimum strength B > of the magnetic 
field such that 

Prob[A 3 , fc (- • •) is loo -good] > 1 - (4)^ (6.9) 
with some 700 > for any B > B and for any k. 
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7 Fractional moment bound for the resolvent 

As mentioned at the beginning of the preceding section, the multi-scale analysis has the 
disadvantage for the decay estimate of the resolvent in the infinite volume. In order to 
compensate for the disadvantage, we rely on the fractional moment method. The key 
points of the method are that the fractional moment of the resolvent is finite due to 
the resonance-diffusing effect of the disorder, and satisfies a "correlation inequality" [29]. 
But the resolevnt itself without taking the fractional moment cannot be evaluated by 
the method, as we already mentioned above. The aim of this section is to obtain the 
decay bound (7.1) below for the fractional moment of the resolvent for the present system, 
following ref. [29]. Further, the decay bound (7.1) so obtained yields the decay bound 
(7.17) below for the Fermi sea projection P F . In the article [29], the authors showed that a 
decay estimate of a resolvent in the multiscale analysis yields a fractional moment bound 
for the resolvent. In this section, we obtain the fractional moment bound more directly 
from the initial decay estimate which was studied in Section 5. Actually, the initial decay 
estimate was the initial data for the multi-scale analysis in Section 6. 

Consider the present system described by the Hamiltonian with the random po- 
tential V u on a finite region A or the infinite plane R 2 . Let xa, Xb be the characteristic 
functions of the sets A, B with a compact support, respectively. Then the fractional mo- 
ment bound for the resolvent is 

supE x^Ef + is-H^-'xb S < Const.e"^, (7.1) 

£^0 

where E is the expectation with respect to the random variables of the potential V u , 
s G (0, 1/3), (j, is a positive constant, and we have written r = dist(^4, £>). 

We denote by sg(u) the square box centered at u = (^1,^2) £ R 2 with the sidelength 
i, i.e., sg(u) = {r = (x,y) G R 2 | max{|x — ui\, \y — u 2 \} < 2/2}. Consider the Hamiltonian 

H SL{ *o) = -^(P + eA ) 2 + ( y o + K,)| Si(zo) (7.2) 

on the square region sl(z ) centered at z G L 2 with the sidelength L, where we impose 
the Dirichlet boundary conditions, and write the resolvent as 



Rr 



R L (E + te) = (h sl{xo) -E-ie) (7.3) 



for E,e G R. We write Xf( z ) fo r the characteristic function of the square box Sf(z) 
with the sidelength f := ^/3a/2 for z G L 2 . Let SAl(zq) = sls? (zo)\ s L-23f ( z o) ? and let 
z' G 5Al(z ) fl L 2 . In order to obtain the fractional moment bound for the resolvent, we 
want to evaluate 

supE \\\xr(z)R L (E + te)xr^o)\\ S ] for s G (0,1/3). (7.4) 

In the same way as in Section 4, one can find a ribbon region 1Z such that the conditions 
(4.10) and (4.14) are satisfied with probability larger than 

ppcrc = 1 _ C pcrc Le - mp L ^ tw0 positive constantS; C pcTC and m p , (7.5) 
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and that the ribbon 1Z encircles the square box s? (z ), and that the following two conditions 
are satisfied: dist(7£, s? (z )) > and dist(72., Sf(z')) > f° r a ^ z ' e ^Al(zo) flL 2 . Further, 
we can find a C 2 , positive cut-off function X i such that Xi | Sf =(z ) — 1 and supp|Vxi| C C € , 
where the region C e near the center C of the ribbon 1Z is the same as in Section 4. Since 
we can choose X i so that X f( z ')Xi = 0> we have 

Xf(z')R L Xf(z ) = Xf( z ') R LXlXf(z ) = Xr( Z ') R LW(Xl)RLXr( Z o)- ( 7 - 6 ) 

In the same way as in Section 4, we can take the C 2 , positive cut-off function X n an d 
the resolvent R-r for the Hamiltonian H n . Therefore the right-hand side can be further 
written as 

Xf(z')R L Xf( z o) = Xf( z )RLW(xi)R L Xf( z o) 

= X f{ z ')RLXnW(xi)RLXf( z o) 

- -Xr( z ')RLW(xk)RnW( X i)RLXr( z o), (7.7) 

where we have used X%i\n\R. f — 1> Xf( z o)XTi — an d the geometric resolvent equation, 
RlXti = XnRn — RiW T {Xn)Rn- Using the bounds (5.22) and (5.23) for the resolvent R L 
instead of Rm,m'i we obtain 

\\Xr(z')R L Xr(zo)\\ < \\Xr( z ')RLW(xk)XARnXBW( X l)RLXr( z o)\\ 

< l|iW(x^)llllx^xsllll^(xi)^|| 

< h{\E\A\R L \\)f 2 {\E\^\R L \\)\\ XA R nXB \\, (7.8) 

where A = supp|VX^| an d B = C e D supp|Vxi|- In the same way as in Section 4, we have 

\\XaR*Xb\\ < Ci n) e~^ (7.9) 

with probability larger than (1 — C pcic Le~ mpL ), where Cj™' ) and (3 are the corresponding 
positive constants. By using the Wegner estimate, the norm of the resolvent can be also 
evaluated as \\Rl\\ < (^) 1 with probability larger than (1 — Cw\\g\\oo8EL 2 ). We choose 
L to satisfy 

C p C rc Le -m p L < (L/ a )-£/2 (7.10) 

with a positive number £ which we will determine below, and choose SE so that 

CwWgW^SEL 2 = (L/a)"V2. (7.11) 

Then we have 

|| Xr ~(z')ite(zo)|| < C{ n) /i(l^|, (5Ey')f 2 (\E\, {5E)" l )e-^ (7.12) 

with probability larger than (1 — (L/a)~^). We denote by D L the set of the events oo 
satisfying the above inequality (7.12). Note that 



E 



\\Xr(z')RLXr( Z 0)\\ S \ < E [|| X f ( z> ) R L Xr ( z o) \\ ' l(D L 

+ E \\\ X r(z')R LX r(z )\\ s I(D c L ) , (7.13) 
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where 1(A) is the indicator function of an event A. The first term in the right-hand side 
is estimated as 



E 



Xr(z')R L Xr(zo)\\ s l(D L ) < h(\E\,(5E)^)f 2 (\E\,(5E) 



-s/3r 3 



(7.14) 



Using the Holder inequality with s < t < 1, the second term is estimated as 



E 



\Xr(z')RLXr(z )\\ S l(D C L ) 



< E 



IMz^Mzo^E^Dl)] 1 -^ 



< Const. B s x (L/a)- (1 - s/t) t, 



(7.15) 



where B is the strength of the magnetic field, and we have used the fractional moment 
bound eq. (3.19) in [29], and E [I(D£)] < (L/a)~^. The factor B s comes from a careful 
but easy calculation in the fractional moment bound. The positive constant depends only 
on the index of the Landau level in the condition (5.5) for the energy E. 

We take L to be a function of B as L = L(B) = Const, x B 1 ^ 2 . Then the argument of 
Section 5 yields (3r^ = Const, x B 1 ! 2 . We choose t and £ to satisfy (1 — s/t)£ > 3 + 12s. 
Combining these, (7.13), (7.14) and (7.15), we obtain that the quantity, 



B 5s L 3 sup E \\ X r(z')RL(E + ie)xf(zo) 



(7.16) 



becomes small for a sufficiently large B of the strength of the magnetic field. This implies 
that the finite- volume criteria 7 of Theorem 1.2 in [29] is satisfied for magnetic fields whose 
strength B is larger than a positive B . The factor B 5s comes from the 5-dependence of 
the constant in eq. (1.18) of the criteria. Thus the fractional moment bound (7.1) for the 
resolvent holds for such a large magnetic field. 

Let Pp be the projection on energies smaller than the Fermi energy Ep, and let xa-, Xb 
be the characteristic functions of the sets A, B with a compact support, respectively. The 
following Lemma 7.1 is due to [29]. In order to make the paper self-contained, we give a 
proof which is slightly different from that in [29] . 



Lemma 7.1 The following bound holds: 

E WxaPfXbW < Const. exp[-(i dist(„4, B)\. 

Proof: Write R(z) = (z — H^ 1 . Using the contour integral, one has 

1 r E F 1 ry+ 

XaPfXb = — L dE XaR(E + iy-)xB + 7^— J id V XaR(E f + iy)xs 



(7.17) 



+ 



2iri Je 



1 r 

—J e dE xaR(E + /.'/+ ) \ b -f 



2ni 
1 

2tH 



idy XaR(E + iy)xs, (7.18) 



y+ 



where E is a real constant satisfying if w > Eq. The integral near the Fermi energy is 
justified because the operator norm limit, lim e j_ xaR(E ± is)XB, exists [29, 38, 39] almost 
surely for almost every energy E G R. We can choose finite E and y± so that 



\XaR(Eq + w)Xb\\ < Const. e Mr for any real y 



(7.19) 



7 We should remark the following: The condition z' e 5Al(zo) H L 2 is slightly different from that in 
Theorem 1.2 of [29]. In fact, our argument relies on Lemma 4.1 of [29]. 
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and 

\\ X aR(E + iy±)xB\\ < Const.e^ r for E 6 [E , E F ] (7.20) 

with the same decay constant u,. See Appendix D.l for details. Therefore it is enough to 
evaluate the second integral in the right-hand side of (7.18). It is written 



lim inf / dy xaR(E f + iy)xB, 

£n^U Jl(s n ) 



(7.21) 



where {e n } n is a decreasing sequence, and we have written I(e n ) = [y_, e n , e n ). 

Using Fatou's lemma, we have 



E 



rv+ 

/ dy XaR(Ef + iy)Xi 



< lim inf E 

£„-+0 



Hen) 



dy xaR(e f + iy)xB 



(7.22) 



This right-hand side is evaluated as 



E 



/ dy xaR(E f + iy)xs 

Jl{e n ) 



< E 



/ dy \\xaR(E f + iy)xs\\ 

< E/ dy \\xaR(E f + iy)XB\\ s \\xaR(E f + iy)xs\ 

Jl{e n ) 

< E / dy \\ X AR(E F + iy)xB\\ s \y\ s - 1 

J Hen) 



1-s 



< Const.s- 1 (|y + | s + |y_| s )e^ r , 



(7.23) 



where we have used Fubini-Tonelli theorem and the decay bound (7.1) for the resolvent. 
This yields the desired result. ■ 



8 Finite volume Hall conductance 

We recall the previous results of the linear response coefficients [8]. The total conductance 
for finite volume and for t > is written 

atot ' sy[t) \lyy ■ t + 5a yy (t), for s = y. 

Our goal is to give the proof of all the statements of Theorems 2.1 and 2.2. Namely, 
when the Fermi level lies in the localization regime, the Hall conductance a xy is quantized 
to the integer as in (2.26), and both of the acceleration coefficients 7 SJ/ vanish, and the 
corrections 8a sy {t) due to the initial adiabatic process are small as in the bound (2.27). For 
this purpose, we first treat the Hall conductance a xy , and prepare some technical lemmas 
for the Hall conductance a xy for finite volume in this section. 

In the following, we write A = A sys for short. The explicit form of the Hall conductance 
a xy for the finite region A is given by [8] 

o xy = -y-^Tr Pf A [P xA , P yA ], (8.2) 
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where Pp,A is the corresponding Fermi sea projection and 

P s ,a = 77— / dzR A (z)v s R A (z) for s = x,y. (8.3) 
Zm J j 

Here v s are the velocity operators, i.e., (v x , v y ) = v(t = 0) for v(t) of (2.21), and R A — (z — 
Hu,tv)~ l with the finite-volume Hamiltonian H wA with the periodic boundary conditions; 
the closed path 7 encircles the energy eigenvalues below the Fermi level E F . 

Take two rectangular regions f2 and A' so that the following conditions are satisfied: 

flcA'cA, dist(Q, OA') = 5L/2 and dist(A', OA) = 5L/2, (8.4) 

where we have taken the width 5L of the boundary regions as 5L = a(L/a) K with k G 
(1/2, 1) and with L = max{L x , L y }. Here a is the lattice constant of the triangular lattice 
L 2 . Clearly we can take f2 satisfying \Q\ — 0(L 2 ) and | A\0| = O(LSL). We decompose 
u X y into two parts as o xy = o x xy + 0"°^ with 

0% = -^Tr X nP F ,A[P x ,A, P y ,A]xn (8.5) 

and 

0% = -y^-Tr x^f,a[P,,a, P y ,A]Xn, (8-6) 
where xn is the characteristic function of f2, and xh — 1 — Xn- We choose k so that 

84a <SL < 84+ia, (8.7) 

where 4 is a length scale in the sequence {4}fc which is determined by the recursive 
equation (6.8). 

Lemma 8.1 Let A,B be subsets of A. If dist(»4, B) > then the following bound 

is valid: 

\\XaRa(z) X b\\ < Const.L K «- 2 )+ 4 exp [-fi^L 2 ^ 3 ] (8.8) 

with probability larger than (1 — Const. L~ 2 ' K ^ +2 )~ 3 ^ 3 ) 7 where Xa-,Xb are > respectively, the 
characteristic functions of A, B, and /Xqo is a positive constant. 

Proof: In order to prove the statement of Lemma 8.1, we rely on the argument of the 
multiscale analysis in Section 6. Therefore we use the same notations, A^(z), Xi( z )i a / ;L( z )> 
etc. From the assumption dist(»4, B) > 74«/2, there is a sublattice Lg of 4L 2 such that 
B C UueL B ^.4(11) an( i that dist(»4, A 3 ^ fc (u)) > for all u G L B . Using the adjoint of the 
geometric resolvent equation, 

RA(z)xleM) = xi»flsfc,»(*) + Ra(z) (^34 ( u ))* R 3 e k , u (z), (8.9) 

we have 

XARA(z)xe k (u) = X AX S 3e k ( u )Rzt k ,u( z )Xe k ( u ) + XaRa(z) (w^ k (0))* R 3 e k:U (z) X e k (u) 

= xaRa(z) Ru k , u (z) X e k (u) (8.10) 
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for u G L B , where we have used XAxle k ( u ) — which follows from dist(^4, A 3 ^ fc (u)) > 0. 
This yields 



\\ XA RA(z)xe k (u)\\ < \\Ra(z)\\\\ (wi ek (u)Y R MkM (u)\\. 



(8.11) 



On the other hand, we can prove the bound which is given by replacing the operator 
W^ e ,(u) with its adjoint in the bound (F.24) in the same way as in the proof of Lemma F.3. 
Therefore the argument of Lemma 6.3 yields 

{W*(U)Y i?34,uX4(u) <e" Wfc (8.12) 



with the probability larger than (1 — 



k j. From the Wegner estimate (A. 38), one has 
-Ra|| < CVi^3||fi , ||o C /! |A| with probability larger than (1 — £^). From (F.37), we have 

,3/2 



5L < 8£k+\a < 16£l /z a. Immediately, (5L/I6a) 2 ^ < £ k . Substituting these inequalities 
into (8.11), we have 



||x^a(^X 4 (u)II < Const, x £\h 2 exp l-^L 2 ^ 3 



(8.13) 



with probability larger than (1 — 2£^), where /Xqo is the corresponding positive constant. 
The set B is covered by the sets A^(u). The number |Lg| is at most 0(L 2 /£\). Let 
M A denote the event that the bound (8.13) holds for all of the site u e £ k L satisfying 
A ik n A ^ 0. Clearly the probability Prob(MA) that the event Ma occurs, is larger than 
(1 — Const. £^L 2 £^ 2 ). From these observations, one can easily show the statement of the 
lemma. ■ 

Let 5 be a small positive number, and define A 5 := {r e A| dist(r, dA) > 5}. Then one 
has A\A 5/2 = {r E A\ dist(r,dA) < 5/2}. Let xi e C 2 {A) be a positive cutoff function 
satisfying the following two conditions: 



xi 



= 1 and 



xi 



A\A<V2 



= 0. 



(8.14) 



Lemma 8.2 The Hall conductance for the bulk region is written 



h L x L y 



Trx^A^A,*], [P F ,A,y]]Xn + O (exp [-^L 2 ^]) (8.15) 



with probability larger than (1 — Const. L 2 I K ^ +2 ) 3 ^ 3 ) ; where u 1 ^ is a positive constant. 

Remark: Since k(£ + 2) — 3 > from their definitions, the Hall conductance a™ y for the 
bulk region in the infinite volume limit is given by 



xy 



-limT(P FjA ;fi), 



(8.16) 



h £T°o 

with probability one if the limit in the right-hand side exists. Here we have written 

2-77-1 

J(P F , A ;fi) = -r^rTr xnP F A[PF,A,x},[P F ,A,y]]xn- (8.17) 
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Proof: Using the contour integral representation as in (8.3), one has 

1 



Tr Xn-fF,A-Px,A-Py,A 



dzi / dz 2 TrxaRA(z 1 )R/ L (z2)v x R/ L (z2)P y ,AXn- ( 8 - 18 ) 



{2nif 

The integrand is decomposed into two parts as 

Tr XnRAiz^ R A (z2)v x R A (z 2 )Py :A Xn 
= Tr XnRA(z 1 )R A (z2)v x x S A RA(z 2 ) P y , AXn 
+ Tr X nRA(zi)R A (z2)v x (l - Xa)^a(^) Py, AXn- 

For the second term in the right-hand side, we have 

XciRa(zi)Ra(z2)v x (1 - x S a)Ra( z 2) 
< llxn^A^OII \\xa'Ra{z 2 )v x {1 - Xa)Ra{z2)\ 
+ ||xni*A(*i)(l-XA')ll Ra(z 2 )v x (1-x 5 a)Ra(z2) 



(8.19) 



(8.20) 



From (5.18), (8.4), (8.7) and Lemma 8.1, this gives the small contribution. From this and 
the contour integral representation, it is sufficient to consider 



Tr XnRA(z 1 )R A (z 2 )v x x 5 A RA(z2)RA(z 3 )v y x 5 A RA(z 3 )xn- 
Using the identity, v x x A = (i/?i)[H u ,x\xA, one nas 

Ra(z 2 )v x x S a Ra(z 2 ) = ^[Ra(z2),x X {} - ^Ra(z2)xW( X {)Ra(z2). 



(8.21) 



(8.22) 



In the same way, the second term in the right-hand side leads to the small correction. The 
statement of the lemma follows from these observations. ■ 

We denote by Z 2 the rectangular lattice {(6ini, 62^2) | (^1, ^2) £ Z 2 } with a pair b = 
(61,62) of lattice constants, and denote by (Z 2 )* the dual lattice, i.e., (Z 2 )* = Z 2 , — (bi, b 2 )/2. 
Let Sfc(u) be the bi x b 2 rectangular box centered at u = (1x1,1*2) £ Z 2 , and x?<( u ) the 
characteristic function of s&(u). When we consider the characteristic function Xfe( u ) 011 
the region A, we restrict Xfe( u ) to A. 

Lemma 8.3 Let u G Z 2 satisfying s&(u) fl A 7^ 0. Then there exists a positive constant C 
which is independent of the location u and of the size |A| such that 



E[|Tr X6 (u)P F , A [P FiA J][P FiA ,ti]x 6 (u)|] < C, 
where Jj is either x or y. 
Proof: Note that 

E[|Tr X6 (u)P F , A [P F , A ,x][P FiA ,y] X6 (u)|] 



5.23) 



< 



E 

v,wGZf: 

s 6 (v)nA^0, s ft (w)nA^0 



E [| Tr Xb (u) P f ,aX6 ( v) [P f ,a , x] Xb ( w) [P F , A , y] Xb (u) 



5.24) 
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In order to estimate the summand in this right-hand side, we introduce the two component 
function (x b ,y b ) of r which is defined by (x b ,y b ) = (ui,u 2 ) for r G s&(u) with u G Z 2 ,. 
Using this function, one has 



[Pf,a, A = [Pf,a, i? - + [ft,A, a; 6 ] and [P F , A , y] = [P F , A , (y - y 6 )] + [P F , A , y 6 
Further we note that, for any bounded operator A, 



1.25) 



Tr Xfe (u)P FiA x b (v)A| < JTr X 6(u)P F ,AX6(u)- JTt AA* X b(y)P ¥AXh (y 



< \\A\\^Tr X b(u)P F ,AXb(u) ■ V Tr Xb(v)P F ,AXb(v) 

< Const.||A||, 



5.26) 



where we have used 
Tr X£ (u)P FXe (u) = Tr Xe (u) 



(H + C )P F (H + C ) 



H + C 



H + C 

Tr Xe ( u )— 1— (P w -V u + C )P F (P W - K, + C ) ! 



Hq + C 



Po + C 



< Const.Trx £ (u) ( — - ) Xe (u) < oo. 



Xe(u) 

(8.27) 



Here we have taken a real number Co satisfying Po + Co > 0, and used ||P W P F || < oo. 
Therefore the statement of the lemma follows from the norm bounds such as 



and 



X&(v)[P F , A ,z 6 ]x&(w)| < \wi - ^i|||x&(v)Pf,aX&(w)||, 
Xb(v)[P F ,A, (a: - x b )]xb(w) < Const. ||x&(v)P F ,aX&(w) 



E[||x 6 (v)P Fi aX 6 (w)||||x6(w)P F) aX 6 (u)||] 

< {E[||x b (v)P F ,AX 6 (w)|| 2 ]} 1/2 {E[||x b (w)P Fi AX6(u)|| 2 ]} 

< Const. e -M|v-w|/2 e - M |w-u|/ 2) 



(8.28) 
(8.29) 



1/2 



(8.30) 



where we have used Schwarz's inequality, ||Xft( v )-PF,AXf>( w )ll < 1> an d the bound (7.17) for 
the Fermi sea projection. ■ 

In the same way, we obtain 

Lemma 8.4 Let u, v G 7? h satisfying s&(u) fl A ^ and s&(v) fl A ^ 0. Then there exists 
a positive constant C which is independent of the locations u, v and of the size |A| such 
that 

E[|Trx 6 (u)P F , A [P F ,Aj[P F ,Ajx 6 (u) •Trx 6 (v)P F , A [PF,A,tt][ J PF,A,tt]x 6 (v)|] < C, (8.31) 
where Jj is either x or y. 
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Using the magnetic translations and the argument in the proof of Lemma 8.2, the Hall 
conductance for the boundary region is written as 



2m 



(jOUt 

h L T L. 



£lYx^, A [[P^y,x], [P^ A ,y]]xn> + O (exp [-/4L 2k/3 ]) (8.32) 



with probability larger than (1 — Const. L~ 2 ^ +2 ^~ 3 ^ 3 ) , where Q' is the translate of a 
portion of the boundary region A\Q, and Pp A is the corresponding translate of the Fermi 
sea projection. From this result, the Hall conductance a™ 1 for the boundary region is 
vanishing in a limit L j oo with probability one: 



Theorem 8.5 There exists a sequence {L n } n of the system sizes L = L n such that 



o~ 



xy 



as L n — ► oo for almost every u>. 



Proof: Write 



I'{L) 



L X Ly n 



J2\^Xn'P^ A m, A ,x],[P^ A ,y}}xn' 



(8.33) 



By Lemma 8.3, we have E[J'(L)] — > as L — > oo. Combining this and the inequality 
E[J'(L)] > eProb(J'(L) > e), we can find a sequence {L n } n of the system sizes L = L n 
such that {L n } n satisfies the following two conditions: 



£Prob(/'(L n )>e n )<oo and £ L -2[*«+2)-3]/3 < QO> 



(8.34) 



where {e n }„ is a sequence satisfying e n — * as n — > oo. The application of Borel-Cantelli 
theorem yields that for almost every a;, there exists a number n (c<->) which may depend 
on oo such that I'(L n ) < e n for all n > n (uj) and that the finite size correction for is 

evaluated by O (exp -/i^I^ 3 ^ for all n > n (uj). By combining this with the expression 
(8.32) of a°y l : the statement of the theorem is proved. ■ 



9 Integrality of the Hall conductance 
— Index theoretical approach — 

In this section, integrality of the Hall conductance is proved by using the index theoretical 
method [10, 17, 18]. 

When we apply the method of [10] using a pair index of two projections to a concrete 
example of a continuous random model such as the present system, there arises a problem 
that we need a decay bound for the integral kernel of the Fermi sea projection whose 
Fermi energy lies in the localization regime. But getting such a decay bound is very 
difficult, and so this problem is still unsolved. Recently, Germinet, Klein and Schenker 
[22] proved the constancy of the Hall conductance for a random Landau Hamiltonian which 
is translation ergodic, without relying on a decay bound for the integral kernel of the Fermi 
sea projection. In their proof, they used a consequence of the multiscale analysis which 
is related to multiplicity of the eigenvalues of the Hamiltonian, for the Hall conductance 
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formula 8 which is expressed in terms of switch functions instead of the position operator 
of the electron. This Hall conductance formula was justified [7] within the linear response 
approximation under the assumption on a spectral gap above the Fermi level. The integer 
of the quantized value of the Hall conductance can be determined under the assumption 
that the disordered-broadened Landau bands are disjoint, i.e., there exists a nonvanishing 
spectral gap between two neighboring Landau bands. 

In our approach, we assume neither the above disjoint condition for the Landau bands 
nor the periodicity of the potentials Vq and Ap which implies a translation ergodic Hamil- 
tonian. But we must require the "covering condition" (2.10) which is not required in [22]. 
This condition is needed to estimate the number of the localized states and to obtain a 
decay bound [29] for a fractional moment of the resolvent. In order to circumvent the 
above problem about a decay estimate for the integral kernel of the Fermi sea projection, 
we introduce a partition of unity which is a collection of the characteristic functions of a 
small rectangular boxes. 

Let s £ (u) be the £i x £ 2 rectangular box ceneterd at u e (Z 2 )* with the pair e — (£i, £2) 
of the sidelengths, and Xe( u ) the characteristic function. We take N 1 e 1 = a and N 2 e 2 = 
\/3a/2 with large positive integers N 1 , N 2 and with the lattice constant a of the triangular 
lattice L 2 so that the set of all the boxes s e (u) is invariant under the lattice translations 
of the triangular lattice L 2 . When we take the limit £1,62 I 0, we keep the ratio £1/52 
finite. We introduce a unitary operator, 

K= E X e (u)exp[i0 a (u)] foraeZ 2 , (9.1) 

ue(z2)* 

where a (u) is the angle of sight from a to u, i.e., arg(u — a) in the terminology of the 
complex plane. Consider the operator, T := P F — U^P F (C/|) 



Lemma 9.1 For fixed parameters e 1: e 2 , E [Tr |T| 3 ] < 00. 
Proof: Note that 

T = P F - UiP F (UlY = £ Xe (u) [P F - KP F (Ul)*} X e(v) 

u,v G (Z§)* 

E [l-e^ (u) -^ (v) ]x £ (u)P F X £ (v). (9.2) 



u,v G (Z§)* 

Define t U)V := 1 - e *a(u)-tf a (v) and T ^ ;= t UivXe (u)P F x e (v). Following the idea of [18], 
we introduce = T UjV 5 u _ bjV . Clearly, one has Eb^S = ^u,v and 

(T^ ^ ^ ^ ' ^w.n^w-b.u^w.v^w-b.v 

' w 

~ J u+b,u i u+b,u°u,v 

= \t u+hyU \ 2 Xe{u)P F Xe{u + b)P F Xe{u)5^ v . (9.3) 



8 The explicit form of the Hall conductance formula using switch functions is given in Appendix I. 
We also discuss the relation between this and the standard Hall conductance formula using the position 
operator instead of the switch functions. 
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Using these identities and Minkowski's inequality, one obtains 



(ETr|T| 3 ) V3 < ^(ETr|T( b f) 



1/3 



< E(ElWb,u| 3 E[Tr \ X e(u)P F Xe(u + b)P F XsH\ 3/2 

b I u 



1/3 



Since the inequality, 



I — g»0a(u)-i0 a (v) 



< 



2|u - v| 



u — a 



holds as in [10], one obtains 



(9.4) 



(9.5) 



ElWb, u | 3 < 2 3 E^T? ^ Const.|b| 3 . 



u — a 



Note that 



Tr |x £ (u)P F x £ (u + b)P FX£ (u; 



1 3/2 



< ^TTXeHP F Xe(n + b)P F Xe(u) 

x y'Tr Xe (u)P F x £ (u + b)P F x £ (u)P FXe (u + b)P FXe (u) 



(9.6) 



< V / TV X£ (u)P fX£ (u) • ||x £ (u)P F x £ (u + b)|| ^Tr X£ (u)P F x £ (u) 

< \\Xs(u-)P F Xe(u + b)|| Tr x £ (u)P F x £ (u). 

From this, the decay bound (7.17) for the Fermi sea projection and (8.27), we have 



E 



Tr |x £ (u)P F x £ (u + b)P FX£ (u)| 3/2 <Const.e^ b l 



Combining this, (9.4) and (9.6) yields 



(E Tr |T| 3 ) 1/3 < Const. ^ |b|e^ |b|/3 < oo. 



(9.7) 



(9.S 



(9.9) 



The result implies that the operator T 3 is trace class for almost every u. Thus we can 
define the relative index [10], 



Index(P F , UiP F {UtY) := Tr (P F - t/|P F (^)*) 3 , 



(9.10) 



for the pair of the projections. This right-hand side takes an integer value as proved in 
[10] - 

Let £p be a large positive integer. Let A LP e C^R 2 ) be a periodic vector potential 
satisfying the periodicity, 



A (r + £ pai ) = A (r + £ P a 2 ) = A (r), 



(9.11) 
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and Vq F a periodic electrostatic potential satisfying the same periodicity, 

V LP (r + £ pai ) = V LP (r + 4>a 2 ) = V LP (r), (9.12) 

where slj are the primitive translation vectors of the triangular lattice L 2 . In order to 
prove integrality of the Hall conductance, we consider the Hamiltonian, 

H]f = 7^~[P + e(A LP + A )] 2 + V LP + V u , (9.13) 

on the whole plane R 2 . Namely this Hamiltonian is obtained by replacing Ap, Vq with 
A LP , Vq P in the Hamiltonian H u of (2.1). We choose the integer i F so that the unit cell 
of the large triangular lattice £pL 2 contains the rectangular region A sys of (2.3) on which 
the present finite Hall system is defined. 

We note that the magnetic translations act on the random potential V w as the corre- 
sponding translation, and that the pair index does not depend on the location a of the 
flux [10] . Therefore the pair index for the Hamiltonian H]f is an invariant function of the 
randomness under the lattice translations of the triangular lattice £pL 2 . Further, since the 
index is measurable [10] and integrable with respect to the random variables, Birkhoff's 
ergodic theorem implies that [18] the value of the pair index does not fluctuate in the 
sense that it takes an integer given by its mean for almost every random potentials. But 
the value of the integer may depend on the period £p. 

As in [10], the relation between the pair index, Index(P F , UlP F (U a )*) , and the Fredholm 
index, Index (P F U a P F ), of the operator P F U a P F in range of P F is given by 

Index(P F [/|P F ) = -Index(P F , U £ a P F (U £ a )*). (9.14) 

Consider another unitary operator, 

u a := * + *y-^ + ^\ (9 . 15) 

\x + iy — (ai + ia 2 ) | 

Clearly, one has Ppf/ a Pp = Ppf/^Pp + P-p{U a — U a )Pp. Since the operator of the second 
term is compact, stability theory 9 of the Fredholm indices implies that PpUaPp becomes a 
Fredholm operator, too, and that the index is invariant under the compact perturbation. 
Thus we have 

Index(P F f/ a P F ) = Index(P F f/fP F ) = -Index(P F , U £ a P F (U £ a )*) . (9.16) 

In consequence, the pair index does not depend on the parameters £i,£2- 

Following [20], we obtain the expression (9.26) below with (9.27) for the pair index. 

The expression leads to the well-known Hall conductance formula [17] which is written in 

terms of the position operator of the electron. 
To begin with, we note that 

Tr(P F -^P F (^)*) 3 

= J2 Tr xsH(p f - u £ p F (u £ a y) X e(v)(p F - u a p F (u £ a y) X sM 

u,v,w 

= E tu,vtv,wtw,uTrx £ (u)P F x e (v)P F x £ (w)P FX£ (u). (9.17) 



9 See, for example, the book [40]. 
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Since the index is independent of the location a of the flux [10], one has 



index(P F f/ a P F ) = E E tu,vtv,wtw,uTrx £ (u)P F x e (v)P F x e (w)P F x £ (u), (9.18) 

aGA^ u,v,w 



where A e = e 1 {-£, -£+ 1, ...,£} x e 2 {-£', -£' + 1, £'} C Z 2 £ , and V e = (2£+ l)(2f + 1). 
We choose so that the ratio is finite. 

Lemma 9.2 There exists a sequence {£ n } n such that for almost every uj , the index is 
written 

Index(P F [/ a P F ) = - lim — V EE *u,v*v,w*w,u'S , u,v,w,u (9-19) 

MooV<„ U6A . v,w aeZ 2 

TY X£ (u)P F x £ (v)P F x £ (w)P FX£ (u), (9.20) 

where the lattice A| gwen fry 
A* = + 1/2, + 3/2, ...,£- 1/2} x e 2 {-£' + 1/2, + 3/2, ...,£'- 1/2}. (9.21) 

The proof is given in Appendix G. Using Connes' area formula [41], 

2-jri 

E Wv,wt w ,u = (v - u) x (w - u), (9.22) 

the index of (9.19) is written 

index(P F [/ a P F ) 

,. 27r 
= — lim - — 

|H|Too 



2ni 

ni 1 ? 1 ToT 51 E( V_U ) X ( W - U ) Tr Xe( U ) P FX £ (v)P F X £ (w)P F X £ (u), 



ueA; v,w 



(9.23) 



where Q is the L\ x Z-2 rectangular box centered at r = with the sidelengths L\ = 1£ n E\ 
and L 2 = 2£' n e 2 . Note that 

(v — u) x (w — u) = (v — w) x (w — u) = (t> 2 — w 2 )(u>i — ui) — (vi — wi)(w 2 — u 2 ). (9.24) 

Further we have 

(v 2 - w 2 )(wi - Ui)Tr x £ (u)P F x £ (v)P F x £ (w)P F x £ (u) 

= Tr X£ (u)P FX£ (v)[X 2 £ ,P F ] X£ (w)[X 1 £ ,P F ] X£ (u), (9.25) 

where the two-component function (Xf,X|) of r = (x,y) is given by (Xf,X|) = (1x1,1*2) 
for r in the e 1 x e 2 rectangular box s e (u) centered at u = (1*1, *u 2 ). From these observations, 
the index is written 

Index(P F [/ a P F ) = lim J £ (P F ;fi,£ P ) (9.26) 

fi|too 
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for almost every u, where we have written 

27T7 

T(P F ;Q,£ P ) = —TT X nP F [[PF,Xl], [P F ,X|]] Xn . (9.27) 

From the proof of Lemma 9.2, we obtain 

E[\Index(P F U a P F ) -T(P F ;tt,£ P )\] -> as \Q\ | oo. (9.28) 

We also write 

27T? 

Z(P F ;fi,£ P ) := — IVx n PF[[PF,^i],[PF,^2]kn, (9-29) 

where (Xi,X 2 ) = (x,y). This right-hand side is nothing but the well-known form of the 
Hall conductance formula [17]. 

Lemma 9.3 The following bound is valid: 

E[\l(P r ,n,£ P ) - T(P F ;tt,£ P )\] < Const. |e|, (9.30) 

where \e\ = \Js\-\- e\, and the positive constant in the right-hand side is independent of Vt 
and £p. 

Proof: Note that 

|Tr X nP F [P F , X x \ [P F , X 2 ] X n ~ Tr X nP F [P F , X{] [P F , X|] Xn | 
< |Tr X nP F [P F , (X, - Xf)][P F , X 2 ] Xn \ + |Tr Xn P F [P F , X[][P F , (X 2 - X|)] X n| , 

(9.31) 

[P F , X 2 ] = [P F , (y - y b )] + [P F , y b ] and [P F , X{] = [P F , {X{ - x b )} + [P F , x b ]. (9.32) 

Therefore we can prove the statement of the theorem in the same way as in the proof of 
Lemma 8.3. ■ 

Since we can apply Borel-Cantelli theorem as in the proof of Theorem 8.5 to this result, 
Lemma 9.3 yields that there exists a sequence {e n = (£i >n , e 2 ^ n )} n of e = e n satisfying 
£ j,n — * as n — > oo, for j = 1, 2, such that 

\i(p F -,n,e P ) - T n (p F -,n,e P )\ -> o as n^oo (9.33) 

almost surely for any fixed large Q. 

Let {S n }^ =l be a sequence of positive numbers 8 n satisfying S n — > as n — * oo, and 
let {p n }^=i be a sequence of positive numbers p n < 1 satisfying J2 n Pn < oo. Relying on 
Lemma 9.3, we choose e = e n = (£i, n , £ 2,n) for each n so that 

Prob[\l(P F ;Q,£ P ) - T-(P F ;Q,£p)\ > 5 n /2]<p n . (9.34) 

Further, for this e = e n , we can choose a sufficiently large f2 = Q, n so that 

Prob [|Index(P F f/ a P F ) - J £ "(P F ; fi n , £ P )| > 5 n /2] < Pn , (9.35) 
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from the proof of Lemma 9.2. The application of Borel-Cantelli theorem yields that for 
almost every u, there exists a number no(u) such that for all n > n (u), the following two 
inequalities are valid: 

\l(Pr,tt n J P ) -T"(P r ,n n ,£ P )\ < 5 n /2. (9.36) 

and 

|Index(P F f/ a P F ) - T"(P F] n n ,£ P )\ < 5 n /2. (9.37) 
These two inequalities imply 

|Index(P F f/ a P F ) -l(P F ;Q n ,£ P )\ 

< |Index(P F f/ a P F ) - l £ -(P F ;n n Jp)\ + \l £ "(P F ;n n ,£p)-l(P F ;n n ,£ P )\ 

< 5 n /2 + 5 n /2 = 5 n . (9.38) 

This result is summarized as the following index theorem: 

Theorem 9.4 For a fixed period £p and for almost every uj, there exists a sequence {Vt n } n 
of the regions Q = Q n going to R 2 as n — > oo such that 

2ni 

Index(P F f/ a P F ) = lim — - Tr X n n P F [[P F , X,}, [P F , X 2 ]} X n n - (9-39) 

For taking the infinite-volume limit A j R 2 , we want to take a sequence {A n } n of 
the finite region A = A n of the present system so that the condition (8.4) is satisfied for 
A = A n D Q = Q n . Then we must take the limit £p j oo together with the limit f2 n j R 2 
so that the unit cell of the lattice £pL 2 includes the region A = A n of the present system. 
In the following, we consider a sequence {A n , Q n , £p >n } n which satisfies this requirement. 

Since both of the key bounds in the proofs of Lemmas 9.2 and 9.3 do not depend 
on the period £ F , we can take the limit £p ] oo together with the limit Q n j R 2 in the 
above argument for Theorem 9.4 so that the above requirement is satisfied. But both of 
Index (P F U a P F ) and X(P F ; £l n , £p) may go to infinity as £p \ oo. First let us prove that this 
case does not occur. From the argument of the proof of Lemma 8.3, one can easily show 
that the expectation value E[|X(P F ; Q, £p)\] is bounded uniformly in Q and £ F . Combining 
this with Fatou's lemma, we have 



E 



liminf \X(P F ;Q,£ P ) 
QTRVptoo 



< liminf E[\l(P F ;Q,£p)\]<oo. (9.40) 

SZTR/^pToo 



This implies that for almost every u, there exists a sequence {Q n (u), ^p,„(^)} n of the 
pair {Q,£ P } such that {Q n (u)} is a subsequence of the sequence {tt n } of Theorem 9.4, 
and that lim n | 00 X(P F ; fl n {u), £ Ftn (uo)) exists. Here we should stress that the sequence 
{Q n (u) , £p t n(u)} n may depend on the random event uj. On the other hand, the inequality 
(9.38) holds for a large pair £ P } = {f2 n (cj), £ P:fl (uj)}. These observations imply that for 
a fixed uj, the index Index(P F ?7 a P F ) converges to an integer as n j oo, too. But, since the 
index does not depend on uj as mentioned above, we can write {£p 5 „}„ for the sequence 
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{£p, n (ui)} n by dropping the u dependence, and obtain the result that the following limit 
exists and is constant for almost every uj: 

Inde Xoo (P F f/ a P F ) := lim Index(P F [7" a P F ). (9.41) 

Newly we choose {Q,£p} = {VL n ,£p n } in the inequality (9.38). Then, since the index 
converges to the integer as n j oo, we obtain 

Theorem 9.5 For almost every uj, there exists a sequence {Q n ,£p,n}n of the pair {Q = 
Q n ,£p = £p,n} such that the following relation holds: 

2iri 

Inde Xoo (P F C/ a P F ) = Km — - Tr Xa P F [[P F , [P F , X 2 ]] X n n . (9.42) 

S2„TR ,«p,„T°° |"n| 

Theorem 9.6 There exists a subsequence {Q n ,£p,n}n of the sequence of the preceding 
Theorem 9.5 such that 

2ni 

Inde Xoo (P F f4P F ) = Km — - Tr X n n P F [[P F , X ± ], [P F , X 2 ]} X n n 

S2„tR- ,tp,nT°o |"n| 

2iri 

= Km 2 T^Tr X n n PF ) A„[[PF,A„,x],[P F , An ,y]] X n n (9.43) 

with probability one. Here we take A = A n and £p = £p tH so that the region A = A„ of 
the present system satisfies the condition (8.4) for each Q = Q n , and the unit cell of the 
lattice ^p, n L 2 includes the region A n . 

Theorem 9.6 follows from the following lemma: 

Lemma 9.7 Under the same assumption as in Theorem 9.6, we have 

E[\I{P F ;Q n ,£p, n )-l(P FAn ;Q n )\] ^0 as n -> 00. (9.44) 

The proof is given in Appendix H. 

Theorem 9.8 There exists a sequence {L n } n of the system sizes L = L n such that the 
Hall conductance a xy in the infinite volume limit exists and is quantized to an integer as 

a xy = ^Inde Xoo (P F f/ a P F ) (9.45) 

for almost every uj. 

Proof: From Theorems 9.5 and 9.6, there exists a sequence of {A n , Q n , £p >n } n of the triplet 
{A,Q,£p} such that the following three conditions are satisfied: (i) the condition (8.4) is 
satisfied for A = A n D Q = Q n , (ii) the unit cell with the period £p^ n includes the region 
A = A n of the system with the linear size L = L' n) and (iii) for almost every u, the 
following formula holds: 

27ri 

Inde Xoo (P F ?7 a P F ) = lim Tr XQn P F , A J[P FiAn , x], [P F;An , y]} X n n . (9.46) 
"Too \\l n \ 

Take a subsequence {L n } n of {L' n } n so that the sequence {L n } n of the system sizes satisfies 
the two conditions of (8.34) in the proof of Lemma 8.5. Then, for almost every uj, the 
contribution of the Hall conductance from the boundary region is vanishing as n | 00, 
and the correction of the Hall conductance cr™ of (8.15) for the bulk region is also vanishing 
in this limit. Combining this with (8.15) yields the desired result. ■ 
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10 Constancy of the Hall conductance 
— Homotopy argument — 

In this section, we prove that the Hall conductance o xy is constant as long as both the 
strengths of the potentials and the Fermi energy vary in the localization regime. 

10.1 Changing the strengths of the potentials 

Consider first changing the strengths of the potentials V LP ,A LP ,K, in the Hamiltonian 
Hlf of (9.13) on the whole plane R 2 . In order to prove constancy of the Hall conductance, 
we extend the homotopy argument of [19] for lattice models to continuous models by relying 
on the fractional moment bound [29] for the resolvent. As a byproduct, we prove that the 
quantized value of the Hall conductance is independent of the period £p of the potentials 
\/ LP , A LP . See Theorem 10.1 below. 

Since all the cases can be handled in the same way, we consider only the case where 
the strength of the vector potential Ap varies. We denote by P F the Fermi sea projection 
for the Hamiltonian H' u with the vector potential A p = A P + SA P with a small change 
|| |<5Ap] Hoo- Since the index does not depend on u, we have 

|Index(P F f/ a P F ) - Index(P F f/ a P F ) | = E [|Index(P F C/ a P F ) - Index(P F f/ a P F ) |] 

< E [|Index(P F f/ a P F ) - X e (P F ; Q,£ P )\] 

+ E[\i £ (p^n,i P ) -i%p r ,n,£p)\] 

+ E[\r(Pr,n,£ P ) -Index(P F C/ a P F )|]. (10.1) 

From (9.28), the first and the third terms in the right-hand side become small for a large 
Q. Therefore it is sufficient to show that the second term become small for a small change 
III^AplHoo of the vector potential. 

Relying on the expression given by the right-hand side of (9.23), let us estimate the 
difference, 

Tr x £ (u)P f x £ (v)P f x £ (w)P fX£ (u) - Tr x £ (u)P F x £ (v)P F x £ (w)P F x £ (u) 
= Tr x £ (u)AP fX£ (v)P f x £ (w)P fX£ (u) + Tr x £ (u)P f x £ (v)AP fX£ (w)P fX£ (u) 
+ Trx £ (u)P FX£ (v)P F x £ (w)AP FX£ (u), (10.2) 

where AP F = P F — P F . The first term in the right-hand side is estimated as 

|Tr X£ (u)AP FX£ (v)P F x £ (w)P F x £ (u)| < ^Tr Xe (u) AP fX£ (v)P f x £ (v) AP fX£ (u) 

x y/Tr x £ (u)P f x £ (w)P f x £ (w)P fX£ (u) 

< Const.||x £ (u)AP FX£ (v)||||x £ (u)P F x £ (w)|| 

(10.3) 

by using Schwarz's inequality and the bound (8.27). Similarly, the second term is estimated 

as 

|Tr X£ (u)P F x £ (v)AP FX£ (w)P F x £ (u)| 
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< jTr X e(u)P FX e(u) 



x ^Tr X£ (u)P^(w)AP FX£ (v)P F x £ (v)AP FX£ (w)P^(u) 

< Const. ||x £ (u)P F x £ (w)AP FX£ (v) || 

< Const.||x £ (u)P F x £ (w)||||x £ (w)AP FX£ (v)||. (10.4) 

The third term can be handled in the same way as for the first term. 
The contour representation of the Fermi sea projection yields 



with 



and 



| X£ (u)[P F - P F ] X£ (v)|| < J« + + J# + J# (10.5) 

|x £ (u)[P'(P F + iy) - R(E F + iy)] Xe {v)\\ , (10.6) 
||x £ (u) [R'(E + iy) - R(E + iy)] X s{v) || (10.7) 



(10.8) 



IW = 




I ' 


u,v 


27T . 




7(2) _ 


1 


/•S/+ 


u,v 


" 2^ 


•A/- 


= 




' C 


u,v 


2^ J 


Eo 



where P'(^) = ( z ~ K)' 1 and = ( z ~ H^)' 1 . 

First let us estimate the last three integrals except for Note that R'(z) — R(z) = 
-R'(z)SH 0J R(z), where 

2 

bH w = ^[5A • (p + eA) + (p + eA) • 5 A] + ^|5A| 2 . (10.9) 

Since all the contributions in the perturbation SH^ can be handled in the same way, we 
consider only 

Xe(u)R'(z)5A s (p s + eA s )R(z) X e(v) 
= £ X e(u)R'(z)x b (u')SA s (Ps + eA s )x 5 b (u>)R(z)xe(v) (10.10) 

u' 

as a typical one. Here {Xb( u )} u is the partition of unity which is given in the proof of 
Lemma H.l in Appendix H, and x&(u) is the characteristic function of the support of 
Xb( u )- The norm is estimated as 

\\Xe(u)R'(z)5A s (p s + eA s )R(z)xe(v)\\ 

< ||M s || 0O ^||x £ (u)P / (z)x 6 (u / )||||(p s + e^)x^(u / )P(^)x £ (v) 

u' 

< Const.||5A|| 00 ^||x £ (u)P / (z)x 6 (u / )|| 

u' 

\x b (u')R(z) X e(v)\\ + Const. ||x6(u')P(^)x £ (v)|| 1/2 ] , (10.11) 

where we have used Lemma H.l for getting the second inequality. Since dist(a(P w ), z) > 
in the present situation, all the norms about the resolvent R(z) decay exponentially at the 
large distance. Therefore we obtain 

/W v < Const.|||5A P ||| 00 exp[-//|u- v|] (10.12) 
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5A P 


|| oo 









with a positive constant //, where (j = 2, ±. Consider the contribution from (10.3) because 
the rest can be treated in the same way. The corresponding contribution is estimated by 

£ £ |u - v||u - w|e-^- v lE[|| X£ (u)P^ £ (w)||] 

uGA| v,w 

< Const. || ItfAplHoo, (10.13) 

where we have used the decay bound (7.17) for the Fermi sea projection. 
Let s G (0, 1/3). The rest of the integrals is written 

= ^/^^iixsH^F+^-^+^wir 73 

x \\x £ (u)[R'(E F + iy)-R(E F + iy)} Xe (y)\\ 1 - s/3 
< - f V+ dy \\xe(u)[R\E F + iy) - R(E F + iy)]* e (v)|| s/3 \y\ s ^'\ (10.14) 

7T Jy_ 

where we have used the inequality (J^jdjY < ^j a j f° r s e (0> 1) an d % > 0, and the 
inequality \\R ( (E F + iy)\\ < \y\~ l for R i = R',R. For the norm of the operator in the 
integrand, the contribution from the term (10.10) can be estimated as 



\Xs(u)R'(E F + iy)SA s {p s + eA s )R{E F + iy)xs{v) 



is/3 



< Co™t]\5AAUY,{\\xs{n)R\E F + iy)U^)\\ sl * \\xMR(E F + iy)xe(y)\\ s/3 

u' 

+ Const. \\ X e(u)R'(E F + iy )x b (u')\\ s/3 fo(u')P(P F + iy) Xe (v) f /6 } , (10.15) 

where we have used Lemma H.l for getting the inequality. 

Consider the contribution from the first term in the summand in the right-hand side 
of (10.15) because the second term can be handled in the same way. The corresponding 
contribution from (10.3) is estimated by 





5A P 


|| oo 









J2 J2J2\ U ~ v ll u ~ w l E f v+ dy\y\ 

ueA;v,w u ' J y- 



|s/3— 1 

'Ml 5 / 3 II ~ I >\Tf(T? I • \ / \ 1 1 s/3 



x \\ X e(u)R'(E F + iy )x b (u')\\ s/ ' i \\x b (u')R(E F + iy) Xe (v)\r ||x £ (u)P fX£ (w)||. 

(10.16) 

Using Holder's inequality, we have 

E [\\ X e(u)R'(E F + ^)x,(u')ir /3 \\xMR(E F + ^y)Xe(v)\\ s/3 \\ X s(u)P fX sM\[ 

< {E [\\xe(n)R\E F + iy)Mu') f] } V3 {E [\\ Xb (vi')R(E F + iy) Xe (v) f] } V3 

x{e[|| X£ (u)P fX£ (w)|| 3 ]} 1/3 

< Const.e^ |u - u ' l/3 e^ |u '- v|/3 e-^ |u - w|/3 , (10.17) 

where we have used the decay bounds (7.1), (7.17) and ||x £ (u)P F x £ (w)|| < 1. Relying 
on Fatou's lemma and Fubini-Tonelli theorem, and substituting the bound (10.17) into 
(10.16), we can obtain the desired result. 
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Thus the index is constant as long as the strengths of the potentials vary in the local- 
ization regime. In order to describe our statement more precisely, we recall the definitions 
of the lower and upper localization regimes which are given by the intervals (5.36) and 
(5.37), respectively, in the case of Ap = 0. In the case of Ap ^ 0, the corresponding 
condition is given by (5.47). See also Section 11. We continuously change the strength 
of all the potentials, V , V u and Ap, starting from the special point, ||Vo||oo = ||K;||oo = 
and HIApllloo = 0, for a fixed Fermi energy Ep. Then, if the Fermi energy Ep is staying 
the lower or upper localization regime, the index is equal to the special case that all the 
potentials are vanishing. This result also implies that the index does not depend on the 
period £p of the potentials. In consequence, Theorem 9.8 is refined as 

Theorem 10.1 Suppose that the Fermi energy Ep lies in the localization regime around 
the n-th Landau energy £ n -i — (n — l/2)ftij c . There exists a sequence \^L x ri) Ly n ^ n of the 
system sizes such that the Hall conductance a xy in the infinite volume limit exists and is 
quantized to an integer as 

— 6 x { n f° T ^ C u PP er l° ca ti za ti° n regime , , 

xy h I n — 1 for the lower localization regime 

with probability one. 

Remark: 1. When the strength of one of the potentials becomes sufficiently large for a 
fixed strength of the magnetic field, the localization regimes become empty in our defini- 
tion. Thus we need the condition that the strengths of the potentials are weak, compared 
to the strength of the magnetic field. 

2. The number of the states in a localization regime will be proved to be of bulk order for 
the weak potentials, compared to the strength of the magnetic field in Section 11. 

3. We do not require any assumption on the tails A G [A m i n , — A_] U [A+,A max ] of the 
coupling constant of the random potential V u . Therefore we allow the possibility that the 
spectral gap between two neiboring disordered-broadened Landau bands vanishes owing 
to the tails of the random potential. 



10.2 Changing the Fermi level 



Next let us prove the constancy of the Hall conductance for changing the Fermi level Ep. 
The Hall conductance a xy of (8.2) is written 



ihe 2 



0~xy 



E 



LxLy m ,n:E m <E F <E n 



(iPm,V x ip n ) (V? 

\E m -En) 



(x <-> y) 



;io.i9) 



in terms of the eigenvector tp n of the single electron Hamiltonian H w with the eigenvalue 
E n , n = 1, 2, . . ., on the box A sys . We take the energy eigenvalues E n in increasing order, 
repeated according to multiplicity. 

Consider changing the number of the electrons below the Fermi level from N to N' in 
the localization regime. Without loss of generality, we can assume N' > N. We denote 
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by E F and E' F the corresponding two Fermi energies for N and N' electrons, respectively. 
The sum in the right-hand side of (10.19) for N' electrons is written as 



E 



<E' F <E.. 



= E E 

m<N n=N'+l 
N' oo 

+ E E 



(<fm,V x <Pn) (<Pn,Vy<Pm) 
(E m — E n ) 2 

(y m ,^y n ) {<Pn,Vy<Pm) 

(E m E n j 

(ip m , V x ip n ) (if 



x <-> y) 
(x <-> y) 



m=N+l n=N'+l 



= E E 

m<A r n=A r +l 
AT' 

- E E 

m<JV n=W+l 
AT' oo 

+ E E 



(E m — E n ) 2 

y^n _ {A .^ !/) 



(E m — E n ) 2 

(iPm,V x fn) (<Pn,Vy<Pm) 
(E m — E n ) 2 

((fi m ,V x lfi n ) (if 



m=N+l n=N'+l 



(E m - E n y 



-(x+-*y) 
— -(x<-*y) 



(10.20) 



The first double sum in the right-hand side of the second equality leads the Hall conduc- 
tance a xy for N electrons. Therefore it is sufficient to estimate the other two double sums. 
These two sums are compactly written as 

(lfim,V x ip n ) (Vn,VyiPm) 



N' 

E E 

m=N+l n<N and n>N' 



(E m — E n y 



(x <-> y) 



;i0.21) 



In consequence, the difference between the two Hall conductances for N and N' electrons 
is written 

iTic 2 

Aa l Z = -^^Tr AP l A oc [P xA , P yA }AP A °, (10.22) 



xy 



L X Ly 



where AP A oc is the spectral projection onto the localization regime, and 

1 



s,A 



2ni 



/ dzR A (z)v s R A (z), 



(10.23) 



with the resolvent R\(z) = (z — H^^) -1 for the present Hamiltonian H wA on the box 
A = A sys . Here the closed path 7 encircles the energy eigenvalues of the "localized" states. 
In the same way as in Lemma 8.2 and Eq. (8.32), we have 

Aa^ c = ^ J loc (AP{ oc ) + S(L) + O (exp [-/i'L 2 ^ 3 ] ) (10.24) 



with probability larger than (1 — Const. L 2 [ R K+2) 3 l/ 3 ) ) where 



2m. 



x ioc (Apr) = Tr XQ APr[[APr,x], [AP A °,y]] X n, 



(10.25) 



and S(L) is the correction which comes from the boundary region A\Q. In the same way 
as in the proof of Lemma 8.5, we can show that E[5(L)] — > as L — > 00. 
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Lemma 10.2 The following bound is valid: 



E 



X loc (AP A oc ) < Const. AE 1/2 , 



10.26) 



where AE = E' F — E F . 

Proof: Using Schwarz's inequality, one has 
E J loc (AP A oc ) 



2tt i 

< — v/E[Tr Xn&P^Xn] • E[Tr Xn A*AP[° c A Xn ], 



(10.27) 



where we have written A = [[AP A oc , x], [AP A oc ,y]]. From the Wegner estimate (A. 38), one 
has 



\n\ 



Tr xnAP^xn 



Tr X aAP| oc xa 



< Const. AE. 



10.28) 



Therefore it is sufficient to show that the quantity : E[Tr X qA* AP A C A X q] is bounded. 
Schwarz's inequality yields 

Tr Xb (u)A*AP A ° c A Xb (u) = ^Tr Xb (u)A* Xb (v)AP A ° c Xb (w)A Xb (u) 



< J2 V Tr X6(u)l*X6(v)AP| oc X6 (v)l X 6(u) 
v,w 

xy/Tr XbHA*x b MAP^ c x b (w)Ax b (u) 

< Const.^\\ Xb (v)A Xb (u)\\-\\ Xb (w)A Xb (u)\\, (10.29) 



where we have used the bound (8.27) for getting the second inequality. One can show that 
the expectation value of the right-hand side is finite in the same way as in the proof of 
Lemma 8.3. ■ 

We denote by M' A the event M A for the Fermi energy E' F in the proof of Lemma 8.1. 
Note that 



E 



Index(P F *7 a P F ) - ( Ve 2 X I(M A n M' A ) 



< E [|Index(P F f/ a P F ) - T(Pr,nj P )\] 

+ E [|X e (P F ; n, £ P ) - X(P F ; SI, i P ) |] + E [| J(P F ; Q, £ P ) - J(P F , A ; O) \] 



+ E 



J(P F , A ; fi) - ( h/e 2 )o% I(M A n M' A ) 



(10.30) 



From (8.15), (9.28) and Lemmas 9.3 and 9.7, all the terms in the right-hand side become 
small for large L and for a small e. Further, from (8.32), the proof of Lemma 8.5 and 
(10.24), we have 



E 



® xy ® xy 



i(M A nM' A ] 



< E 
+ E 



Aa l ° c \I(M A nM A )} +E 



out' 

xy 



i(M A nM' A ) 



out 

xy 



i(M A n M' A ) 



J loc (AP A oc ) + (small correction). 



;io.3i) 
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From these observations and the fact that the indices are constant for almost every u, we 
obtain 



|Index(P F C/ a P F ) - Index(fW a P F )| E [I(M A n M' A )} 
< E [|Index(P F f/ a P F ) - Index(P F f/ a P F ) | I(M A n M' A )] 



< E 



l loc (AP i 



(small correction) 



< Const. AP 1//2 + (small correction), 



;i0.32) 



where we have used Lemma 10.2 for getting the last inequality. This implies that the index 
must be constant for a small change AE of the Fermi energy in the localization regime 
because the index is equal to an integer for almost every uj. 



11 Widths of the Hall conductance plateaus 

In this section, we prove that the widths of the Hall conductance plateaus are of the bulk 
order under certain conditions for the potentials, by estimating the number of the localized 
states. The conditions are realized for weak potentials as we will see in this section. 

Consider first the case with Ap = 0. To begin with, we note that, when the strength of 
the random potential V u continuously increases from A = to A G [— A_, A + ] C [A min , A max ], 
the energies E of the n + 1-th Landau band are broadened into the interval, 

-Halloo - A_ Ml < E - S n < HVIU + X +Ul . (11.1) 

Let 5 be a small positive parameter. For the lower region of the band, we choose 
A + = A'° w < A max , A_ = A'° w and 5_ = 5 in the condition (5.36) so that the pair 
(A+,A_) = (A!£ w , A'° w ) satisfies the condition (2.7) with a small A_ = A 10 ". Then the 
condition (5.36) for the energy E leading to a localized state becomes 

Sn-i + HK+IU + \ l ° w Ul + 6hw c <E<S n - Halloo - A!^! - AS. (11.2) 

We call this interval the lower localization regime. For the upper region of the band, we 
choose A + = A+ p , A_ = A" 15 > — A min and 5+ = 5/2 in the condition (5.37) so that the pair 
(A+,A_) = (A+p.A" 15 ) satisfies the condition (2.7) with a small A+ = A+ p . The condition 
(5.37) for localization is 

£n + HK+H + A^ui + AS < E < S n+1 - ||V^-||oo - \ up Ul - 5hu c . (11.3) 

We call this interval the upper localization regime. We require that the positive constants, 
6, \ l ° w and A" p , satisfy 

\\V + \\oo + ||^o"lloo + (A+ w + A^ p ) Ul + 26hu c < hcu c (11.4) 

so that the lower and upper localization regimes overlap with each other. This condition 
is satisfied for a large strength B of the magnetic field for fixed strengths of the poten- 
tials. We stress that we allow the possibility that the spectral gap between two neigh- 
boring disordered-broadened Landau bands vanishes owing to the tails A G [A m i n , — A" p ] U 
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[\ l ° w , A max ] of the coupling constants. In this situation, all the states in the n + 1-th broad- 
ened Landau band are localized except for the energies E satisfying —5E- < E—£ n < SE + , 
where 5E_ = \\V ~\\ oo + A'° w «i + AS and 5E + = ||V r + || oo + A+ p ui + A£. In other words, the 
number of the extended states can be bounded by the number of the energy eigenvalues 
E satisfying -5E_ < E - S n < 5E + . 

In order to obtain the upper bound for the number of the extended states, consider 
first the special case with Vq = in the Hamiltonian Hq of (2.2). Namely, the Hamiltonian 
Hq is equal to the simplest Landau Hamiltonian H^ of (3.1), combining with the present 
assumption Ap = 0. In this case, the number of the extended states can be estimated 
with probability nearly equal to one for the sufficiently large volume 10 by using the Wegner 
estimate (A. 38). When V ^ 0, the deviation of the energy eigenvalues is bounded from 
above by HV^Hoo and from below by — 1 1 Vq - 1 1 oo - From these observations and the min-max 
principle, we can estimate the number iV ex t of the extended states which appear only near 
the center of the band as 

iVext < C$Ki 0) \\g\\ oo (6E + + 8E_ + HVo+IU + ||y "IU)|A| (11.5) 

in the case with Vq ^ 0, where C$ and are the positive constants for the case of 
Vq = 0. Since the total number of the states in the n + 1-th Landau band is given by 
M = \A\eB/(2nh), the number N\ oc of the localized states in the band is evaluated as 

/ ^(0) \ 

N loc — M — N cxt > B\A\ ^ - dV-^WgW^SE] (11.6) 

with 5E = 2 (Halloo + HVlloo) + (A 10 ™ + A+ P )«i + 2A£. We note that K^/B ~ Const, 
for a large B from the remark below Theorem A. 2. Thus, if the strength of the potential 
Vo is sufficiently weak, we can choose the parameters X l ° w , A!J_ P , AS so that the right-hand 
side (11.6) is strictly positive for any large magnetic field. This implies that the number 
Aioc is of order of the bulk. In order to discuss the case for a strong random potential 
which behaves like II 

^||oo ~ B for a strong magnetic field, we recall u\ — ^ 1 1 1 1 oo - 

We also 

have Kf ] = 0(1) which was already obtained at the end of Section 5. From these and the 
same argument, we can also get the lower bound for the number of the localized states, 
i.e., the width of the Hall conductance. 

Let us see that the above estimate for the widths of the plateaus gives the optimal 
value in the limit B j oo for V = 0. From the above bounds, we have 

^ < Const. [(A 10 " + A^ p ) Ml + 2A£] . (11.7) 

From the argument about the initial decay estimate for the resolvent in Section 5, we can 
take the three parameters, A'° w , A" p , AS, so as to go to zero in the strong magnetic field 
limit B | oo. Thus the density of the extended states in the Landau level is vanishing in 
the limit. 

Next consider the case with Ap ^ 0. The method to show the existence of the Hall 
conductance plateau with the width of bulk order is basically the same as in the above 



10 



See, for example, Chap. VI of the book [28]. 
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case with A P = 0, except for considering strong potentials. We assume that the bump u 
of the random potential V u is written u = hu c u with a fixed, dimensionless function u, and 
that the vector potential Ap satisfies HlAplH^ < uqB 1 / 2 with a small, positive constant 
ao- Instead of the condition (11.4), we require that the corresponding positive constants, 
6, \ l ° w and A" p , satisfy 

IIVIU + || VIU + (AT + A^H + 5Tiw c 
+ ■^|||Ap||| 00 (v^+ v ^) + -^|||Ap|||^ <^ c . (11.8) 



Similarly, for the lower localization regime, we choose A + = A 10 ™ < A max , A_ = \ lmv and 
$+ = $- = 5/2 in the condition (5.47) so that the pair (A+, A_) = (X l ° w , X l ° w ) satisfies the 
condition (2.7) with a small A_ = A'° w . Then all the states in the n + 1-th broadened 
Landau band with the energies E < S n — 5E_ are localized, where 

.low, . If. , v/2e 



5E_ = H^Tlloo + A'° w Ml + -6hw c + ^—11 1 ApHUy^. (11.9) 

2 ,/mT v 



For the upper localization regime, we choose A+ = A+ p , A_ = A 1 ? > — A m i n and 8+ = 
5- = 5/2 in the same condition (5.47) so that the pair (A + , A_) = (A+ p , A" p ) satisfies the 
condition (2.7) with a small A + = A" 15 . Then all the states in the n + 1-th broadened 
Landau band with the energies E > S n + 5E + are localized, where 

V2e nr e 2 



5E + = HVIU + \f Ul + ^-5huo c + -^HIApIIU^ + /HHApHIL (11.10) 

2 x /m P v lm P 



The corresponding 5E in (11.6) is given by 



5E = 2 (HFo+IU + Halloo) + (XT + A lm >i + 6hu c + -^|||A P ||| 00 ^+ — |||A P ||^. 

(11-11) 

Consequently there appears the Hall conductance plateau with the width of the bulk 
order for a fixed potential Vo, for a strong magnetic field, and for small parameters, 
A^ p ,Ai_ ow ,<W 



12 Corrections to the Linear response formula 

The aim of this section is to prove that both of the acceleration coefficients ^y U y, u — %-,y, 
in the linear response formula (8.1) are vanishing in the infinite volume limit, and that 
the corrections 5a uy (t), u — x,y, due to the initial adiabatic process in (8.1) satisfy the 
bound (2.27). 

We recall the expression of the acceleration coefficients [8], 



luy 



L X Ly 



N 



m, 



5 u , y + Tr v u (P y , A PF,A + Pf,a P va) 



for u — x,y. 



[12.1) 
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Using the partition of unity, {xb( u )}u, which was introduced in Section 8, we have 

Tr v u (P yA P FA + P F , A P y ,A) = ]T [Tr v u P yAXb (u)P FA + Tr ^Pf.aX&^P^a] . (12.2) 

u 

First let us consider the first term Trv u P yy \Xb( u ) Pf,a in the summand in the right-hand 
side. Since we can shift the location of the box s&(u) by using the magnetic translations, 
we can assume dist(u, <9A) = 0(L). Write R A (z) = (z — P^a) -1 - Note that 

Tr v u P y! AXb(u)P F ,A = 77— f dzTr v u R A (z)v y R A (z)xb(u)P F , A 

= -^—j dzTrv u R A (z)v y x S A RA(z)Xb(u)P FA 

+ ~hi I dzTl v M z ) v v^ ~ Xa)Pa(^)X6(u)P FiA , (12.3) 

where Xa l $ the C 2 , positive cutoff function which was also introduced in Section 8. By the 
same argument as in the proof of Lemma 8.2, the absolute value of the second term in the 
right-hand side has a stretched exponentially decaying bound as in Lemma 8.1. Since the 
number of u for the summation in the right-hand side of (12.2) is of order of the volume 
L x L y = 0(L 2 ), the corresponding contribution is vanishing in the infinite volume limit 
L f oo. Using the identity, v y x A = (i/h)[y, z — P W) a]aa> one nas 

^- J dzR A (z)v y x A R A {z)xb{u) 
= ^[P F , A ,yx A }Xb(u) - ^ ^dzR A (z)yW( X S A )R A (z)xb(u) (12.4) 

for the first term in the right-hand side of (12.3). This second term in the right-hand side 
also gives a small correction. In consequence, only the first term in the right-hand side of 
(12.4) may lead to a nonvanishing contribution in the infinite- volume limit. 

Since the second term in the summand in the right-hand side of (12.2) can handled in 

the same way, we get 

Tr [v u P y , A Xb(u)P F , A + v u P FA Xb{u)P y , A ] 
= ^Tr [v u [P F)A , yx A }Xb(u)P F)A + v u P F)A Xb(u)[P F)A , yx A }\ + corrections 

= ^-Tr {-v u yx{PF, A Xb(u)P F , A + v u P FA Xb(u)P F , A yx{} + corrections 

= — —5 Uty Tr Pf,aX&(u)P f ,a + corrections, (12.5) 

where we have used v u y = — {ih/m e )5 u ^ y + yv u for getting the last equality. Substituting 
this and (12.2) into the expression (12.1) of 7 UJ/ , we obtain 

lim 7™ = for u = x,y with probability one. (12.6) 

A.sys| R 2 
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Next consider the corrections 5cr uy (t) due to the initial adiabatic process in (8.1). We 
begin with recalling the expression [8], 



6* uy (t) = (*§>,v™[l " WW - <A>r*o,A) + cc, (12.7) 



7.W - ffW^,W<fiW 



with 



iT 



h 



e- r > T e l£T/h + 



ft 



h 



s 2 (s + thny 



jEt/h 



'12.t 



where <I?oa * s the N electron ground state vector with the energy eigenvalue Eqa \ v ^ the 
N electron velocity operator, and Pq2 the projection onto the N electron ground state. 



Since all the contributions can be handled in the same way, we consider 



^0,A ) u u - r cx 



L X Ly 



±j x ±j v \ L ' 'J 

(12.9) 

with = (E^2 — H^l)t/% as an example. Here we have written P cx = 1 — pfff for short. 
In order to eliminate the factor e ld , we use Schwarz's inequality. As a result, we obtain 

Wuy{t)\ 2 <N u {t)N y {t) With 



N u {t) := 

and 
M y {t) : = 



L x L y 



+ 



4 JV) -<A ) +^ 



y u V 0,A 



(12.10) 



^0,A ) u j/ -'ex 



L X Ly 



p (JV) „(JV) 
-^O.A _ ^o;,A 



p (AT) „(jV) fi 



u 2/ V 0,A / • 



Further the application of the inequality \l ah < (a + 6)/2 for a, 6 > yields 



uy\ 



< 



V s / 4 Af u (t) + r s/ *tiy{t)\ /2 for s G (0, 1/3). 



(12.11) 
(12.12) 



Since the present system has no electron-electron interaction, N y (t) is written as 



X- 



z - if w>A z - if w A + iftr? 



L x L y a 27T2 J 7 2; — I 



1 



J y v 



1 



^ - -£L,a ^ - -£L,a + ^ 



1 



z' - H^a z' - Hu,a - ihfj 



y z'-H L 



12.13) 
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Here we have introduced the partition of unity {xf>(a)}. In the same way as the above, 
the summand can be further written as 

X Xs(a)^ / dz'- -P- —1(2'-^)-',!/]+ correction. (12.14) 

ZTCZ J~f Z — ii w A — liXTj 

Here the correction vanishes almost surely in the infinite volume limit by taking a suitable 
sequence {L x , L y } of the system sizes as in the proof of Theorem 8.5. If z of the resolvent 
R\(z) — (z — H^^aJ^ 1 is not near to the spectrum <j(H Wi a), then the resolvent is bounded 
and decays exponentially at large distance. Therefore we consider only the contributions 
from the paths near the Fermi energy E F in the first term in the right-hand side of (12.14). 
As a typical example of the corresponding contributions, let us consider 

4 := E r dt dt '^Xb(vL)P FA x b (v)iy,RA(E F + it)]x b M 

u,v,w,z Jt - 

x r]R A (E F + i(t + hr]))xb(a.)r]R A (E F + i(t' - hr})) Xb (z)[R A (E F + it') , y] Xb (u) . 

(12.15) 

Note that 



|Tr Xb (u)P F ,AX 6 (v)A| < x /Tr X6 (u)P FiAX6 (u)^TrA* Xb (v)P F)AX6 (v)A 

< Const. \\A\\ (12.16) 

for any bounded operator A, where we have used the bound (8.27), and that 

\\x b (v)[y,R A (z)} Xb (w)\\ < (Const. + \v 2 -w 2 \)\\ Xb (v)R A (z) Xb (w) || , (12.17) 

where we have used the decomposition y = y — y b + y b in the proof of Lemma 8.3. From 
these observations, we have 

< Const. E (Const. + \v 2 — u> 2 |)(Const. + \z 2 — u 2 \) 

u,v,w,z 

x f t + dt ||x 6 (v)i2A(^F + iOX6(w)||-77||x6(w)i2 A (^F + i(t + ^))X6(a)|| 
x f t + dt' 77 ||x6(a)P A (P F + i(t' - h V )) Xb (z)\\ \\ Xb (z)R A (E F + it') Xb (u)\\ . 



/ dt--- = dt--- + dt---. (12.19) 

Jt- Jt- J-%n/2 



(12.18) 

The first integral is decomposed into two parts as 

dt---+ I 

-hri/2 

The second part of the integral is estimated as 

f t+ dt \\x b {^)RA{E F + tt)x b M\\-n\\x b {w)RA(E F + t{t + hn)) Xb {^ 

< Const, x rj s / 4 [ t+ dt \\ Xb (v)R A (E F + it)xb(vr)\\ 8/A l^ /4_1 

J-hri/2 

x || X6 (w) J R A (P F + z(t + ^)) Xb (a)|r /4 , (12.20) 
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where we have used the inequality ||xb(w)i? A (-^F + i{t + hr))) Xb (a) || < 2/(^77) for t > 
—Tin/ 2. For the first part of the integral, we can obtain a similar bound by using 
\\Xb{v) R\(E-p + it)xb(w) || < 2/(hn) for t < —Tin/2. Clearly the second integral about 
t' in the right-hand side of (12.18) can be treated in the same way. Combining these 
observations with Holedr inequality 

E [|| Xb (v)it> A (£ F + ^)x6(w)|| s/4 \\ Xb (w)R A (E F + i{t + ^)) X6 (a)|| s/4 
x \\ Xb (a.)R A (E F + i{t' - hn)) Xb (z)\\ s/ * \\x^)Ra(E f + it') Xb i^)\\ s, \ 
< E [\\ Xb (v)R A (E F + it) Xb (w)\\ s ] 1/4 E [\\ Xb (w)R A (E F + i{t + hr))) Xb (a 



x E 



\ Xb (a.)R A (E F + i(t' -%v))Xb(z) 



1/4 



E 



\ Xb (z)R A iE F + it') Xb (u) 



1/4 
1/4 



;i2.2i) 



we obtain 



E|/ a | < Const, x n s/2 (Const. + \v 2 - w 2 |)(Const. + \z 2 - u 2 \) 

u,v,w,z 

Xe -M|v-w|/4 e -^|w-a|/4 e -/i|a-z|/4 e -/u|z-u|/4 

< Const, x n s/2 , 



(12.22) 



where we have used Fatou's lemma, Fubini-Tonelli theorem and the fractional moment 
bound (7.1) as in the proof of Lemma 7.1. Combining this bound, (12.13), (12.14) and 
(12.15), the expectation of Af y (t) of (12.12) is bounded by Const, x n s/2 . 
Since the bound, 



S- 1 + iS + ihr})- 1 



<4£" 



(12.23) 



holds for SeH, the expectation of Af u (t) of (12.12) can be proved to be bounded in a 
easier way. As a result, we obtain that there exists a sequence {L Xtn , L y ^ n } n of the system 
sizes such that the bound, 



lim N uy {t) < Const.?f /4 , 



(12.24) 



holds almost surely, where the positive constant may depend on uo. 

Since the rest of the contributions for Sa uy it) can be handled in the same way, we 
obtain the desired result that the bound, 



\Sa uy (t)\ < IdH + ^HTle^ + CaH^ 4 , 



(12.25) 



holds almost surely for s G (0, 1/3) and for u — x,y. Here the positive constants, Cjiu) < 
00, j = 1, 2, 3, may depend on 00. Choosing s = 4/13 < 1/3, we get the bound (2.27). 



A Wegner estimate 

In this appendix, we study the Wegner estimate [35] for the density of states for a single- 
electron Hamiltonian in a general setting. For this purpose, we follow the argument by 



Widths of the Hall Conductance Plateaus 



52 



Barbaroux, Combes and Hislop [15] . However, the following argument is slightly simplified 
compared to their original one because we are interested in two dimensions only. 

Consider a single spinless electron system with a disorder potential V u in d dimensions. 
The Hamiltonian is given by H u = Hq + V^ on L 2 (R d ) with the unperturbed Hamiltonian, 

H = ^-(p + eA) 2 + V . (A.l) 

We assume A G C 1 (R d ,R d ) and V G L°°(R d ), and assume that the Hamiltonian H is 
essentially self-adjoint with a boundary condition. As a disorder potential V u , we consider 
an Anderson type potential of impurities, 

V u (r)= E A a H«(r-a), (A.2) 

aeL d 

for r = (xi, x%, ■ ■ ■ , x d ) G H d . The constants {A a (cj)| a G L d } form a family of independent, 
identically distributed random variables on a rf-dimensional periodic lattice L d C R d . 
The common distribution has a density g > which has compact support and satisfies 
g G L°°(R) fl C(R). We assume that the single-site potential u is non-negative and has 
compact support. We write the sum of the single-site potentials u as 

U{r) : = ]T «(r - a) for r G R d . (A.3) 

aeL d 

We assume ||W||oo < +oo. Clearly this implies ||m||oo < +oo. From these assumptions, 
we have ||K,||oo < < +oo with some positive constant i>r which is independent of the 
random variables. 

For a bounded region A C R d , we denote by H^^a = H 0j \ + V^a the Hamiltonian H w 
restricted to A with a boundary condition. Here V^a = K;|A, i.e., is the restriction 
of Vu to A. We assume 

U A := W|A > W mi „XA, (A.4) 

where ZY m i n is a positive constant which is independent of the bounded region A, and xa is 
the characteristic function for A. Namely there is no flat potential region satisfying u — 0. 
Further we assume that 

Tr (# ,a + £ mi „r 2 < K \n\ n " (A.5) 

for a finite region Q C A, with a positive constant £ m i n > HV^Hoo, where Tr stands for 
the trace on L 2 (A); K and no are the positive constants which are independent of the 
volumes |A| , \Q\ of the finite regions A, il. If the vector potential A satisfies the additional 
assumption A G C 2 (R d , R d ), then the inequality (A.5) is valid in the dimensions d < 3. 
See ref. [15] for details and also for the treatment in the case of higher dimensions d > 4 
in which they require a stronger assumption than the above assumption (A.5). Consider 
the present system with the unperturbed Hamiltonian (2.2). From the inequality (3.9) for 
the lower edge of the Landau band, we have 



TY(#o,A + £min) 2 Xn < Const, x \b\/B 



(A.6) 
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for strong magnetic field strengths B. Here we have assumed that the region VL is contained 
in a rectangular box b such that the volume \b\ of the box satisfies the flux quantization 
condition \b\/(27v£ B ) G N. Namely 

n — 1 and K ~ Const, x B^ 1 for a large B. (A. 7) 

Let <5cj,a(A) denote the spectral projection for H u \ with an energy interval A C R. 
Let ipE be an eigenvector of the Hamiltonian H^a, i-e., H^^e = EipE with an energy 
eigenvalue £eR. The Schrodinger equation is written as 

(# ,A + ^rnin) i=E = (-K,,A + £ + £min) ^- (A.8) 

Since if ,A + £min > from the assumption £ min > 1 1 Vq - 1 1 co j one nas 

V'-E = TT 77; (~K;,A + E + £ min ) ^£ = — (-K,,A + # W ,A + £ min ) V'E- 

-"0,A + ^min -"0,A + ^min 

(A.9) 

Using this identity, one obtains 
Tr Q WiA (A) 

= Tr (i?0,A + £min) 2 ( — K»,A + -#w,A + £min) <3uj,a(A) ( — V^A + Huj,A + £min) 
2 



Tr (# ,A + ^minK K, iA Q Wi a(A)K >)A 



— Tr (if ,A + ^min) 2 K),A (#u>,A + £min) <?w,a(A) 



— Tr (i?Q,A + £min) 2 (-f^j,A + £min) <3u;,a( A) V^A 



+ Tr(# ,A + £min) (#u,,A + £min) 0,,a(A). (A.10) 

Let us evaluate the first term in the last line of the right-hand side. Substituting the 
expression (A. 2) of into the term, one has 

Tr (# 0> A + £minr 2 V uA Q uA (A)V uA 
= ^aA b Tr (Hq A + ^min)" 2 U &A Q ujA ( A)w bj A 



a,b 

^ A a A b Tr (#0,A + ^min)~ 2 uj^ uJ%Q uA (&)v% A , (A.ll) 

a,b 



where we have written u a (r) = tt(r — a) and u &A — u & \A. Since the operator 

Tgi := <a (#o,A + £ mi „r 2 <1 (A.12) 

is compact, there exist a pair of orthonormal bases, {(p£ ^}^Li and {^ 1 - ) }^ =1 , and nonneg- 
ative numbers {/4^}$£Li such that 11 

CO 

T&=E^V?WV")- (A.13) 

n=l 



11 See, for example, Chapter VI of the book by Reed and Simon [42]. 
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The numbers 11$ are the eigenvalues of |T^ a |. Using this representation (A. 13), one has 
|TrT«^iQ w , A (AK(i 



n=l 
■y oo 



Z n=l 



(1) .1/2 
n ) "a,A 



. (A.14) 



Therefore the expectation value of the left-hand side of the first inequality can be bounded 
from above as 



Ea [I^TW uigQ WiA (A)«i 



/2 



1 



(A.15) 



where Ea[- • •] stands for the expectation with respect to the random variables on a region 
A C R d . The right-hand side can be evaluated by using the following Lemma A.l which 
is essentially due to Kotani and Simon [43]. In order to make this paper self-contained, 
we give the proof of Lemma A.l in Appendix B, following from ref. [37]. 



Lemma A.l Let v be a nonnegative function satisfying v < u a \. Then 



[ dA a <?(A a y/2Q WiA(A y/2 



< IMIco|A|. 



From the bound (A.16) and the inequality (A.15), one has 



E, 



1/2 



< 



T (l) 
1 b,a 



(A.16) 



(A.17) 



where 



|i := Tr [• ■ -|. Moreover, combining this bound with (A. 11), one gets 



E, 



Tr (H 0A + £ min ) 2 V Ui aQ Ui a(A)V uA 



= E Ea [A a A b Tr T« A)u^" 

a,b 

< M 2 £ E A [|Tr ^5Q.,a(A)<1 

(A.18) 



a,b 



< -M 2 |klloo|A|5:||T£ 

a,b 

where .M := sup Ae suppg |A|. 

Next consider the second term in the right-hand side of the second equality in (A. 10). 
From the assumption (A. 4), one can easily find a set of nonnegative functions {wa,A}a 
satisfying the following two conditions: 



u a ,A < Ua,A for any lattice site a, 



(A.19) 



Widths of the Hall Conductance Plateaus 55 
and 

E ^a,A = WminXA- (A. 20) 

a 

Using the identity (A. 20), one has 

Tr (H 0t \ + £ m in) 2 V Uj A.(H Uj A + £min)<5cj,A(A) 
= ^ A a Tr (iJ ,A + £ min )~ 2 M a ,A(^,A + ^min)<5a;,A(A)xA 



1 



^]A a Tr (# ,A + £min) Ma,A(^,A + ^mm)<5^,A(A)M b ,A 



= 7-7- E A aTr Tg a ^(^, A + £: min )g a; , A (A)45, (A-21) 

where 

TgL := <1 (#o,A + ^ mi „)- 2 u l a % (A.22) 

In the same way, 

|Tr u VJ(ff WiA + £^ n )Q WiA (A)fii2 

< E I W 2) , uS(^,a + £ m in)Q.,A(A)<i^ 2 )) 

n=l 

< E /4> 2) [ W 2) , <1^,a(a)<M 2) ) + K 2) <iQ.,A(A)<i^))] , 

Z 71=1 

(A.23) 

where {y^j-^Li and {V^ln^i are orthonormal bases such that 

00 

T&=E^?W\"0 (A.24) 

71=1 

with the eigenvalues /i^ 2 ) of |T b 2 a |, and 

£ max (A) = sup|£ + £: min |. (A.25) 
EeA 

Combining (A. 21), (A.23) and Lemma A.l with the condition (A. 19), one has 



E, 



Tr (H 0A + £ min )- 2 V^ A (H^ A + £ min )Q^ A (A)\] < M £ ^ A) \\g\U Aj ^ ||Tg 



a,b 



In the same way, 



(2) 
>,a 

(A.26) 



E; 



Tr (^o,A + £ m i„)" 2 (^,A + £ m in)Q.,A(A)K,, A l <M S ™; {A) \\g\UA\J2\\r h „ 



(3) 
b,a 

(A.27) 
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and 



E, 



Tr (#0,A + £min) 2 (^,A + ^min) 2 <5aj,A(A) 



< 



^max(^) 

1Arr\\r\ 



\9\ 



(4) 



a,b 



where 



~(3) 1/2 /rj . c \-2 -1/2 , ~,(4) -1/2 /tt . c \-2 -1/2 

T b,a : = <AW + M <A alld T b,a := <A (#0,A + £min) <A- 



(A.28) 



(A.29) 



Next let us estimate X) a b 



T 



(0 

b,a 



, i = 1,2, 3, 4. For simplicity, we consider only the case 
with i = l because all the other cases can be estimated in the same way. We decompose 
the sum into two parts as 



£ll T 



a) 

b,a 



E Ki , + E 



a,b 



a,b 

overlap 



a,b 
non — overlap 



T 



(i) 

b,a 



(A.30) 



where the first sum is over the lattice sites a, b such that the corresponding two potentials 
«a,A, Wb,A overlap with each other, and the second sum is over those for the non-overlapping 
potentials. Note that 

(A.31) 



T 



CO 



b,a 



< 



1 b,b 



T 



(0 



where we have used the inequality < || A|| 2 ||i?||2 for bounded operators A,B. Here 

the norm || ■ • • || 2 is defined as 1 1 ^4 1 1 2 := \^TrA*A for a bounded operator A if the right-hand 
side exists. Using the inequality (A.31) and the assumption (A. 5), one has 



E | T bi 1 < Const, x XolklU Isupp u\ no |A| , 



(A.32) 



a,b 

overlap 



where the positive constant depends only on the lattice L d and on the support of the 
potential u, and so the constant is finite from the assumptions on the lattice and the 
potential. 

Using Proposition C.l in Appendix C, the second sum in the right-hand side of (A.30) 
can be evaluated as 



E 



a,b 

non — overlap 



T 



(i) 

b,a 



< 



< 



Const, x \\u\ 



(1 ft 2 )(£min 11^0 I loo 
Const. X \\u\\nn 



K |supp u\ no/2 Yl rV 



a,b 

non — overlap 



\V f 



■AVG(a)|supp M| no/2 |A|, (A.33) 



lloo; 



where a is given by (C.2), r = dist(supp u a , supp Ub), and the positive constant K 2 (a) 
satisfies the bound: 

K 2 (a) < Const. xY,r n e- ar (A.34) 

b 

with the positive constant and with a fixed lattice site a. Combining this result with the 
bound (A.32), one has 



E T b,a 1 < Const, x AVfilMUsupp M| no/2 |A| 

a,b 



(A.35) 
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with the positive constant 

Kl = |supp u^ 2 + f {Q) (A.36) 

1^1 tv )\O mm \\Vq | |oo ) 

Clearly this constant K\ is independent of the strength B of the magnetic field and of the 
random potential V^. From the condition (A. 19) for m &i a, the following bounds also hold: 

2 l^af < Const - x ^o^i||M||oo|supp u\ no \A\ (A.37) 

a,b 

for all i = 1,2, 3, 4. Combining this, (A.10), (A.18), (A.26), (A.27) and (A.28), one obtains 
the following theorem: 

Theorem A. 2 Assume the conditions (A. 4) and (A. 5). Let A = [E — SE, E + SE] be an 

interval of the energy with SE > 0. Then the following inequality is valid: 

Prob [dist((T(^ WiA ), E) < 5E] < E A [Tr Q W)A (A)] < C w ^ 3 ||ff||oo^|A| (A.38) 

TOi/i some positive constant Cw wra^/i 

K 3 = [M + £ max (A)/U TDin ] 2 K K 1 \\u\\ oo \supp u| n °. (A.39) 

Here E A [- • •] £/ie expectation with respect to the random variables on a region A C R '. 

Remark: Consider the present system with the unperturbed Hamiltonian (2.2). The 
positive number K behaves as K ~ Const, x B^ 1 for a large B as in (A. 7). From the 
definition (A. 25), £ max (A) ~ Const, x B for a large B. Substituting these into the above 
expression (A.39) of K 3 , one has K 3 ~ Const, x B for a large B. From this observation 
and the above theorem, one notices that the upper bound for the number of the states in 
the energy interval with a fixed width SE is proportional to the strength B of the magnetic 
field for a large B. 



B Proof of Lemma A.l 



Following Combes and Hislop [37], we give the proof of the Kotani-Simon lemma [43]. 
We begin with preparing the following lemma: 

Lemma B.l Write H^^a = H^ A + A a w aiA . Let v be a nonnegative function satisfying 
«a, a > v, and define 



K(X, E-iS) = v 1 ' 2 (H' uA + Xu aA -E + i6) ^ v 1 ' 2 
for E G R and S > 0. Then 

X 2 



L dX vTxl K ^ E -^ 



< IT 



(B.l) 



(B.2) 



for A > 0. 
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Proof: Since K(X, E — iS) is holomorphic in A in the upper half-plane, one has 



L 



An 



A 2 + A 2 



K(X, E-iS)= ir\ K(i\ Q , E-iS). 



(B.3) 



Note that 



-ImX^Ao, E - iS) 

> 1/2 (K A -i\ oUaA 



-E-i6 



-1 



(Ao«a,A + 5) (H' uA + i\ u aA -E + iS) 1 v 1/2 



> \ K(i\ ,E -i5)*K(i\ ,E -iS), 



(B.4) 



where we have used the assumptions u a jy > v and 5 > 0. This implies HfT^An, E — i8)\\ < 
Xq 1 . Combining this with the above (B.3), one has the desired result (B.2). ■ 



Proof of Lemma A.l: Let Ad A and A. Using Stone's formula, one has 

(<p,v 1/2 Q uA ( A ) vl/2( p) < -limlm I dE U,v 1 / 2 (H„ >A - E + iS)' 1 v^tp) 
for any vector <p. Further, 

\ 2 



(B.5) 



< — lim Im / dE 

7T 8io 



\ 2 

d\ 

R A 2 + Aq 



(</?, X(A a , £ - i<$)p) < |A||M 



:b.6) 



by using Fubini's theorem and Lemma B.l. Here An > 0, and K(X,E — i8) is given by 
(B.l). Since g G Loo(R) with compact support from the assumption, one has 



R 



d\ a g{K)v l l 2 Q w ^)v 1 ' 2 
for any A > 0. This proves the bound (A. 16). 



A 2 + A 2 



< sup#(A) 

A /v 



An 2 



(B.7) 



C A decay estimate of T^ a 

In this appendix, we follow Barbaroux, Combes and Hislop [15], in order to estimate T^ a 
that appear in Appendix A. The result is summarized as the following proposition: 

Proposition C.l Let v,w be bounded functions with a compact support, and suppose 
dist(supp v, supp w) — r with a positive distance r. Then 

Const. X IHUHU T , f\ mm „ino/2 M , _ ,n /2\ r n-ar 



v(Ro,a) 2 ™ < " Ullb 2 ' " Ko (|supp vr' 2 + |supp wr' 2 ) r" e — (C.l 

1 ~ K Jv^min ~~ || Vq Woo) 



with 



\\V Woo) 

a = — k with k £ (0, 1) (C2) 

Oil 



and with some positive number n, where Ro,a = (-^o,a + £min) 1 , and || • • • ||i := Tr| 
The constants K and no are given in (A. 5). 
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Let O be a region in R d . Then we denote by dfl the boundary of fi, and define the 
subset f2i n of Q as f2; n = {r G fl| dist(r,<9f2) > 5} with a positive 5. With a small 5, one 
can take three regions i = 1, 2, 3 satisfying the following conditions: 



supp v C fi£> C fiW C n£> C f2< 2 > C fig? C fi (3) C A, 



< Const, x r d , 



dist(supp w,f2 {3) ) > r/3, 



and 



dist(supp w,r (1) ) > r/3, 

where T« = and we also write r« = for i = 2, 3. 

Let Xi £ C 2 (A), i = 1,2,3, be three nonnegative functions satisfying 



XiUo = 1 and Xi 



A\n« 



0. 



(C.3) 

(C.4) 
(C.5) 

(C.6) 



(C.7) 



In the following, we denote by x\ the characteristic function Xr« °f the region rW for 
i = 1,2,3, and write i? ,i = -Ro,n« — (-ffo.fiW + ^min)" 1 for i = 1,2,3. Next introduce the 
geometric resolvent equation, 



XiRo,A = Ro,iXi + Ro,iW(xi)R 



0,A 



for i — 1,2, 3, where W(- ■ •) is given by (5.7). 
Using (C.8) and v\2 = v, one has 

V (Ro,A) 2 W = V X2 (Ro,A) 2 W = vR ,2X2Ro,A W + vRq, 2 W(x 2 ) (#0,a) 2 W. 

The first term in the right-hand side can be rewritten as 



(C.8) 



(C.9) 



>Ro,2X2Ro,A W = V ( R 0,2) 2 W(X2)R0,AW 

= V (Ro,2) 2 W(x 2 )X3Ro,A W 

= v(R ,2) 2 W(x 2 )R 0i3 W(x 3 )R 0A w 



(CIO) 



by using X2~w = X3 W — 0, W(x 2 )x 3 — W(x2) and the geometric resolvent equation {CI 
Therefore 



\vR ,2X2Ro,A W \\l 



< 
< 

< 



v (Ro,2) 2 W(x 2 )R , 3 W(x3)Ro,a™ 
\vR 0:2 \\ 2 \\Ro,2W(x2)Ro,3W(x 3 )R 0jA w\\ 2 
\vR , 2 \\ 2 \\wR ,aW*(x3)Ro,3W*(x 2 )R , 2 \\ 2 

\VR ,2\\ 2 \\WR 0A X' 3 W*(X3)R0,3W*(X2)R0,2\\2 

\VR , 2 \\ 2 \\™Ro,AX3\\ \\W*(X 3 )R0,3W*(X2)R0,2\\2 1 ( C - U ) 



where we have used X3W*(x 3 ) = W*(x 3 ), the equality ||A*|| 2 
||AB||i < ||A|| 2 ||-B||2 for bounded operators A, B. The norm || 
\ / TrA*A for a bounded operator A. 



2 and the inequality 
2 is defined as \\A\\ 2 : = 
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The second term in the right-hand side of (C.9) can be written as 

vR 0:2 W(x2) (Ro,a) 2 w = vjdRo^W^) (R ,a) 2 w 

= vR ,iW(xi)Ro,2W(x2)(Ro,a) 2 w 

by using vxi = v, XiW(x2) — 0, and the geometric resolvent equation 

xiRo,2 = Ro,iXi + Ro,iW(xi)Ro,2- 

Therefore the norm can be evaluated as 



(C.12) 
(C.13) 



vRo,2W(x2) (Ro,a) 2 w 



vR ,iW(xi)Ro,2W(x2) (Ro,a) 2 ™ 



vRo,ix'iW(xi)R ,2W(x2) (Ro,a) 2 ™ 

< \\vR 0A x[\\ \\w(xi)R ,2W{x 2 ) (R ,a) 2 w 

< WvRo^x'A ||W(xi)i^2W(x 2 )i2o^|| 2 ||i2o^«;|| 2 ,(C.14) 

where we have used the identity XjW(Xi) — W(Xi)- 
Combining (C.9), (C.ll) and (C.14), one has 

^(.Ro.a) 2 ^ i < \\vR 0>2 x 2 Ro,a w \\i+ vR 0>2 W(x 2 ) (R ,a) 2 w 1 

< \\vRo,2h IMo,aX 3 II \\W*(x 3 )R , 3 W*(x2)Ro,2\\2 

+ II^A^Ibll^ix'ill \\W(xi)R , 2 W(x 2 )Ro,a\\ 2 - (C.15) 

In order to estimate the right-hand side, we use the following lemma: 

Lemma C.2 Let ip be a vector in the domain of the Hamiltonian H Q (2). Then 

\\W*(x3)Ro,sW*(xM\ < Const, x \\<p\\, (C.16) 

where the positive constant in the right-hand side depends only on the cut-off functions, 
X 2 and xs- 



Proof: Note that 

\\W*(x 3 )R , 3 W*(xM\ < ||w*(x3)<3 2 || |< 3 V*(x 2 V 

Using (5.9), one has 



W*{xz)Rl{ 3 



< 



h 



(Vxs) ' (P + eA)< 3 2 



+ 



ti 2 



2m e 



1/2 



(C.17) 



(C.18) 



The first term in the right-hand side can be estimated as follows: Using the Schwarz 
inequality, one has 

< 3 2 (p + eA) • (VxsXVxs) • (p + eA)< 3 V) 



< 



x 



,<3 2 (P + eA)2< 3 2 



, < 3 2 (p + eA) • (VX3)|VX3| 2 (VX3) • (P - ' A 
< v^lHVxallL (Vx3)-(p + eA)<3 2 ^ 



1/2 



(C.19) 
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for any vector ip, where we have used 



<3^(P + eA)X3 <1 
which is derived from the assumption £ min > 1 1 Vo — 1 1 oo - As a result, one obtain 

|(Vxs) • (P + eA)< 3 2 | < V2hTe || | Vxsl IL ■ 
Substituting this into the right-hand side of (C.18), one gets 



^2 D i/2 



n 2 



|Axs| 



2m * y/S^-WVo- 



(C.20) 



(C.21) 



(C.22) 



Using (5.9) again, one has 



Rl'lw*{x 2 )tp 



< 



< 3 2 (p + eA)-(Vx 2 )^ 



+ 



h 2 



2m e 



<3 2 (AX2)^ 



(C.23) 



The norm of the first term in the right-hand side can be evaluated as 

(<p, (V&) • (P + eA)i?o,3(p + eA) • (V&M 
< IIIVX2IIL \\<p\\J(<P, (V&) • (p + eA)i? 0)3 (p + eA) 2 i2o, 3 (p + eA) • (V&M 



< v^lllV^I 



< 3 2 (p + eA)-(Vx2V 



where we have used the Schwarz inequality and the inequality (C.20). Therefore, 



< 3 2 (p + eA) • (VX2V| < v 7 ^ || I Vx 2 | 
Substituting this into the right-hand side of (C.23), one gets 

h 2 \\A X 



<3V*(x 2 V||<|^^|||Vx2 



+ 



2 00 



2 m e , If ■ — \\v~w 

Y mm II Ik 

The bound (C.16) follows from (C.17), (C.22) and (C.26). ■ 
From this lemma, immediately one gets 

||^*(X3)i?o,3^*(x2)i?o,2|| 2 < Const, x ||x 2 #o, 2 || 2 . 

In the same way, 

\\W(xi)Ro,2W(x2)R ,Ah < CollSt - X IIX2^0,A|| 2 . 

Substituting these two bounds into (C.15), one gets 



(C.24) 



(C.25) 



(C.26) 



(C.27) 



(C.28) 



v (Rq,a) 2w 



< Const. 



\vRo,2h 11X2-^0,2 || 2 lk^0,AX 3 ll + \\R0,AW\\ 2 ||x 2 -Ro,a|| 2 \\vRo,lXl 

(C.29) 
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Note that the ground state energy E of the Hamiltonian H 0>A satisfies E > —\\V 
Taking 

\j2m e (£ min -\\Vq~ 



h 



-k with k G (0, 1) 



(C.30) 



in the bound (D.32) in the next Appendix D, one has 



\vRonw\\ < 



(1 /t 2 )(£min 1 1 Vq I loo) 



e l3r for a region £7, (C.31) 



where r is the distance between the supports of the two functions v and w. Combining 
this inequality, (C.5), (C.6) and (C.29), one has 



v (Ro,a) 2 w 



< 



Const, x IMInJIw 



oo oo 



i - (1-k 2 )(£, 



i u i i j " || ^0 | oc i 



X 



k-Ro,2||2 n , D n H-Ro.A^lh || , D n 

-fn 11X2^0,2 1| 2 H — n — n 11X2-^0^112 



\w\ 



(C.32) 



where a = (3/3 with the above (3 of (C.30). Further, 
v (R 0:A ) 2 w 



< 



Const, x ||v||oo||w||oo 



1 - (i-«2)(£: min -||v||oo) 



x K (|supp v\ no/2 + |supp w\ no/2 ) 



)(2) 



no/2 



(C.33) 



by using the assumption (A. 5). Thus one gets the desired bound (C.l) from the bound 
(C.4) on the region 



D Decay estimates of resolvents 

In this appendix, the exponential decay bound for the resolvent (H w — z)~ l is obtained by 
using the Combes-Thomas method [36]. 

D.l The general case 

For the general case with A P 7^ 0, we use the improved version [15] of the Combes-Thomas 
method [36]. The results are given by Theorem D.2 and the bound (D.32) below. 
Consider the <i-dimensional Hamiltonian, 

H A = ^-(p + eA) 2 + V A , (D.l) 

with a general vector potential A e C 1 (A,R d ) and a general electrostatic potential V\ G 
L°°(A). Let js G C^(R d ) satisfying j$ > 0, supp js C {r| |r| < 5} with a small 8 
and f R d js(r)dxi • • • dxd = 1. Let Q be a bounded region with smooth boundary, and 
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define p(r) = dist(r, fi). Following [12], we introduce the smooth distance function p(r) = 
js * p(r). Note that, for f3 > 0, 

e~ pp H A e pp = (p + eA - ihf3V pf + V A 

2m e 

= ^A-^-(Vp) 2 -^[Vp-(p + eA) + (p + e A).Vp]. (D.2) 

We write 

e-P p H A eP p = H A + i(3J (D.3) 

with 

H A = H A -^-{Vp) 2 and J = -A- [Vp • (p + eA) + (p + eA) • Vp] . (D.4) 
We take C > satisfying 

-||VII|oo-|^|||Vp|||L + C'o>0, (D.5) 

where = max{±V A , 0}. Then one has 

#A + Co > Co > (D.6) 

with some constant Co. We define 

X E+l£ = Hk ~~ E ~~ l£ for E,e <E~R, (D.7) 

-"A + (so 

and define 

Y = (H A + C,)- l ' 2 J{H A + Co)- 1 / 2 . (D.8) 

Let us estimate the norm of this operator. From the expression (D.4) of the 
operator J, one has 

ll^ll < A [ll^ 1/2 Vp • (p + eA)/? 1 / 2 ! + \\R 1/2 (p + eA) • VpR 1/2 \\] , (D.9) 

where we have written R = (H A + Co) -1 . Take tp = R 1 / 2 (p with <p e L 2 (A). Using the 
Schwarz inequality and (D.6), one has 



(V>, (p + eA) • VpRVp ■ (p + eA)^) 



< Jty, (p + eA) 2 ^) x (p + eA) • WpR\Wp\ 2 RWp ■ (p + eA)^) 



< A/2mr|||Vp|||oo||¥»||V(^,(P + eA) • Vpi? 2 Vp • (p + eA)^>) 
< 



l2m e 



Co 



III Vp| HoolMIV^ (P + eA) • VpRVp ■ (p + eA)V). (D.10) 
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Therefore 



Similarly 



(V, (p + eA) ■ VpRVp ■ (p + eA)V) < ^|| |Vp| ||^|M| 2 . (D.ll) 



(V>, Vp • (p + eA) J R(p + eA) • Vpip) 



< II I Vp| llooll^ll (ip, Vp • (p + eA)i2(p + eA)2 J R(p + e A) • Vp^ 

< V^TeW | Vp| lUHVll y Vp • (p + eA)£(p + eA) • Vp^) . (D.12) 

This implies 

, Vp • (p + eA)i?(p + eA) • Vp^) < 2ro e ||| Vp| ||^f < ?^|| 1 Vp| || 2 JM| 2 - (D.13) 
Substituting these bounds into (D.9), one has 

imi<^^ll|Vp|||oo. (D.14) 

Consider the situation that the Hamiltonian H\ has a spectral gap (E-,E + ), and we 
take E € (E_,E + ). We define d± := dist(<r(X E ) D R^O), and w ± = P±m, where P± is 
the spectral projections for X E onto the subspaces corresponding to the sets a(X E ) n R ± , 
respectively. We take /3 satisfying E + — E > h 2 f3 2 \\\Vp\\\l /(2m e ). Then one has 

^^flWOH),^ (m5) 

P+ + Co 

and 

> E ~ E Z =; <$_. (D.16) 
P- + Co 

Lemma D.l Suppose that the parameter f3 satisfies the condition, 



Then 

\\{X E+i£ + ipY)u\\ > ^min{d+,cL}||u|| (D.18) 

for e E R. 

Proof: From the bound (D.14) for ||Y|| and the assumption on f3, one has 

1 r~— 1 



\Y\\ < -y/6 + 6- < -^d + d_. (D.19) 
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Note that X E+ie = X E —ieR with R = (H A +C ) 1 . Using this and the Schwarz inequality, 
one gets 

i\\\\(X E+i£ + ipY)u\\ > Re((u + -u_),(X E -ieR + ipY)(u + + u_)) 



u 



> d + \\u + \\ 2 + cL||w_|| 2 - 2/3Im(u + ,Yu-) 

> 7; ( d +\\ U +\\ 2 + d -\\ U -\\ 2 ) ■ 



(D.20) 



y/2m e 



This implies the desired bound. 
We write 

^lllVpH 

in terms of the parameter k G (0, 1). Substituting this into (D.15), one has 



Ea. — E x k 



6, 



(D.21) 



(D.22) 



E + + C 

Further, by substituting these into the bound (D.17), the maximum value of k satisfying 
the bound is obtained as 



k = 



C (E - E. 



\ C (E - E_) + 16(£ + + C )(£_ + C ) 



< 1. 



As a result, we can take 



y / 2m e 



C (E+ - E)(E - E. 



M\ |Vp|||oo \ C {E - E_) + 16{E + + C )(£_ + C ) 



(D.23) 



^D.24) 



Theorem D.2 Let E be in the spectral gap (E_,E + ) C R of the Hamiltonian Ha- Let 
v,w be bounded functions with a compact support. Suppose that the boundary of supp t> is 
smooth. Then 



v (H\ — E — ie) -1 w < Ci\\v\\ 00 \\w\\ 0D e~ l3r foreEH, 
where r = dist(supp v , supp w), f3 is given by (D.24), and 

1 E + + Cq E_ + Cq I 



C 1 = Const, x C 1 max 



(1 -k 2 )(E + -E)' E-E_ 



(D.25) 



(D.26) 



with k given by (D.23). The two real constants C and C satisfy the conditions H + C > 
C > and (D.5). 

Proof: Let ip be in the domain of the operator Ha- Then 



(h a + i(3J -E-ie 



[H k + C ) 1/2 (X E+te + i f3Y)(H A + C ) 1/2 </? 



> C 1/2 {X e+1£ + iI3Y){Ha + C g ) 1/2 

> ±Cl /2 nw{d + ,d-}\\(H A + Co) 1/2 <p\ 



> -C niin{d + ,d_}||^||, 



(D.27) 
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where we have used the inequality (D.18) and if a + Co > Co- Taking 

1 



ip — e 



-Pp. 



H A -E-ie 



one has 



e- f)p (H A -E-ie)- 1 e f>p u < d\\u 



(D.28) 



(D.29) 



where we have used (D.3), (D.15) and (D.16). Choosing Q = supp v in the definition 
of the distance function p(r) and using this bound (D.29), the desired bound (D.25) is 
obtained as 



Const. 



v (H\ — E — ie) 1 w 



< 



ve 



-to (Hk - E 



ie w 



< v 



< CiHulloole-^iul 

< Const, x CAv lloollwllooe - ^. 



(D.30) 



Next consider the case with H A > E. Then the Schwarz inequality yields 



(H A — E + i/3J)(p > Re((p, (H A - E + i/3 J)<p) 
= (<p, (H A - E)<p) 



> 



E -E-h 2 f3 2 \\\Vp\\\ 2 J(2m e 



(D.31) 



for tp in the domain of H A . Here E is the ground state energy of H A . Therefore, in the 
same way as in the proof of Theorem D.2, one has 



1 



Hk-E 



-w 



< 



Const, x |M 



w 



oo oo 



-e'to for0</5< 



^2m e (E - E) 



E -E-h 2 pmVp\\U(2m e ) ^ h\\\Vp\l 

(D.32) 

where v, w and r are the same as in Theorem D.2. 

Finally we consider the decay of (H A — E — iy)~ l with ?/^0. To begin with, we note 
that 



H A - {h 2 (3 I /{2m e )){VpY -E-iy + i(3J 



> 



(if, [H A - (h 2 p 2 )/(2m e ))(y P ) 2 -E-iy + t{3j] V ) 



[H A - (^ 2 /(2m e ))(Vp) 2 -E]p)\ + \(<p, (y - (3J) V )\ 2 . (D.33) 
Let E' be a real number satisfying 

(D.34) 



h 2 3 2 



2m P 



C 2 



with a positive constant C 2 . Then if the vector (p satisfies (ip,H A (p) > E'\\ip\\ 2 , the right- 
hand side of (D.33) is bounded from below by C2 1 1 1 1 2 - 
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Thus it is sufficient to consider the case that ((p,H A <p) < E'\\(p\\ 2 . Note that 

\(<P,(y-0J)<p)\ > \\v\\\<p\\ 2 - 0\(<P, J<p)\ 
The expectation value in the right-hand side is evaluated as 

\(<p,J<p)\ < — \(<p,Vp-(p + eA)<p)\ 
m e 

< — IHVpllUMlVfo (p + eA)V) 
nip 



(D.35) 




< hd — ll|Vp|| 



< ^|||Vp||| 0OV /2(^ + +||VA|| 0O )/m e || V 9|| 2 . 
For a given y, we choose a small (3 and a small C 2 to satisfy 



\y\ > ^HlVpllUv^ + +\\V A \U/m e + C 2 . 
Then these bounds yield 



(D.36) 

(D.37) 
(D.38) 



\(<p, ( y -pj) v )\>c 2 y\\ 2 . 

This implies that the right-hand side of (D.33) is bounded from below by the same C2||y?|| 2 
Thus, in both of the cases, one obtains 



H K - {h 2 (3 2 /{2m e )){Vp) 2 -E-iy + if3j] <p\\ > C 2 \\(p\ 
In the same way as the above, this leads to the decay bound, 



v(H\ — E — iy) l w < Const, x C 2 1 ||w|| 00 ||tf|| 00 e 



(D.39) 



(D.40) 



Clearly, for a small \y\, both of the parameters (3 and C 2 must be small. But the resolvent 
always decays exponentially at large distance for any y ^ 0. 

D.2 The Landau Hamiltonian with Ap = 

Following [12, 13], we obtain the exponential decay bound (D.85) for the resolvent in 
Theorem D.7 below. 

Consider the two-dimensional single electron in the uniform magnetic field and with 
a electrostatic potential V\. The Hamiltonian has the form, H\ = + V\, on the 
rectangular box A with the periodic boundary conditions as in Section 2, where H L is 
the Landau Hamiltonian (3.1). We also impose the condition of the flux quantization, 
\A\/(2tt£ 2 b ) G N. We assume V\ G C 2 (A). We denote by Q^l the spectral projection onto 
the n + 1-th Landau level whose energy eigenvalue of is given by £ n = (n + l/2)hu c 
with uj c = eB/m e . We introduce one-parameter families of operators as 



H L ((3) = e^H^r, H A ((3) = H h {(3) + V A and Q<$(/3) = e^Q^e^ (D.41) 
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for (5 e R. Here the distance function p(r) is given in the preceding subsection. We write 
Hl(P) = H L + i(3J, where the operator J is given by (D.4) with A = A , and 



h 2 B 2 



(D.42) 



Lemma D.3 Let z be a complex number satisfying dist(a(H^), z) > huj c /4, where c(H^) 
is the spectrum of the Hamiltonian H^. Then there exists k(z) G (0, 1) which depends only 
on z such that the following bound is valid: 



< 



8 

TlUJr 



for any real f3 satisfying < K[z)i^ . 



(D.43) 



Proof: Note that, for a vector ip in the domain of the Hamiltonian and for a small (3, 



[H L -z + t(3J)<p 



> 
> 



l + ipjfa-z)- 1 ] (H L -z)<p 
I-I/^IIj^l-z)- 1 !] \\(H L -z)<p 

J{H h -z)- 1 dist(<7(ff L ),z)|M 



1 



- 4 



1 - 



J(Hl-z)- 1 



(D.44) 



Therefore it is sufficient to show \{3\ J(H L — z) 1 < 1/2 for a small (3. Since 

J 

one has 
J(H L - z)' 1 



iti 2 A ^ n / A \ 

Ap Vp- p + eAo , 

2m» m P 



(D.45) 



h 2 

< l|Ap|| c 



2m P 



1 2/i 

T , , + — IHVp||Umax||(p s + eA , s )(H L - z)'^]. 

dist(o-(-HL), z) rn e s 



The norm of the operator in the second term is evaluated as 
(<p, (H L - z*)-\p s + eA 0yS ) 2 (H L - z)- l V ) 



(D.46) 



h 2 l3 2 
2m e 



< 2m e if, (H L - z 



< 2m, 



IIIVpHIl) (£l-*)V 



_dist(a(H L ),z) 



+ 



z\ + h 2 (3 2 \\\Vp\\\ 2 J(2m e 
dist(a(H L ),z) 2 



(D.47) 



Combining this, (D.46) and the assumption on z, the desired condition for f3 in (D.43) 
can be obtained. ■ 



(n) 

We introduce the integral representation of the projection Qq ^ as 



Qo 



(n) 



2m 



dz' 



z'-h,: 



(D.48) 
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where the closed path 7 encircles the spectrum of the n + 1-th Landau level. Further we 
choose the path 7 such that the length of the path is bounded as I7I < 3hu c , and that the 
distance from the spectrum of the two Hamiltonians H^, satisfies 

dist (7, <r(i7 L )) > hu c /4 and dist(7, a(H L )) > hu c /A (D.49) 

for any real (3 satisfying \(3\ \\ | Vp| ||oo < ■ Then, from the above lemma, there exists 
K n e (0, 1) which depends only on the index n such that the representation, 



2ni 



dz' 



(D.50) 



is well defined for any (5 satisfying \(3\ \\ | Vp| ||oo < K 



-1 



Lemma D.4 Assume the above condition \ j3\ \\ | Vp| ||oo < ^n^B 1 - Then the following bound 
is valid: 



— ^0,0 < 



with 



r (n) _ 48 



?b|| A^aIIoo + 2(4 + V8n + 29) max \\djV K \\< 



(D.51) 
(D.52) 



where we have written V = (d x , d y 
Proof: Note that 



1 

2tH 
1 

27ri 



-Va - Va 



1 



The commutator in the right-hand side is computed as 



z'-H L ((3) <z>-H L (P) 

1 [HMv A \- 1 



l z'-H L (py 



[Hl(P), v a ] = -^-av a - — vv A • n, 

2m e m e 



(D.53) 



(D.54) 



where we have written LI = (II X , H y ) = p + eA — ih(3V p. From these observations, one 
has 



\[Q^M,v A 



< 



< 



3h 2 u c 
2nm e 



lll^fHAVAlU + ll^llmaxll^VAlU £ 

j=x,y 



48 



7T 



ft 



+ — - maxll^lUmaxlin^)!! 

Zm e £=x,y 3=x,y 



(D.55) 



where we have written R(/3) = {zl — H L (fl)) 1 for short, and used the bound (D.43). 
Therefore it is sufficient to estimate ||rijit!(/3)||. Note that 



E iiiwvi 



3=x,y 



j=x,y 

2m e (R(P)<p,H L (J3)R(P)<p) + 2ihp £ (R(P)<p, (djp^R^) . 

(D.56) 
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Further 

max || UjR((3)ip || 2 < 2m e (\\<p\\ \\R(P)<p\\ + \z'\ \\R((3)tp\\' 

j=x,y V 

+ ^iiiVpHUll^^iimaxiin^^)^!! 

^ ¥^( 8 ™ + 13 )IMI 2 + -7SIIMI max ll^WVIh (D.57) 

where we have used (D.43), < dist(£ n ,z') + £„ and /3|||Vp|||oo < t~B~ • Solving this, one 
has 

WUjR^W < —=^( y 4 + ^8n + 29) for j = x,y. (D.58) 
VheB 

Substituting this into (D.55), the bound (D.51) is obtained. ■ 



Lemma D.5 For any given e G (0, 1), there exists K n>e G (0, 1) such that k Uj€ depends only 
on e and the index n of the Landau level, and that, for any real (5 satisfying \(5\ || | Vp| ||oo < 
K n ,e£~B~, the following bounds are valid: 



(n) 



< e 



and 



[H h ((3) - z) Q§1((3) - Qo"i < e^c for z satisfying dist(z, £„) < A£ n 
with a positive constant A£ max . 
Proof: Note that 



(D.59) 



(D.60) 



Q { oXm-Q { oX = i-I dz> 



z'-H h ((3) z'-H L 



Hi 



1 



dz' 



1 



z' - H h {(5) 



2n2 



2m, 



{Vp) 2 + i(3J 



z> - H, 



(D.61) 



for any f3 satisfying \f3\ < K n £ B . Therefore the norm of the left-hand side is evaluated as 



< 



3hu r 



sup 

2tc z '(zr 

12 



2«2 



z' - H L ((3) 
< -(2^ 2 |||Vp|||L + e ') ) 

71 v 7 



2m f 



•IHVpl 



|2 

I oo 



z 1 - Hi 



J- 



z> - Hi 



(D.62) 



where we have used the bound (D.43) and \\(z' — H L ) 1 || < A/(hcu c ), and we have chosen 
P so that, for a given small e', 

J , X „ <e' (D.63) 



z' -Hi 
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which can be proved in the same way as in Lemma D.3. The resulting bound (D.62) with 
a small j3 implies the desired bound (D.59). 

In order to obtain the second bound, we note that 



[MP) - z) [Qg(/3) - Q$\ =^-J dz' 



-1 + 



z — z 



z' - H L (P) 



(H L (f3)-H L )- 



In the same way as in the above, the norm is estimated as 



i(n) 



< 



3hu c 
2ir 



9 + 8- 



[H L (P)-H h )- 



Z' - Ey 



z'-H L 
(D.64) 



(D.65) 



Here we have used \z' — z\ < dist(£„, z') + dist(£ n ,z) < Tiu c + A£ max . The norm of the 
operator in the right-hand side is already estimated in the above. ■ 

We write z = E + ie with E, e e R. 

Lemma D.6 Suppose £ n -i + IIV^H^ + bJhio c < E < £ n+1 — ||Va~||oo — 25 + hu c for n = 
0, 1, 2, . . . with some positive constants 5± and with £_i = — oo. Let tp be a vector in the 
domain of the Hamiltonian. Then the following bound is valid: 



(H A (P) -z)[l- Qi>|| > cQhu c ||[1 - Q<> 
for any (3 satisfying \(3\ \\ | Vp| ||oo < ^ n ^li > where 

5 + , forn = 0; 

, min{<5 + ,<L}/2, for n = 1, 2, . . ., 



(D.66) 



r {n) _ 
°0,1 — \ 



(D.67) 



and 



28 + , 



mm 



for n = 0; 
-.,^\, forn = 1,2,.... 



(D.68) 



^||Ap|| 00 /|||V P ||| 00 + 2 v ^r^T ; 
Proof: We write ip — [1 — Q^\](p, and decompose the vector ip into two parts as 



l P+ = J2 Qo aV and <p_ = ^ Q^W- 

j>n j<n 



(D.69) 



Then one has 



(H A (P)-z)^\\ > Re\v + -<p-,\H A -E 

> huj c 5 + \\ip + \\ 2 + huj c o~-\\ip- 



2n2 



2m. 



(Vp) 2 -ie + ifiJ) y?+ + y?- 



2/9Im(v?+, J<p_), 



(D.70) 
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where we have used 

£n+l — \\V, 

and 



A lloo 



-E-h 2 f3 2 \\\Vp\\\ 2 J(2m e )>5 + hu c 



E-Sn-!- ||y A + |U >L?WJ C 



(D.71) 



^D.72) 



which are easily derived from the assumptions. Clearly, in the case with n — 0, one has 
ip- = 0, and so the statement holds. For the rest of the cases, we use the following bound: 



\(<P+,J<P- 



< 



< 



< 



h 2 



2m f 



ll A Plloo||^+|| \\<P-\\ + — \{<p+, Vp- (p + eA )^ 



rn. 



h 2 



2m, 



■||Ap|| 



n 



»ii¥>+ll \\<P-\\ + — V(^+> |Vp|V+)(<P-, (P + eA )V- 

?7l e 



2m, 



■||Ap||oo + ^ 



'2£ n -i . 



m P 



|Vp| 



IV+IIIIV- 



(D.73) 



where we have used (D.45) and the Schwarz inequality. Combining this, (D.70) and the 
assumption on f3, the desired bound is obtained. ■ 

Let us estimate \\(H A ((3) — z)ip\\ for a vector ip in the domain of the Hamiltonian and 
for z E C. We take (3 satisfying \f3\ \\ | Vp| ||oo < ^b 1 niin{K„ ie , K,' n } for a given e. Note that 



Q { &(H A ((3) - z)<p 



> 



e-l»{H A - z)Q^ V 



The first term in the right-hand side can be evaluated as 



(D.74) 



e-f»(H A - z)Q<ftef»<p 



(£ n + V A - z)e-^Q^ p <P 



> AE 

> AE 

> AE 



Q$14-ae\\[q$1(p)-q%1 



Q<> -eAEM, 



where 



A£ = inf \S n + V A -Rez\ 



(D.75) 



(D.76) 



and we have used the inequality (D.59). The second term in the right-hand side of (D.74) 
can be evaluated as 

Q { A(Ha(P) - *)<P - e~ Pp {H K - z)Q$e? p <P 
Q$(H A (P) - z)<p - (H A ((3) - z)Q$(J3)<p 
Q { A - Q { &(P)} (H A (Jf) - z)tp\ + 

QtX - Q^l(P)} (Ha(P) - z)<p\ + 
(H A (tf-z)<p\\+c§j>e B \\<p\\, 



< 
< 
< 



[Qixl(P),v A 



(D.77) 
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we have used the inequalities (D.51) and (D.59). Substituting these bounds into (D.74), 
one has 

(1 + e)\\(H A ((3) - z)<p\\ > AE \\Qto<p\\ - (eAE + C$£ B ) |M|. (D.78) 



Note that 

(n) 



I - Q%{)(H A ((3) - z)<p 



> 



> 



> 



[1 - Q$(P)](Ba(P) - z)<p 
(H A (f3) -z)[l- qM(M4 - ||[QSS(^), H A (0)]<P 
{H A (P) - z)[l - Q$\<p\\ 



iQ ( A(P)-Q ( A)](H A ((3)-z)v\\. (D.79) 
Using the inequalities (D.51), (D.59), (D.60) and (D.66) for this right-hand side, one has 
(1 + e)\\(H A (P) - z)<p\\ > C$hu c ||[1 - Q$]<p\\ " (etwc + e||^A|U + C$i B ) \\<p\\ (D.80) 
Combining this with the above inequality (D.78), one has 



> 



:i+e) 



1 



+ 



[H A (I3) - z)<p\ 



ri( n ) 
°o,i 



1 + 



\\VaI 



— U 0,0 l B 



+ 



AE cSitHv, 



\(p\\. (D.81) 



Assume that the energy E satisfies the condition in Lemma D.6 with the positive 
constants 5± which are independent of the strength B of the uniform magnetic field. 
Under this assumption, C^} and n' n in Lemma D.6 can be chosen to be independent of B 
except for a small B. Further assume that AE satisfies AE > AS > with some constant 
AS which is independent of B. Then there exists B$ such that, for any B > B$, 



-l 



< 1/3. 



(D.82) 



C$£ B [(AE)' 1 + (c$hu, c y 

Moreover we can choose e satisfying 

e{l + (C$y\l + \\V A \U(hu c )]}<l/3. 

Substituting these into the above bound (D.81), one has 

C^iASy'UH^P) - z)<p\\ > M for any B > B$, 

where the positive constant Cq^ depends only on the index n. Therefore a similar decay 
estimate for the resolvent is obtained in the same way as in Theorem D.2. We summarize 
the result as 



(D.83) 



(D.84) 



Widths of the Hall Conductance Plateaus 



74 



Theorem D.7 Let v,w be bounded functions with a compact support, and suppose that 
the boundary of the region supp v is smooth. Write z = E + ie with E, e G R. Assume that 
the energy E satisfies the condition in Lemma D. 6 with the positive constants 5± which are 
independent of the strength B of the uniform magnetic field. Further assume that AE of 
(D.7 6) satisfies AE > AS > with some constant AS which is independent of B. Then 
there exist B^l and k n which depend only on the index n of the Landau level such that 



Q{n) 

v(H A -E-ie)~ 1 w < -^\\v\\ 00 \\w\\ 00 exp[-R n £'] 3 1 'r} for any B > B$ ■ (D.85) 



Here r = dist(supp v, supp w) and the positive constant Cq$ depends only on the index n. 

E Proof of Lemma 5.1 

The first inequality (5.18) can be obtained as 

Waiifa + eAJItyf = (^^^Pi + eA^a^ipt + eA^R^) 



< |H|^#*(Pi + eA) 2 i^) 

< 2m e |H|L {||i?|| + [\E\ + \\(V - + VJ)\U\\R\\ 2 } 



2 

(E.l) 



for i — x,y and for any vector ip. Here we have used 



u> ) Moo 



R 



Y, R*(Pi + eA t ) 2 R < 2m e R* [H u + || (V ~ + V' 

i=x,y 

< m e {(R + R*)+2[\E\ + \\(V - + Vj)\UR*R}. (E.2) 
In order to get the second inequality, we first note that 
\\( Pi + eAi)R(p + eA) ■ aupf 



'ip, ol ■ (p + eA)R*{pi + eAi) 2 R(p + eA) • ctip) 
< m e (tp, ol ■ (p + eA) (R + R*) (p + eA) • oup) 

+ 2m e [\E\ + || (V - + K7)lloo] (V, a • (p + eA)R*R(p + eA) • cup) . (E.3) 
Using the Schwarz inequality, one has 

(tp, a ■ (p + eA)R*R(p + eA) ■ atp) 



< \\\ct\\\ ^WtpW )J (tp, a ■ (p + e A) R* R(p + e A) 2 R*R(p + e A) -atp). (EA) 

Combining this with the inequality (E.2), one obtains the bound (5.20). Similarly, 

\(tp, a ■ (p + eA)R(p + eA) • atp)\ 
< HlalllooHI^/^a- (p + eA)i?*(p + eA) 2 J R(p + eA) • a<p). (E.5) 
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Combining this with the inequality (E.2), one obtains 

\(ip, a • (p + eA)R(p + eA) • ctip)\ 2 

< 2m e |||a||| 2 J|^|| 2 {](</?, ex ■ (p + eA)R(p + eA) • atp)\ 

+ [\E\ + || (V - + K7)lloo] a ' (P + eA)i?* J R(p + eA) • «y?)} 

< 2m e |||c*||| 2 JM| 2 |(p,a- (p + eA)i?(p + eA) -a^)| 

+ 4m 2 |||a||| 4 0O / £iii (l + / B , i? )||^|| 4 , (E.6) 

where we have used the bound (5.20). Solving this quadratic inequality, one has 

\(tp, a ■ (p + eA)R(p + eA) • a<p)\ < 2m e |||a||| 2 (l + /^)|M| 2 . (E.7) 

Substituting this and (5.20) into the right-hand side of (E.3), the desired bound (5.19) is 
obtained. 

F Proofs of Lemmas 6.2 and 6.3 

For the purpose of this appendix, we prepare the following three lemmas: 

Lemma F.l (i) Let £,£' be odd integers larger than 1 such that £' is a multiple of '£. Let 
A good be the event that no two disjoint y-bad boxes of size 3£ with center in Ye n As^(z) 
exist. Assume Prob[A 3 £(- • •) is y-good] > 1 — r\ with a small i] > 0. Then Prob(A good ) > 

l-(5f/£)V- 

(ii) Assume that the event A good occurs. Let u, v G T e such that A^(u) C A^(z) and that 
A £ (v)n(A 3 ^(z)\A^(z))^0. Then 

\\ Xe (u)R 5e/ , z (E + ie)xe(v)\\ < (Se^)' 7 "" 4 \\R 5 e, z (E + ie)\\ (F.l) 

Proof: The statement (i) follows from elementary combinatorics, 
(ii) Using the geometric resolvent equation, 

Xu{n)R U ',, = R3t,uxle(u) + V5(u)^^ (F-2) 

one has 

= x^)Ru,vW^{u)R ul ^xiiy) 

= x^)Ru,uW 5 u {n) X«(u)^',-xKv), (F.3) 

uer £ n(A 3< (u)\A £ (u)) 

where we have written Ri z for Re >z (E + is) for short. We can choose Ui from the set of u 
so that ||x^(u)i?5^ jZ x^(v)|| becomes maximal. Thus one has 

\\xi{u)R u ,, zX i{v)\\ < 8e-' e WxeWRwXeMW (FA) 
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when A 3 ^(u) is 7-good. Since the norm of the operator in the right-hand side can be 
estimated in the same way, one can repeat this procedure and construct the points, 
Ui, U2, . . . , Ufc 6 r^, as long as A3^(u&_i) is 7-good and does not hit A^(v) or <9As^(z). 

The same type of estimate can be applied to v as a starting point as follows: Using 
the adjoint of the geometric resolvent equation, 

^,.Xa(v) = A(v)Rse, v + Rw,* (W&(v))* Ru,v, (F.5) 

one has 

Xe(uk)R5t'rXe(v) = Xt{uk)Rw,zxlt{v)xt{v) 

= Xe(^k)R5i',z (W^v))* Rsi^xeW 

= xtMRw^ J2 Xe(v) (Witty))* Ru tV xe(v)- 

ver £ n(A 3£ (v)\A £ (v)) 

(F.6) 

Thus 

\\xe(uk)R5t>,*Xe(v)\\ < 8e~ 7 ' WxeMRst'pXtMW (F-7) 

when A 3 £(v) is 7-good. In the same way, the procedure yields the points, Vi, v 2 , . . . , Vj, 
and one obtains the bound, 

\\xi(u)R5t>,*Xe(v)\\ < (Se-^) k+1 WxeMR^XeMW < (8e-^) fe+J \\R 5i ,J . (F.8) 

This process moves in steps of £. The assumption that A good occurs implies that there 
may be only one cluster of overlapping 7-bad boxes. The diameter of such a bad cluster 
is at most 5£. From these observations, one has that k+j > |u — v\/£ — 4 iterations 
can be performed before the process stops on the both sides. Since |u — v| > £' from the 
assumption, the desired bound (F.l) is obtained. When the process hits the boundary of 
A 5 £/(z) without hitting a 7-bad box, we have Ik > 21' or £j > £'. Therefore the bound 
(F.l) remains valid. ■ 



Lemma F.2 Let £,£' be odd integers larger than 1 such that £' is a multiple of £ and 
satisfies £' > 4£. Assume that the event A good given in the preceding lemma occurs. Then 

IxWzfez (W&,(z))*| 
< 6 (T) (8e^Y 7 '"^/4(I^IJI^z||) + /5(|SUI^z||)]||^,z||, (F.9) 
where the functions, f\ and f$, are given by 

/4(|£|,||i2||) = ^ + ^f max {\\d m ct> n , w {z)\\J{l + f E , R ) 

file ''bfib — x,y 

+ 2v^f^) 3 max{||A0 n , 3 Hz)IL}(l + /^) 1/2 ||i?ir /2 (F.IO) 
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and 

M\E\,\\R\\) 



2V2h 



= C IW R W + ^= C ^ |Vx^(z)| +2max{||9 mX ^(z)|| 00 } (1 + f E>R )^\\R\\^ 



Ah 2 



+ ^max{||9 mX ^(z)||oo} |V X ^(z)| Jl + /*,*). 
Here 



S.UJ 



2m f 



+ 



(F.11) 



(F.12) 



4>i,w{^) = xie>(z)dixli>(z), 
and the function fE,R is given by (5.21). 

Proof: Using the adjoint of the geometric resolvent equation (F.5), one has 

= xt(*)K** (wSeiz))* + X e(z)R 5e , z (w&(z))* R 3e , z (w^(z))* . (F.13) 

Therefore 



< 



Xt>{i)Ru>,.Xw{*) (W$e(z)y<p\\ + \\xe>(z)R5i',* (w^z))* R M ,, Z (w&,(z))' 



(F.14) 



for any vector ip in the domain of the operator p + eA. 

Let us estimate the first term in the right-hand side. Using the expression (5.9) of 
W(- ■ •), one has 

X£'( z )^5£',zXs£'( z ) f^'( z ))V| 



< 



Xe'(z)R 5 e',zX S 3e'(z)v3e'(z)ip + — J2 Wxt'^Rw^wi 7 )^ + eA i)v\\ > 



i=x,y 



where we have written 



% 



(F.15) 



(F.16) 



Using the bound (F.l), the first term in the right-hand side can be estimated as 
Xe>(z)R5e>,zxie>(z)v3e>(z)¥ < J2 \\Xe(u)R 5e ,, z xe(v)\\ ||w 3 £'( z )llc 



< 12 



/ e\Z'/e-4: 



(F.17) 
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The summand in the right-hand side in (F.15) can be written as 

\\xe'(z)R5e',z</>i,M'(z)(pi + eAi)ip\\ 
= ||xf'( z )^',zW^(0 ii 3^(z)) J R 5 ^ iZ (pi + eAi)(p\\ 



< 



h 2 



2m. 



||X£'(z) J R 5 ^ iZ (A0 ii 3£/(z)) J R 5 £/ iZ (p i + eAi)<p\\ 



h 



+ — llx^"( z )^5£',z(9 i 0i i3 <;'(z))(Pi + eAj)R 5 e> jZ (pi + eAi)<p\\ 



m 



< 6(f (Se-')""- 1 !!^, 



h 2 

— maxiHA^^^z)!!^} \\R 5e , z (pi + eA^W 

m e x )f 



Ah 



+ — max {||9 m ni3 ^( z )lloo} life + e-^ifo'.-CPi + eA)<^|| 



, (F.18) 



where we have used (5.9) and (F.l) again. From (5.19) and (5.20) with ol = (1, 0) or (0, 1), 
one has 

\\R( Pi + eAVII < v^(l + f E , R ) 1/2 \\R\\ 1/2 y\\ (F.19) 

and 

IKp^ + e^CPi + e^MI <2m e (l + / B)fl )|M|. (F.20) 
Combining these, (F.15), (F.17) and (F.18), we have 



(F.21) 



Next let us estimate the second term in the right-hand side of (F.14). Note that 



(# 3 *,(z))*it^, z (w&(z))* 



h 



< Cs,u WRst'A IMI + —Cs,u, PWP + eA) • (V X3 V(z)V 



m, 



n 

+ ~ Cs *> ™^,{ll 9 ^'( z )lloo} £ ||(Pi + eA)^', 2 



i=x,y 



h 



+ — max{||,9 m x^( z )lloo} E (Pi + eA)i^',-(p + eA).(Vx^(z)^ , (F.22) 



t=x,y 



ml m=x,y 

where we have used (5.9). Thus we have 

X e(z)R 5e , z (w&(z))* R 3e ,, z (Wi e ,(z 



< 6 



/_ / 5 (|£|,||i^,z||)||i^ 



(F.23) 



in the same way. Substituting this and (F.21) into (F.14), the desired bound (F.9) is 
obtained. ■ 



Similarly, one has the following lemma: 
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Lemma F.3 Let £,£' be odd integers larger than 1 such that £' is a multiple of £ and 
satisfies £' > A£. Assume that the event A good given in Lemma F.l occurs. Then 

W 5 u ,{z)Ru>,.Xe{z) 

< 6 /T) (8 e -^)' 7 ^ 4 [/ 6 (|EU|i? 5 ,, z ||) + / 5 (|EU|i?3^z||)]||i?5£',z||, (F.24) 



where the function f 6 is given by 



8ti 2 



fe(\E\A\R\\) = 2C*, W + — max{|||V0 m , 3 Hz)|IL}(l + /^) 



+ 2V2 



>m P / n=x 'V 



rrmx{||A0 n , 3 Hz)L}(l + f E , R ) 1/2 \\R\\ 1/2 . (F.25) 



Proof: Using the geometric resolvent equation, one has 

^£'( z )X3£'( z )^5£',zX£'( z ) 

= W^{z)R 3e>zXf j{z) + W s 3e {z)R u , tZ W^{z)R be>z xe{z). 



(F.26) 



Therefore 



W£ e ,(z)R 3e ,, zX e>(z) 



< 
+ 



w£ e ,(z)x s 3e ( z )R5e>,zXe>(z) 
W£ i ,(z)R 3e , tZ W^(z)R 5e , :zX Az) 



;f.27) 



Since the second term in the right-hand side can be estimated in the same way as in the 
proof of the preceding lemma, it is enough to estimate the first term. Using (5.9), one has 



^'( Z )X3>( Z )^',zXW Z ) 



< 
+ 



^'( Z )X^'( Z )- R 5^',zX^( Z ) 

h_ 

m, 



\\(Pi + eA 'i) ( d iXle'( z )) xie>( z ) R 5e',zXe>{z) 



i=x,y 



(F.28) 



The first term in the right-hand side can be estimated by using the bound (F.l). The 
operators in the sum in the right-hand side can be written as 

(pi + eAi)<t) ij3l >(z)R u ^ z xt'{z) = (pi + eAi)R u , tZ W(<\)i M \z))R u ^ z xt<{z) 

h 2 



2m, 
ih 



(Pi + eAi)R 5£lyZ (A0 ii3£ /(z)) Rw,2.Xe{ z ) 

(pi + eAi)R 5 e :Z (p + eA) • (V(f)i^(z))R 5e ^ z xe'(z), 

(F.29) 



where we have used (5.9) again. Similarly the norm of the operators in this right-hand 
side can be estimated by using (5.19), (E.l) and (F.l). ■ 
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Proof of Lemma 6.2: Assume that the event A good given in Lemma F.l occurs. Then, 
from the preceding three lemmas, one has 



Xi'Ru,. (W&,(z))*|| < Const, x (£') 3 exp K'{ 7 (1 - U/t) - 3\og2/£}} 



x \E\(\\R. 



U,x\ 



\R, 



5£,z | 



\R. 



5£,z | 



;f.3o) 



and 



Wit, (z)Rsi, z Xi> I < Const, x (£'f exp [-^{ 7 (1 - 4^/0 - 3 log 2/£}} 

X 1^1(11^11 + ||i?5£,z||)||^,z|| (F.31) 

with the probability larger than 1 — (5£'/£) 4 rj 2 for a large \E\ and for large ||i?3£, z ||, H-Rs^zll- 
In the Wegner estimate (A. 38), we choose 



(5E)- 1 = CwKMUAsel xA(£'f. 

Then, for q = 3,5, one has 

\\R q eJ < (SE)' 1 
with the probability larger than 1 — (£')~^/A. Clearly, one has 

||^,z|| < Const, x Kz(£'f +2 



(F.32) 
(F.33) 

(F.34) 



for a large £'. Since the probability that each event occurs is larger than 1 — (£')~^/A, the 
probability that the two events simultaneously occur is larger than 1 — (£')~^/2. 

From these observations, the right-hand side of (F.31) can be bounded from above by 

exp [-£'{i(l - U/£') - 31og2/£ - log(Const. x K$\E\)/£! - (2f + 7) log £'/£'}] (F.35) 

with the probability larger than 1 — (5£' /£Yrf — (£')~^/2 for a large \E\ and for a large £'. 
Since one has 



31og2 log(Const. x K 2 \E\) \og(c K 2 \E 



+ 



< 



£ £' ~ £ 

with a positive constant Co, the proof of the lemma is completed. 



(F.36) 



Proof of Lemma 6.3: Take £' = 4+i and £ = £ k in Lemma 6.2, and assume that A 3 ^ fe (- • •) 
is a 7fc-good box with the probability larger than 1—n with rj = (4)~^- From the definition 
(6.8) of £k+i, we have 



£ 



fe+i 

4 



{,1/2 



odd 



£l /2 + 5£ k with < 5£ k < 2. 



(F.37) 



Using this identity, rj' of (6.6) can be written as 



r] = r} k+1 = 5 



-k+i 



) (4)- 2 ^ + ^(4 + i)^ 



= (4 + i)- ? 
= (4 + i)- ? 



4+C 



(4)" f + 



2 

4+5 



(F.38) 
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Therefore, if the initial scale 4 satisfies 



5 4 (4 



(l + 24 1/2 ) 4+? <V2, 



(F.39) 



then we have r/ k+ i < (4+i)~ 5 - Actually this inequality holds for a large 4 because of the 
assumption, £ > 4. 

Next we define 7^ inductively according to (6.7) as 



with 



One can easily find 



7fc+i 



7fc+i = 7fc (1 - 44/4+i) - 4 
log(c K!\E\) , (2s + 7) log 4+i 



4 



+ 



44 



-k+i, 



-d k 



- 1- 



44 



44-i N 

4 , 

4-i 



Kfc+i 



44' 



(F.40) 
(F.41) 



-fc+i, 



44 



'fc+i 



1 r— - I 4-2 



4 



44 



fc+i 



> 



44 

4- 



Note that 
and one has 



i+i/ j=o 

4 > (4-i) 3/2 > (4- 2 ) (3/2)2 



> (4) {3/2)S 



1 

< — — < exp 



4+i £ 



1/2 



^(3/2) fe log4 



(F.42) 



(F.43) 



;F.44) 



Using this inequality, the product in the right-hand side of (F.42) can be evaluated as 



n 1 



44 



3=0 



> I] {l " 4exp [-- (3/2)'' log 4] } • (F.45) 



This right-hand side is uniformly bounded from below by some positive constant if 4 > 16. 
The sum of 4 in the right-hand side of (F.42) can also be evaluated as 

k k 

J2dj < log(co^3 2 |E|)^exp[-(3/2yiog4 
3=0 3=0 

k 

+ (2s + 7) log 4 E (3/2) J+1 exp [- (3/2) J+1 log 4 
3=0 



(F.46) 



This right-hand side becomes small for B large enough because of K 3 = O(B) or 0(1), 
\E\ = 0(B) and 4 = 0(B 1 / 2 ). As a result, 7^ is uniformly bounded from below by some 
positive constant 700 . ■ 



Widths of the Hall Conductance Plateaus 



82 



G Proof of Lemma 9.2 

The difference between (9.18) and (9.19) is estimated by 
1 / 

AI := — I ^ E + E E E Ku,v^v,w^w,u| |5u,v,w,u| J • 

V t \aGA £ ue(Ze 2 )*\A* v,w a eZ?\A e ugA* v,w 



(G.l) 



In the same way as in the proof of Theorem 8.5, it is sufficient to show E[AJ] — > as 
To begin with, we note that 

| Trx £ (u)P F x £ (v)P F X £ (w)P FX£ (u)| 

< V /Tr X £ (u)PFX £ (v)PFX £ (v)PFX £ (u)VTrx £ (u)PFXe(w)P F x £ (w)P FX£ (u) 



< || X£ (u)PfX £ (v)|| ||x £ (u)PfX £ (w)|| v /Trx £ (v)P FX£ (v) v /Trx £ (w)P FX£ (w) 

< Const.||x £ (u)P F x £ (v)||||x £ (u)P F x £ (w)||, (G.2) 

where we have used (8.27). Further, Schwarz's inequality yields 



E[||x £ (u)P FX£ (v)|| ||x £ (u)P fX£ (w) 



< E 



|x £ (u)P F x £ (v) 



1/2 



E 



|x £ (u)P F x £ (w) 



1/2 



< E [\\x £ (u)P F x £ M\\] 1/2 E [|| X£ (u)P FXe (w 



,1/2 



(G.3) 



where we have used ||x £ (u)P F x £ (v)|| < 1 for any u, v. From these observations and the 
decay bound (7.17) for the Fermi sea projection, it is sufficient to estimate 



1 



T>(E E E+ E E E |tu,v*v,wfw,u| e -Hu-v|/2 e - M |u-w|/2 



(G.4) 



iGA £ ue(Z2)*\A* v,w aeZ 2 £ \A e uGA* v,w 

Consider first the case with |u — a| < E\i & with 5 G (0, 1). Using the bound, 

\t t t I < 2' 
which is derived from (9.5), we have 

5^ |^u,v^v,w^w,u| e ^ ^ e ^ ^ < Const. 



u - 


- V 


u 


— w 




u - 


- a 


2 



(G.5) 



u — a 



(G.6) 



Therefore the corresponding error is estimated by 



1 



I 



E + E 

a6A^,u6(Z|)*\A|: aeZf\A £ ,ueA*: 
\ |u-a|<si^ |u-a|<ei^ 



u — a 



Consti • ^(Const. + Const, log I) 
< j 2 • (G.7) 
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This vanishes as £ — > oo. 

When |u — a| > e\£ s , we further decompose it into two cases: (i) both v and w fall into 



u < u - 


- a and 


w 


— u 


< 


u — a 


- u > u 


— a or 


w 


— u 


> 


u — a 



The latter contribution is exponentially small in £ . Actually, one has 

1+ f f I p -M|u-v|/2 -m|u-w|/2 



E 

v,w satisfy (ii) 



< 2' 



E + E 

v:|v— u|>|u— a|, w:|w— u|>|u— a|, 



e -M|u-v|/2 e - M |u-w|/2 



< Const.e-^'l u - a l (G.t 



with a positive constant //. 

Finally consider the former case (i). To begin with, we note that 

^ u ,v^v,w^w,u = 2i {sin Z(u, a, v) + sin Z(v, a, w) + sin Z(w, a, u)} 



(G.9) 



We write a = Z(u, a, v), (3 = Z(v,a, w) and 7 = Z(w,a, u) for short. In this case, one 
notices that a, (5 G (— n/2, ir/2) and a + (3 + 7 = 0. From these, one has 



sinct + snip + sin 7| < 2 I |sina|sin — + ] snip] sin — 

2 Q 1 I _• /Ol -2 



< 2(| sin a I sin /3 + | sin /3| sin a) 
2 

(|u-v||w-u| 2 + |v-u| 2 |w- u l) , (G.10) 



< 



u — a 



where we have used 



v — u w — u 

sinal < r and |sin/3| < - 



|u — a| |u — a| 

for getting the third inequality. From these observations, we obtain 

E |fu,v*v,wfw,u| e -M|u-v|/2 e -^|u-w|/2 < Const 1 
v,w 

The corresponding contribution is estimated by 



u — a 



(G.ll) 



(G.12) 



1 



\ 



E + E 



\ |u-a|>ei^ 

This vanishes as £ — > oo. 



aeA^ue(Z2)*\AJ: aeZ|\A £ ,ueA £ 



*. u — a • 



|u-a|>ei^ 4 



< 



Const. 



(G.13) 
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H Proof of Lemma 9.7 



In order to prove Lemma 9.7, we introduce a partition {x^(u)} u of unity and prepare 
Lemma H.l below. Here Xb( u ) are C 2 positive functions with a compact support such 
that J2uXb( u ) = 1- Let Xf>( u ) denote the characteristic function of the support of X&( u ). 

Lemma H.l The following bound is valid: 

(Ps + eA s )x 5 b (u)R(z) X e(v)\\ < [Const. + {2m e \z\f /2 } ||x 6 (u)i?(z)x £ (v)|| 

+ V2hT e X6Hi?(^)x £ (v)|| 1/2 , (H.l) 



where the positive constant depends only on the strengths of the potentials, Vq, K, and on 
the cutoff functions Xb( u )- 

Proof: Let (p be a vector on R 2 . Then one has 

(^,x £ (v)i?*^)x^(u)(p s + eA s ) 2 x ^(u)i?(z) X£ (v)^) 

< 2m e (||y ||oo + IlKHoo) ||^(u)^)x £ (v)^| 2 

+ 2m e (tp, X e(v)R%z)xt(u)H„xi(^R(z)Xe(v)<p) (H.2) 



by using the inequality (p s + eA s ) 2 / '(2m e ) < H u + ||Vo||oo + 
term in the right-hand side is evaluated as 



j. Further, the second 



2m e (tp, X e(v)R%z)xi(u)H„xi(u)R(z)Xe(v)<p) 
< 2m e \z\ \\x 5 b (u)R{z) X e(v)<p( 

+ 2ft ]T ||[^(u)]i?(^)x £ (v)^||||x^(u)(p s + eA s )i?(z)x £ (v)^ 

s=x,y 

+ h 2 \\xl{u)R{z)xe{v)y [A x s b (u)}R(z)x £ (v)<p 



+ 2m e 



xt{u)R{z)xs{v)y xlHXeWv 



by using 



ih 



ti 2 



^xt(u) = --Vxi(u).( P + eA) ^ 



AxJ(u)+x?(u)fl- w . 



(H.3) 
(H.4) 



Combing these two bounds, we obtain 



(p s + eA s )x 5 b (u)R(z) X e(v) 



< 



2m, 



oo 



+ 



+ \z\) + 2h 2 £ ||9.xf(u)||2o + ft 2 ||Ax2(u) 



s=x,j/ 



x||x 6 (u)^) X£ (v)|| 2 
+ 2ft 2 H^xJ(u)||oo||X6(u)i2(z)x e (v)|| (p. + eA s ) X *(u)it:(z)x £ (v) 



+ 2m e \\x 5 b (u)R(z)xe(v)(p x£(u)x.(v) 



(H.5) 



Widths of the Hall Conductance Plateaus 



85 



Solving this quadratic inequality and using the inequality y/a + b < yfa + Vb for a, b > 0, 
the desired bound is obtained. ■ 

Proof of Lemma 9. 7: Note that 

E[|X(P F ;Q,4>)-X(P F , A ;fi)|] = E [|X(P F ; QJ P ) - X(P F)A ; n,£ P )\ I(M A )] 

+ E[\l(P F ;n,l P )-l(P FX ,n,l P )\l(M c A )], (H.6) 

where M A is the event which was introduced in the proof of Lemma 8.1, and 1(A) is the 
indicator function of an event A. The second term in the right-hand side is vanishing in 
the limit L j oo as 



E[\l(P F ;nj P ) -1(P F:A ;Q)\I(M C A )} 



< E 



\1{P F ;Q,£ P ) - J(P F , A ;fi)| 2 l 1/2 E[I(Ml)] 1/2 < Const.L' 2[K( « +2) - 3l/6 , (H.7) 



where we have used Schwarz's inequality for getting the first inequality, and we have used 12 
Lemma 8.4 and Prob(M^) < Const. L _2 ' K ^ +2 )~ 3 ^ 3 for the second inequality. 

In order to estimate the first term in the right-hand side of (H.6), we first note that 

Tr X nP F [[PF,x], [P F ,y]]xn = Tr X n(P F xP F yP F - P F yP F xP F ) X n- (H.8) 

We want to rewrite this right-hand side. We have 

TrxnP F xP F yP FX n = Tr Xn P F x( X { + 1 - X s A )P F y( X { + 1 - X {)P FXn 

= Tr X nP F x X s A P F y X s A P FX n + corrections. (H.9) 

The contributions from the corrections decay exponentially in the distance between fl and 
the support of 1 — Xa by the bound (7.17) for the Fermi sea projection P F . The rest is 
written 

Tr Xn P F x X 5 A P F y X s A P FXn = ^-j dzTi Xn R(z)x X 5 A P F y X 5 A P FXn . (H.10) 
Note that one has 

XnR(z) = X nxiR(z)=XnRA(z) XA + X nRA(z)W( X s A )R(z) (H.ll) 

from the geometric resolvent equation, Xa R(z) = R A (z) X s A + R A (z)W ( Xa ) R(z) , where 
^(xi.) = H^(x A ) + (Vu,a — Kj)x A ) an( i Kj,a is the slightly modified potential near the 
boundary of A. (The precise definition of V^ jA is given in Section 2.) The contribution 
from the second term in the right-hand side of (H.ll) is 

^- J dzTi Xn R A (z)W( XA )R(z)xx 5 A P F yxiP F Xn- (H.12) 
The integrand is estimated by 

J2 E |Trx6(u)P A (z)W(xi)P(^)xxix 6 (v)P FX6 (w)y X iftX6(u)| . (H.13) 

u:s b (u)n!l^v,w 



12 The bound (8.31) of Lemma 8.4 holds also for A = R 2 . 
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For any bounded operators A, B, 



|Tr A Xb (y)P FXb (w)B\ < ^Tr A Xb (v)P FXb (v) A* ■ ^Tr B* Xb (w)P FXb (w)B 

< Const. PIUI^II, (H.14) 



where we have used the bound (8.27). Using this inequality, one has 

Tr Xb (n)R A (z)W(x 5 A )R(z)xx S A XbMPFXb(^)yX 5 APFXb(n) 

< Const. \\ Xb (u)R A (z)W( X 5 A )R(z)x X s AXb (v)\\ \\ Xb (w)y X 5 A P FXb (u) 

< Const. ^2 L 2 \\ Xb (u)R A (z) Xb (u')\\ \\w( XA ) X s b (u')R(z) Xb (v)\\ || X6 (w)P fX6 (u)|| , 

(H.15) 

where the sum is over u' such that supp X f(u') fl supp | V% A | ^ 0. Because of the existence 
of the indicator function I(M A ) in (H.6), the first factor in the summand is estimated as 

\\ Xb (u)R A (z) Xb (u')\\ I(M A ) < Const.L K «- 2 ) +4 exp[- /Uoo L 2K / 3 ] (H.16) 

from Lemma 8.1. The second factor can be estimated by using Lemma H.l. Therefore it 
is sufficient to estimate the following quantity near the Fermi energy: 



E 



dy \\ Xb (v.')R(E F + iy)xb(v)\\ ||xf>(w)P F x 6 (u) 



Note that 

ry- 



dy \\ Xb (u')R(E F + iy) X b(v)\\ < I' dy \\x b (u')R(E F + iy) X b(v)\\ s/2 \y\ s ' 2 ~\ 



Substituting this into the above, we obtain 



(H.17) 



(H.18) 



E 



y+ 



dy \\ Xb (u')R(E F + iy)xb(v)\\ IIXb(w)P F x&(u) 



< liminf/ ^|y| s / 2 - 1 E[||x fe (u')P(P F + ^)x6(v)ir /2 ||x 6 (w)P FX6 (u) 



< liminf / dy\y\ s l 2 - 1 ¥ 1 \\\ Xb {^)R{EY + iy)Xb^) 



1 1/2 



E 



IXfe(w)P F x&(u) 



< Const. exp[— fi\u — v|/2] exp[— /i|w — u|/2], 



nl/2 



(H.19) 



where we have written I n = [?/_, e n , e n ), and we have used Fatou's lemma, Fubini- 

Tonelli theorem, Schwarz's inequality, the bounds (7.1) and (7.17). Thus, the correspond- 
ing contribution is vanishing as L j oo. 

Consequently, it is enough to consider Tr X nP FtA x( XA ) 2 P F y XA P FX n which comes from 
the first term in the right-hand side of (H.ll). Using the adjoint of the geometric resolvent 
equation, R(z) X & A = X s A R A (z) — R(z)W( X 5 A )R A (z), we have 



Tr X nP F ,Ax( X s A ) 2 P F y XA P FXn = Tr X nPF,Ax( XA j 2 P F y( XA j 2 P FAXn + correction (H.20) 



8 \2 



,S\2. 
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in the same way. The correction is vanishing as L j oo. Using the geometric resolvent 
equation again, the first term in the right-hand side is written 

Tr X nP F ,Ax(x S A ) 2 P F y(x S A) 2 PF,AXn 
= Tr XnPF,AXx{P F ,Ay(XA) 3 PF,AXn 

+ ^~l dzTT XnPF,A^XARA(z)W(xi)R(z)y( XA ) 2 P F ,AXn. (H.21) 
The integrand in the second term in the right-hand side is written 

Tr Xn P FA xx S A R A (z)W(x A )R(z)y( X S A ) 2 PF,AXn 
= Tr X ^P F , A x X { XA/ R A (z)W( XA )R(z)y( XA ) 2 P FAXn 

+ Tr X nP F , A x X {(l - XA ,)R A (z)W( XA )R(z)y( XA ) 2 P F , AX n, (H.22) 

where xa 1 is the characteristic function of the region A' which satisfies the conditions of 
(8.4). This right-hand side can be shown to be vanishing as L j oo in the same way. 
Consequently, we obtain 

Tr x^P F ,A x X A P F ,Ay(x & AfPF,AXn = Tr X uP F , A xP FjA yP F:AX n + correction. (H.23) 

The first term in the right-hand side is nothing but the desired form. ■ 



I The index formula for the switch functions 

The aim of this appendix is to give a proof of the following theorem: 

Theorem 1.1 For a fixed period £p of the potentials A LP and V^ LP in the Hamiltonian 
Hlf of (9.13) on the whole plane R 2 , the following relation is valid almost surely: 

Index(P F f/ a P F ) = 2«Tr P F [[P F , A 1>a ], [P F , A 2 , a ]], (1.1) 

where Xj jSL , j = 1, 2, are two switch functions given by 

' v ' 1 0, for x - a 1 < ' v 1 1 0, fory-a 2 <0 v 1 

with the locations a = (a±, a 2 ) G R 2 o/ the steps. 

Remark: 1. The right-hand side of (1.1) is equal to the form of another Hall conductance 
which was discussed in [10]. Elgart and Schlein [7] justified this Hall conductance formula 
within the linear response approximation under the assumption that the Fermi energy lies 
in a spectral gap. They also proved that the value of (1.1) takes the desired integer under 
the same gap assumption. As mentioned above, Germinet, Klein and Schenker proved 
the constancy of (1.1) in the localization regime, for a random Landau Hamiltonian with 
translation ergodicity, by using a consequence of the multiscale analysis . 
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2. From Theorem 9.4, we obtain that the Hall conductance using the position operator is 
equal to that using the switch functions. 

We write the index as X S (P F ; £p) = 27riTr P F [[P F , Ai ;a ], [P F , A 2>a ]]. First, we shall show 
that the index Z S (P F ; £ P ) is well defined for almost every uj. Note that 

Tr |[P F ,A lia ][P F ,A 2 , a ]| 

< ]T Tr lx £ (u)[^F,A 1)a ] X£ (v)[P F ,A 2)a ] X£ (w)| 
u,v,we(z2)* 

< £ l A i-aH - A lia (v)||A 2 , a (v) - A 2 , a (w)|Tr \ Xe (u)P FXe (v)P FXe (w)\, 
u,v,we(z2)* 

(1.3) 

where Xs( u ) is the characteristic function of the E\ x e 2 rectangular box s e (u) centered at 
u, and we have chosen the set (Z £ )* of the centers u of the boxes s e (u) so that a becomes 
a vertex of a rectangular box, i.e., a e Z 2 . Using Schwarz's inequality, we have 

E [Tr \ X s(u)P F Xe(v)P F XeM\] 



< yjE [Tr X£ (u)P FXe (v)P FX£ (u)] y/E [Tr X £ (w)P fX£ (v)P fX£ (w)] (1.4) 

and 

Trx £ (u)P FX£ (v)P FX£ (u) < ^Tr x £ (u)P FXe (u) • Tr x £ (u)P F x £ (v)P F x £ (v)P F x £ (u) 

< Const. ||x £ (u)P FX£ (v) ||, (1.5) 

where we have used the bound (8.27). From these bounds, we obtain 



E [Tr | X£ (u)P fX£ (v)P fX£ (w)|] < Const.yE [||x e (u)P F Xe(v)||]yE [||x £ (w)P F x £ (v)||] 

< Const.e^ |u - v|/2 e~ Hw - v|/2 , (1.6) 
where we have used the decay bound (7.17) for the Fermi sea projection. Note that 

|A, a (u)-A ja (v)| = (°' ! 0r ^t r t! ; (1-7) 
1 h v 1 ] ' y 71 | 1, for (uj - aj)(vj - dj) < 0, v ' 

and 



yj\xi - yi\ 2 + \x 2 - V2? > \xi - yi|/2 + \x 2 - y 2 \/2. (1.8) 
Combining these, (1.3) and (1.6), we obtain 

E [Tr | [P F , Ai >a ] [P F , A 2)a ] |] < Const. ^ e -^K-»il/4 e -^i-ai|/4 e -Ml«2-^|/4 

u,v,w 

Xe -/i|w 2 -02|/4 e -^|t;2-02|/4 e -;i|tui— ui|/4 < ^ ^ 

Thus the operator [P F , Ai i& ][P F , A 2>a ] is trace class for almost every uj. 

Next we show that the index X s (P F ;£p) is independent of the locations ai,a 2 of the 
steps of the switch functions A Jia . Let a', a e R 2 . Then we have 

Tr P F [[P F , A 1;a ,], [P F , A 2 , a /]] - Tr P F [[P F , A 1>a ], [P F , A 2 , a ]] 
= Tr P F [[P F , (A 1>a , - A 1>a )], [P F , A 2 , a ,]] + Tr P F [[P F , A 1>a ], [P F , (A 2 , a , - A 2 , a )]]. 

(1.10) 
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We will prove that the first term in the right-hand side is vanishing because the second 
term can be handled in the same way. We choose e = (£1,62) so that both of a' and a 
satisfy a'.aeZ^b) :=Z^-b with some b e R 2 . We denote by (Z £ (b))* the dual lattice 
ofZ 2 (b). 

Lemma 1.2 For u, v e (Z 2 (b))* ; the following bound is valid: 

E [Tr | X£ (u)P FX£ (v)|] < Conste-"''^' (1.11) 
with some positive constant //. 
Proof: Note that 

E[Tr | X£ (u)P fX£ (v)|] < Yl E[Tr | X£ (u)P fX£ (w) X£ (w)P fX£ (v)|] 

we(z2(b))* 



we(z2(b))* 



Xv /E [Tr X£ (v)P FX£ (w)P FX£ (v)]. (1.12) 

Further, we have 

Tr X£ (u)P F x £ (w)P FX£ (u) < ^Tr Xe (u)P F x e (u)^Tr x £ (u)P f x £ (w)P fX£ (w)P fX£ (u) 

< Const. ||x £ (u)P FX£ (w) II, (1.13) 

where we have used the bound (8.27). Combining this, the decay bound (7.17) for the 
Fermi sea projection, (1.12), we obtain 

E[Tr |x £ (u)P F x £ (v)|] < Const. £ e -d«-w|/2 e -A«|w-v|/2 < C onst.e^' |u - v| . (1.14) 

we(z2(b))* 

■ 

Now let us consider the first term in the right-hand side of (1. 10). We write AA for 
Ai, a ' — Ai ia for short. 

Lemma 1.3 We have 

E [Tr |AA[P F ,A 2 ,a']|] < 00 and E [Tr |AAP F [P F , A 2 , a ']|] < 00. (1.15) 

Proof: Without loss of generality, we can assume a[ > a±. Then we obtain 

E [Tr |AA[P F , A 2 , a ']|] < £ E [Tr | AAx £ (u)[P F , A 2 , a ']x £ (v)|] 

u ,ve(z2(b))* 

< E l A l,a'(u) -Ai,a(u)||A 2) a'(v) ~ A 2 , a '(u)l 

u ,ve(z2(b))* 

xE[Tr|x £ (u)P FX£ (v)|] 

< ^ E e _M ' |ui ^ l|/2 e _M ' |u2 ~ a 2 |/2 e^' |l,2_a2l/2 < 00, 

ai<ui<a[,U2 u i' u 2 

(1.16) 
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where we have used (1.7), (1.8) and Lemma 1.2. 
Similarly, we have 

E[Tr |AAP F [P F ,A 2 , a ,]|] < ]T E [Tr |AA X£ (u)P FXe (v)[P F , A 2 , a ,] X£ (w)|] 

u,v,w 

< l A l,a'( u ) - Ai,a(u)||A 2 ,a'(w) - A 2 , a /(v)| 
u,v,w 

xE[Tr| Xe (u)P F x £ (v)P FX£ (w)|] 

< Const. l A i,a'( u ) - Ai, a (u)||A 2ia /(w) - A 2 ,a'( v )l 

u,v,w 

xe ~ M |u-v|/2 e - M |w-v|/2 

< Const. ]T e -H«i-«l|/4 e -*i|«l-ii;i|/4 

ai <ui <a' 1 ,ui ,wi 

X ^ g-/x|u 2 — U2|/4 e -/i|t)2-a^|/4 e -^|tu 2 -oi|/4 < ^ (1.17) 

where we have used the bound (1.6). ■ 
Relying on this Lemma 1.3, we have 

TrP F [[P F ,AA],[P F ,A 2 , a ,]] = TrP F AA(l-P F )[P F ,A 2 , a ,]+Tr[P F ,A 2 , a ,](l-P F )AAP F . (1.18) 
Further, the first term in the right-hand side is written 

TrP F AA(l-P F )[P F ,A 2 , a ,] =TrAA(l-P F )[P F ,A 2 , a ,]P F = TrAA(l-P F )[P F ,A 2 , a ,], (1.19) 
where we have used (1 — P F )[P F , A 2 , a '](l — Pf) = 0. The second term becomes 

Tr[P F ,A 2 , a ,](l-P F )AAP F = Tr[P F ,A 2 , a ,](l-P F )xsu PP AAAAP F 

= TrAAP F [P F ,A 2 , a ,](l-P F )xsu PP AA 

= TrAAP F [P F ,A 2 , a ,](l-P F ) 

= Tr AAP F [P F , A 2>a /], (1.20) 

where Xsu PP aa is the characteristic function of the support of AA, and we have used 
P F [P F , A 2ja /]P F = 0. As a result, we obtain 

Tr P F [[P F , AA], [P F , A 2ja ']] = Tr AA[P F , A 2ja /]. (1.21) 

This right-hand side is decomposed into two parts as 

Tr AA[P F , A 2 , a ,] = Tr A\ Xe [P F , A 2 , a ,] + Tr AA(1 - X e) [Pf, A 2 , a '] (1.22) 

with the characteristic function xe of the square box centered at r = with a sufficiently 
large sidelength t. Since we have 

E [Tr | X ,P F |] < £ E [Tr lx^FXe(v)|] < oo (1.23) 
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from Lemma 1.2, the first term in the right-hand side is vanishing by cyclicity of the trace. 
The second term can be evaluated in the same way as in the proof of Lemma 1.3. In 
consequence, it vanishes as £ j oo. Thus we obtain Tr P F [[P F , AA], [P F , A2, a ']] = 0. Since 
the second term in the right-hand side of (1. 10) can be handled in the same way, the index 
Z S (P F ; £p) is independent of the locations a±, a 2 of the steps of the switch functions. 
Using this property, the index is written 

27T? 

UPf; M = yE E Pi2(v, w, u; a) - D 21 (v, w, u; a)] S(u, v, w, u) (1.24) 

Vl a£A e u,v,we(Z 2 )* 

with 

£> 12 (v,w,u;a) : = [Ai,a(v) - A lia (w)][A 2 , a (w) - A 2 , a (u)] (1.25) 

and 

L> 21 (v,w,u;a) := [A 2 , a (v) - A 2>a (w)][A 1)a (w) - A 1>a (u)], (1.26) 

where both and Ve are the same as in (9.18), and 5(11, v,w,u) is given by (9.20). We 
also write 

l £ s (Pr, n,£ P ) = — £ E E l D Mv, w, u; a) - D 21 (v, w, u; a)] S(u, v, w, u), 

V (- ueA* v,we(Z 2 )* aGZf 

(1.27) 

where A| is given by (9.21). 

Lemma 1.4 The following holds: E [\l s {P F ;£p) - T £ s (P F ; ^, £p)\] -> as |Q| ] 00. 

Proof: To begin with, we note that 

|E[Tr X£ (u)P F x £ (v)P F x £ (w)P F x £ (u)]| < E [Tr |x £ (v)P fX£ (w)P fX£ (u)|] 

< Const.e^ |v - w|/2 e^ |w - u|/2 (1.28) 

which is derived from (1.6). Using this, (1.7) and (1.8), we have 

E[\l s (Pr,£p)-l £ s (P F ;n,£p)\} 

< Const - XI E + E E E e~^ Vl ~ ai ^ 4 e~^ Wl ~ ai ^ 4 e~^ Wl ~ Ul ^ 4 

^ aGZ2\A^ueA* aGA^ ue(Z 2 )*\A* 

x ^2 e -f*\v2-w 2 \/<i e -fi\w2-a2\/4: e -n\u2-a2\/4: 
V2,W2 

< E E + E E e-^- a ^e-^- a ^ 4 (1.29) 
^ a ez|\A £ ueA| aeA^ ue(z|)*\A| 

with a positive constant //. This right-hand side is easily shown to vanish as £ j 00. ■ 
Using the identity, X) a ez 2 [-Di2(v, w, u; a) — P 2 i(v, w, u; a)] = — (v — w) x (w — u), one 



has 



2ni 

l!(P F -,n,£p) = -— X X(v-w)x(w-u)Tr X£ (u)P FX£ (v)P FX£ (w)P FX£ (u) 

^ uGAJ v,w 

= T(Pr,n,£ P ), (1.30) 
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where we have used the expression (9.23) of T e (P F ;{l,£ P ) and (9.24). Combining this, 
(9.28) and Lemma 1.4, we obtain Theorem LI. 

Acknowledgements: I would like to thank Shu Nakamura, Hermann Schulz-Baldes and 
Hal Tasaki for helpful discussions. 
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